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Résumé. — La première partie de cet article est consacrée à l’étude théorique des variations de 
densité électronique produites par un potentiel perturbateur central dans un gaz d’électrons libres. 
Elle met en évidence le fait que, même si la portée du potentiel est courte, les variations de densité 
peuvent permettre des interactions à distance assez grande dans un métal, car leur amplitude 
décroît seulement comme l’inverse du carré de la distance au centre diffuseur. La charge d’écran 
est presque toute concentrée dans la région où le potentiel est notable ; à l'extérieur de cette région, 
les variations locales de densité peuvent être importantes même si leur intégrale dans l’espace est 
petite. On justifie enfin, en dégageant sa signification, l’emploi d’un potentiel en puits carré pour 
représenter un atome dissous dans un métal, sous réserve de satisfaire à la règle des déphasages. 

Ces résultats sont appliqués dans la seconde partie pour expliquer les variations du déplacement 


de Knight de la matrice dans un alliage dilué. La théorie développée rend compte de façon quanti- 


tative de la proportionnalité de ces variations aux concentrations d’alliage, ainsi que de l’élargis- 
sement des raies de résonance. 


Abstract. — The first part of this paper is devoted to a theoretical study of the changes of 
electron density in a free-electron gas due to a central perturbing potential. It shows that, 
even for short-range potentials, the density variations allow for fairly long range interactions in a 
metal, for their amplitude decreases only as the inverse square of the distance from the diffusing 
center. Almost all thescreeningcharge is concentrated in the region where the potentialis important, 
outside this region, the local variations of density may be important even if their integral over the 
space is small. Finally, the use of a square-well potential to represent an atom dissolved in a 
metal, as long as it obeys the phase-shift sum-rule, is justified and its meaning is explained. 

These results are used in the second part of the paper to explain the variations of the Knight 


- shift of the matrix in a dilute alloy. The theory accounts quantitatively for these variations being 


proportional to the concentration of the alloy, as well as for the broadening of the resonance lines. 


PREMIÈRE PARTIE 


PERTURBATION D'UN GAZ D’ÉLECTRONS 
PAR UN POTENTIEL A SYMÉTRIE SPHÉRIQUE 


- À. Introduetion. — Les phénomènes que nous 
étudions dans la suite sont caractérisés par la 
situation physique suivante : dans un métal À,ona 
introduit soit des particules chargées ponctuelles 
(positon, muon, trou positif dans une couche 
couche interne), soit des atomes d’un autre métal B 
qui forme avec À un alliage de substitution. Que ce 
soit une particule ou un atome B dissous, l’impu- 
reté introduite perturbe la structure électronique 
du métal À. Nous nous proposons de dégager cer- 
tains aspects de cette perturbation qui sont essen- 


tiels pour l'interprétation de diverses expériences : 
déplacement de Knight dans les alliages, annihi- 
lation du positon dans les métaux, etc. À cet 
effet, nous commençons par schématiser la situa- 
tion à l’extrême de la façon suivante : nous assi- 
milons les électrons de conductibilité du métal À à 
un gaz d'électrons libres pour lequel limpureté 
joue le rôle d’un potentiel diffuseur. Nous suppo- 
sons que ce potentiel a la symétrie sphérique et 
nous étudions les variations de densité qu'il crée 
dans le gaz d'électrons. Nous rappelons d’abord 
brièvement la méthode de Thomas-Fermi linéarisée 
pour montrer qu’elle est insuffisante, puis nous expo- 
sons plus en détail une solution meilleure basée sur 
la méthode des ondes partielles ; enfin, nous ter- 
minons cette première partie en indiquant les 
raisons qui peuvent déterminer le choix d’un poten- 
tiel simple selon le problème particulier considéré. 
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2. Méthode de Thomas-Fermi linéarisée [1], [2], 
[3]. — Nous supposons que le potentiel pertur- 
bateur Vp(r) ne dépend que de la distance 7 au 
centre Ô de l’impureté. En son absence, et si on 
néglige les corrélations, la densité du gaz d’élec- 
trons est uniforme. Compte tenu des deux direc- 
tions possibles du spin de chaque électron, la den- 
sité des états quantifiés du gaz non perturbé est [4]: 
g(E) — kr? par unité de volume et unité d’éner- 
gie (!) pour les électrons d’énergie comprise dans 


1 ; 
l'intervalle (E, E + dE), avec E = 3h. SOIt 


1 , : Ne ; 
Er — = kè l’énergie de Fermi du gaz non perturbé ; 


2 
sa densité par unité de volume vaut : 
: kg rer À.  {) 
a\En) = [side FE 


Dans l’approximation de Thomas-Fermi, on 
admet que si V> varie peu sur une distance compa- 
rable à la longueur d’onde des électrons au niveau 
de Fermi, la densité électronique vaut, à la dis- 
tance r dans le gaz perturbé : 


q'\Er, r) = qlEr — V(r)] 


Ex désignant le niveau de Fermi du gaz perturbé. 
On suppose qu’un fond continu de charge positive 
de densité uniforme qg(£%) assure la neutralité élec- 
trique du système pris dans son ensemble. L’équa- 
tion de Poisson à laquelle satisfait le potentiel V» 
s’écrit : 

AVp = — 4n(q(E; — Ver) — q(Er)] (2) 


où À est l’opérateur laplacien. 

Pour les valeurs de r telles que V» soit petit 
devant Er, on obtient, en développant le second 
membre de l’équation (2) au premier ordre en PV, : 


AVp = X(Vp — Er) (3) 
avec : 


M. __ jp re ùq __ kr 
dEp — Er— Er et: = dr (SE) nn, = 


Appliquons ceci au cas où le potentiel V, est dû 
à l'introduction d’une charge ponctuelle Z en 0. 
Les électrons du gaz se réarrangent autour de cette 
charge de façon à lui faire écran à grande distance. 
Les seules solutions acceptables de l’équation (3) 
sont celles qui s’annulent à l'infini. 


(2) Sauf indication du contraire, nous utilisons dans tout 
le texte les unités atomiques (U. A.) de Hartree : charge de 
l’électron : [e| — 1 ; masse de l’électron au repos : m — 1 ; 
constante de Planck réduite : 4 = 1. 

En particulier, l’unité de longueur est égale au rayon de 
Bohr de l’atome d’hydrogène : as — 0,529.10—8 em ; 
l’unité d’énergie vaut deux fois l’énergie d’ionisation de 
l'hydrogène, soit 27,2 eV. 
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Elles sont de la forme : 


e—Àr 
Vr = const. X A avec ÔEp = 0 (4) 
indiquant que le niveau de Fermi du gaz reste 
inchangé. Au voisinage de l’origine, les hypothèses 
de l’approximation de Thomas-Fermi ne sont plus 
vérifiées ; l'écran y devient négligeable et V» tend 
vers — Z fr. Tenant compte de la forme (4) de la 
solution à grande distance, on prend : 


Pr D (5) 


dans tout l’espace. L’équation (2) donne alors pour 
accroissement de la densité de charge électrique par 
unité de volume, à la distance r de 0, la valeur : 


e—Àr = LE CAT (6) 


En intégrant dans tout l’espace, on vérifie que la 


charge d’écran totale | 3g dv est bien égale à Z. 


La formule (6) donne une intensité infinie à l’ori- 
gine, même pour les électrons d’énergies comprises 
dans un intervalle Æ, £ + dÆ très étroit. Ceci est 
en contradiction avec les mesures du déplacement 
de Knight de la résonance nucléaire [51 et celles du 
temps de vie du positon dans les métaux [6], [7]. 
De plus, d’après cette formule, l'excès de densité 
électronique décroît exponentiellement avec la dis- 
tance au point 0, en gardant un signe constant, 
celui de Z. Avec Z positif et égal à la différence de 
valence entre l’argent d’une part, et le cadmium, 
Pindium, l’étain ou l’antimoine, d’autre part, on 
obtient un résultat du signe contraire à celui qui 
convient pour expliquer les variations observées du 
déplacement de Knight de l’argent quand on y 
dissout l’un de ces métaux. 

En résumé, la méthode de Thomas-Fermi ne 
donne qu’une image assez grossière de la polari- 
sation d’un gaz d'électrons par un centre chargé. 
Elle fournit, au moyen du paramètre à, un ordre de 
grandeur du rayon de l'écran électronique : 
quelques unités atomiques. Inutilisable au voisi- 
nage du centre, elle ne révèle pas les détails de la 
variation de densité électronique avec la distance. 


3. Méthode des ondes partielles. — Pour faire 
mieux que l’approximation de Thomas-Fermi, il 
devient nécessaire de résoudre l’équation de Schrô- 
dinger pour le potentiel perturbateur. Ceci se fait 
avantageusement par la méthode des ondes par- 
tielles quand le potentiel a la symétrie sphérique. 
Nous obtiendrons ainsi une densité électronique 
finie en tout point et, en particulier à l’origine, et 
des variations locales de densité s’étendant bien 
au delà du potentiel pour les électrons d’un petit 
Intervalle d'énergie. 

Considérons un gaz d’électrons libres dont nous 
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négligeons les corrélations. Au point de rayon vec- 
teur r, la fonction d’onde d£ d’un électron d’éner- 


gie E — 5 k? et d’impulsion k, normalisée dans un 


grand volume Y, s'écrit : 


L — LE ik.r 
ÿk(r) Free (2) 

Introduisons dans ce gaz un potentiel pertur- 
bateur Vr(r) à symétrie sphérique autour de l'ori- 
gine et nul pour r > a, le rayon a étant, pour fixer 
les idées, de l’ordre de quelques unités atomiques, 
comme le suggèrent les résultats de la méthode de 
Thomas-Fermi. 

Nous allons montrer que la diffusion des élec- 
trons par ce potentiel crée, même à grande distance, 
des fluctuations dans la densité du gaz d’électrons. 

SOL : 


We = D loir + Jf] (8) 


la fonction d’onde d’un électron d’impulsion k dif- 
fusé par le potentiel Vr. Cette fonction a la même 
amplitude que dx à grande distance ; en vertu du 
théorème de von Laue [8], nous conservons la 
même densité d'états g — k/r? par unité de volume 
et unité d’énergie que dans le gaz non perturbé. 
Dans l’espace réciproque, ceci correspond à une 
densité d'états & par unité de volume telle que : 


ornnade ad redR Ir esoite ot AT 


En l’absence du potentiel Vr, la densité au 
point r des électrons d’énergie comprise entre Æ£ 
et E + dE était : 


dl 
LT 


4 
p(E) dE = V ie dk [7 u* Vf dQ = g dE 
(dQ est l’élément d’angle solide dans l’espace réci- 
proque). En présence de Vr, elle est devenue : 


k2 dk jar 
Re Een 
d£ 4T : 

Ge D da. 


L'introduction du potentiel Vr produit donc au 
point r une variation relative de densité des élec- 
trons d'énergie E égale à : 

SE, r) __ V 


b(EL en [Uk Uk — dk” UK] dQ. (9) 


(LE) + So(E, r)] dE = 


— 4} 


La symétrie du potentiel Vr nous permet d’ana- 
lyser les fonctions d’ondes (7) et (8) en composantes 
sphériques. Soit 0 l’angle des vecteurs k et r. 
Pour r > a,ona[9]: 
eik.er — D ü V2 + 1) jifkr) PA (cos 0) 

1=0 
(10) 
eiker + Jf — Z üV2(21 + 1) eëm fi(kr) P, (cos 6) 
1 


=0 / 


EFFETS DES IMPURETÉS DANS LES MÉTAUX 


IGN 
avec : 
Atkr) = jifkr) cos vw — nfkr) sin y. 


Dans ces formules, 7,et r,sont respectivement les 
fonctions sphériques de Bessel et de Neumann 
d'ordre entier {, P,le polynôme de Legendre normé 
d’ordre let n, le déphasage introduit par le poten- 
tiel Vr dans la le composante ‘sphérique de la. 
fonction d'onde. 

En portant les développements (10) dans l’éga- 
lité (9), il vient, après intégration sur toutes les 
directions de k : 


In£ (kr) — jE(kr)] sin? 9 — jifkr) m(kr) sin 27 L. (11) 


Quand r croît, cette expression prend des valeurs 
alternativement positives et négatives. Pour les 
grandes valeurs de r, amplitude de ses oscillations 
décroit comme 1/r?; on s’en rend compte clai- 
rement sur la forme asymptotique [3] : 


Ô 1 à 
RASE Ge ape (2e + 1) sin UE 


quand r — co. 

Quand on connaît les déphasages », en fonction 
de k, on peut calculer, à l’aide de la relation (11), 
la quantité Q de charge d’écran qui se trouve à 
l'extérieur de la sphère de rayon a : 


APTE 
Q _ dE [_ Ge(E, r) dv 
= D | (21 + 1) Es k2 dk  {trétkr) — ji (kr)] 


sin? (x) — jifkr) mtkr) sin 2m(k) {r? dr 


de ED 1! re = Loy(ka) sin? my + Bulka) sin 2n] (12) 
avec [10] 
God) ET 1. [nè(u) — jé(u)] u? du 
= — 97 sin 2x 
Bolx) = — 47 1 Jolu) nou) u? du = rx cos 2x 
au(x) — 4r(2l + 1) de [nf(u) — jJ£(u)] u? du 
= Qn(21 + 1) a°[jf(x) — f(x) \ (13) 


Es J1—1(x) j1+1(&) n m1 (x) ny +1(x)] 
Bale) = —ar (20 + 1) JO duo mu ut qu 
1x0 


= In(2 + 1) x) ji(x) mx) 


5 Dale) ml 4 dr ae mu | 
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Nous allons montrer sur un exemple que les 
variations locales de densité dans le gaz d’électrons 
peuvent être importantes bien que la charge 
d’écran soit presque entièrement localisée dans la 
sphère d’impureté, c’est-à-dire bien que la charge Q 
extérieure à cette sphère soit très petite. 


4. Distribution de charge dans le cas d’un puits 
carré de potentiel. — Nous montrerons au para- 
graphe suivant que la perturbation produite dans 
un métal par une impureté de substitution peut être 
représentée avec une bonne approximation par un 
puits de potentiel de rayon et de profondeur conve- 
nablement choisis. Puisque cette étude nous ser- 
vira par la suite, nous en profitons pour illustrer les 
résultats du paragraphe précédent sur l’exemple 
simple où Vr est un puits de potentiel de rayon 4 


F1G. 1. — Potentiel en puits carré. 


et de profondeur k$/2 (fig. 1). Nous traiterons plus 
loin, à propos de l’annihilation du positon, le cas du 
potentiel coulombien avec écran. 

Les déphasages se calculent aisément sur ce 
modèle [11]. En posant : #& + 4? = k'?, leurs tan- 
gentes s'expriment au moyen des fonctions de 
Bessel et de Neumann sphériques par les formules : 


mA ka j1(ka) — Bo j0(ka) 
18 0 — Fa na(ka) — Po nolKe) 
, Jatk'a) 
ec = k'a- 
SR) 
__ ka j—(ka) — y jifka) 
BU Ko ny ilta) — Pirate) 


’ ji—1(k'a) 
avec : = ka —— 
: n'a) 


? 


(14) 


pour L > 0. ) 


La figure 2 représente les déphasages que nous 
avons calculés de cette façon pour le puits de poten- 
tiel caractérisé par : k, a = 1,4 et kr à = 2, avec 
kr — 0,64. Au delà de ! — 2, les déphasages sont 
négligeables. Ce puits déplace une charge d’écran 
d'environ 0,9 U. A. ; il n’introduit pas d’état 
lié. Nous verrons au chapitre suivant qu'il cor- 
respond approximativement, compte tenu de 
l'effet de taille, à un atome de cadmium dissous 
dans l’argent. La figure 3 montre que la charge 
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d’écran ne déborde pratiquement pas du puits car 
l'aire totale comprise sous la courbe jusqu’au 
niveau de Fermi est nulle à la précision des calculs 
soit à deux ou trois centièmes d’unités de charge 


près. 
ÿ 
1 


0 L 2-k@ 


F1c. 2. — Déphasages introduits par le potentiel 
en puits carré caractérisé par k.a = 1,4. 


F1c. 3.— Répartition de la charge d’écran Q extérieure 
au puits de potentiel. 


À l’intérieur du puits, c’est-à-dire pour r < 4, 
chaque fonction j, de la formule (10) est remplacée 
par À, j{(4'r). Les constantes À, sont déterminées 
par les conditions aux limites sur la sphère r — 4: 


A ji{k'a) = fi(ka). (15) 


Les déphasages autres que no, 1 et n° étant 
pratiquement nuls, on a, à l’intérieur de cette 
sphère : 


Ôo(E, r F 
ee. PE 2 (27 + 1) A? jf(kr) — 1. 

. La figure 4 montre comment Sp/o varie en fonc- 
tion de r au niveau de Fermi. La densité de l’écran 
à l’origine est finie, contrairement à ce que donnait 
la méthode de Thomas-Fermi. Pour 4er — 3,5, ce 
qui correspond sensiblement à la distance entre 
noyaux voisins dans l’argent, on observe une dimi- 


(16) 


à 


e 
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” nution de la densité des électrons au niveau de 


Fermi. Nous verrons que ceci explique la dimi- 
nution du déplacement de Knight de l'argent 
quand on y dissout du cadmium. 


Les électrons contribuent à la formation de 
3 . : 
l'écran dans un rapport variable suivant leur éner- 


F1G. 4. — Variation relative de densité électronique. 
au niveau de Fermi. 


gie et la distance au centre diffuseur. À l’origine, 
seules les composantes sphériques de type s ont une 
amplitude non nulle ; l'égalité (16) s’y réduit sim- 
plement à : 
Se(E, 0) 
p(E) 


Les équations (14) et les identités : 


— Ag(k') — 1. 


j)x) = sin xfx not) = — cos x/x 


permettent de lui donner la forme simple suivante : 


(17) 


La figure 5 représente les variations de Sp/o à 
l’origine en fonction de k. On constate que les élec- 
trons de faible énergie contribuent à former l’écran 
d’une façon relativement très forte par rapport à 
ceux qui sont proches du niveau de Fermi. En ce 


qui concerne les contributions absolues, les pro- 


portions sont modifiées du fait que la densité des 
états par unité de k : 


desk? 
8Gk 7? 


croît de zéro pour les électrons d’énergie nulle 
à 2E-/r? au niveau de Fermi. La figure 6 montre 
qu’à une distance du centre égale à une unité ato- 
mique la situation est très voisine de ce qu'elle est à 
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l'origine. Ceci tient à deux faits : d’une part les 
composantes s des fonctions d’onde varient très peu 
pour les petites valeurs de r, puisqu'elles ont la 


sin X'r 
forme ; d'autre part les déphasages d’ordre 


do (0) 


e 


Lt 


20 


Poe k 


F 


F1G. 5. — Accroissement relatif de densité 
au centre du puits carré. 


1 > 0 ne commencent à devenir importants qu’au 
voisinage du niveau de Fermi, comme le montre la 
figure 2. 


Ô@(r=1) 


20 


10 


RER 


Fic. 6. — Variation relative de densité. 
en un point intérieur au puits, près du centre (r = 1). 


En un point extérieur à la sphère du rayon 4, 
l'étude de la structure de l’écran se fait au moyen 
de la formule (11) dans laquelle on fait varier k en 
maintenant r constant ; dans cette formule, les 
déphasages dépendent de # comme l'indique la 
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figure 2. La figure 7 montre ce que l’on obtient à la 
distance r définie par Arr — 3,5. Les électrons de 
faible énergie sont attirés par le puits de poten- 
tiel (8o/e > 0), tandis qu’on observe une légère 


Do(r-35) 
0 


0 SN PER 

RE À 

Fic. 7. — Variation relative de densité 
en un point extérieur au puits (7 — 3,5). 


répulsion (30/9 << 0) des électrons d'énergie voisine 
du niveau de Fermi. Ce dernier résultat traduit 
directement la diminution de densité des électrons 
au niveau de Fermi à cette distance, que montre la 
figure 4. 


o. Choix du potentiel diffuseur. — Les servitudes 
du calcul obligent à choisir pour Vr» la forme de 
potentiel la plus simple qui satisfasse aux condi- 
tions physiques du problème particulier que l’on 
veut traiter. 


A) CONSIDÉRATIONS GÉNÉRALES. — Supposons 
que les électrons de conductibilité du métal À pur 
se comportent comme des électrons libres dont la 
fonction d’onde soit simplement modulée par celle 
du bas de la bande de conductibilité. Autrement dit 
admettons que la surface de Fermi soit sphérique 
et que la masse effective des électrons soit égale à 
leur masse normale. Dans ces conditions, la fonc- 
tion de Bloch d’un électron d’énergie 


FRERE : k? 
dans le métal À s'écrit : 


il 
PART rietkire 18 
) y a(r) (18) 
La fonction uA(r) à la périodicité du réseau cris- 


tallin du métal À ; elle est solution de l’équation 
de Schrôdinger : 


Aus + 2(E 4 — V4) ua = 0 (19) 


pour le potentiel du réseau VA et l'énergie £4 d’un 
électron au bas de la bande de conductibilité dans 


à 
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le métal À. La normalisation de YA dans le 
volume V du métal implique que w4 soit norma- 
lisée dans l’unité de volume. 

Si on pose : 


PAlk, r) — war) Yk(r), 
l'équation de Schrôdinger : 
AP 4 + QE — Va) Ya = 0 


montre, compte tenu de l’équation (19) pour wa, 
que 4 doit satisfaire à : 


ua[A UK + k2 Jk] + 2Au4.Aÿk = 0. 
Prendre pour WA la fonction (18), c’est-à-dire 
dk == he eiker, 
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revient à négliger le terme Aua.Adx dans l’équa- 
tion (20). Nous discutons cette approximation dans 
l'annexe I ; elle se justifie dans une certaine mesure 
car, dans les métaux purs, les fonctions d’ondes de 
bas de bande ne présentent d’oscillations que dans 
un volume limité autour des noyaux, tandis qu’elles 
sont très constantes au delà, en raison des condi- 
tions imposées à la surface des polyèdres atomiques 
par la périodicité du réseau [1], [12]. De plus, 
l'amplitude de la fonction de bas de bande, dans la 
région où elle est pratiquement constante, ne dif- 
fère pas beaucoup de celle qu’aurait un électron 
libre uniformément réparti dans le volume ato- 
mique. Ces remarques vont nous guider dans la 
détermination du potentiel diffuseur. 


(20) 


B) CAS D’UNE IMPURETÉ DE SUBSTITUTION. — 
Considérons un premier type de problème, celui où 
l’on étudie les modifications de densité électro- 
tronique produites dans un métal À monovalent 
dont on a remplacé un, atome par un atome de” 
métal B de valence Z + 1. 

On sait que, dans un alliage AB à faible concen- 
tration de métal B, l’écran électronique est très 
concentré autour des atomes B dissous et ne s’étend 
guère au delà des polyèdres atomiques corres- 
pondants [1], [2]. Le calcul de la quantité Q au 
paragraphe précédent le montre d’ailleurs sur un 
modèle particulièrement simple. La forme de la 
fonction de bas de bande dans la matrice n’est donc 
pas beaucoup modifiée par la présence du soluté. 

Considérons alors un réseau analogue à celui de 
la matrice, mais dont le potentiel périodique serait 
engendré par la répétition de celui que subissent 
les électrons de conductibilité de l’alliage dans le 
polyèdre atomique de l’atome B dissous. Dans ce 
potentiel périodique V3, un électron aurait, au bas 
de la bande de conductibilité, l'énergie Æ3 et la 
fonction périodique w8 solution de l'équation de 
Schrüdinger correspondante : 


AuB + 2(E8 = V3) UB —= 0, (21) 


de a EU 


VE Part 
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Les conditions aux limites imposées par la pério- 


_dicité du réseau obligent la fonction uy à avoir un 


gradient nul sur la surface du polyèdre atomique. 
Comme ua dans la matrice, elle est donc assez 


_ constante dans la majeure partie du volume ato- 


mique de B ; ses oscillations ne deviennent impor- 
tantes qu’au voisinage immédiat du noyau B. 
Considérons maintenant l’alliage formé d’un 
atome B dissous dans la matrice À. On sait qu’à 
dilution infinie le niveau de Fermi de l’alliage reste 
celui du métal À pur [2]. Soit Ÿ la fonction d’onde 


d'un électron d’énergie Æ — E\ +5 dans 


l’alliage. Suivant la méthode de Wigner-Seitz, 
remplaçons le polyèdre atomique de l’atome B 


dissous par une sphère de même volume. Soit a le 
rayon de cette sphère. Posons : 


PNG 


l’origine des coordonnées étant prise au noyau B. 
D’après ce que nous avons dit plus haut, le poten- 
tiel perturbateur dû à l’atome B est pratiquement 
nul pour r > a. Par conséquent, dans les conditions 
de validité de l'expression (18), la fonction dx satis- 
sa à l’équation de Schrodinger pour les électrons 
hbres : 


Axe pour (22) 


A%k + k? Üx = 0 pour r > a. (23) 


Posons maintenant, à l’intérieur de la sphère 
atomique de B : 


P = up Ÿk (r < a). (24) 


La fonction ug étant solution de l’équation (21), 
nous trouvons que dans cette sphère dx doit vérifier 
l'équation suivante : 


(25) 
up $ AUk + [k? + 2(F4 — Eg)] ÿx } + 2A%k. Aug = 0. 


Ce que nous avons vu précédemment au sujet 
des fonctions de bas de bande des métaux purs 
indique que vu est pratiquement nul quand la 
distance r au noyau dépasse une unité atomique 
environ. Dans ce domaine, nous pouvons donc 
négliger le produit vdx.Aus de l’équation (25) 
pourvu que la fonction dX ne varie pas trop vite. 
Par contre, au voisinage immédiat du noyau B, 
où la composante s de dX est la plus importante, 
c’est le facteur y dx qui s’annule. En effet, pour les 
petites valeurs de r, l'équation (24) admet, si l’on 
néglige le produit ÿ Ÿx Vus, la solution : 


de = sin VR? + 2(E4 — Eg).r 
— RUB DL AT 0 0 ? 
régulière à l’origine et pour laquelle y®x(0) — 0. 
Nous pouvons donc, avec une bonne approximation 
prendre pour ÿ les solutions de l’équation : 


Aÿk + (X-+ 3) Jk 0 pour r<a, (26) 
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en posant : 
2(E 4 TT F3) = He 


Les équations (23) et (26) montrent que la fonc- 
tion d’onde d4 est celle d’un électron libre d’éner- 


5 : 
ie LA et d’impulsion k, diffusé par un puits de 


potentiel de rayon a et de profondeur ; Ho eMhne 


Nous discutons cette approximation plus en détail 
dans l’annexe I. 

Les fonctions d’onde u1 et un normalisées dans 
l’unité de volume n’ont pas, en général, strictement 
la même valeur pour r = a. La continuité de la 
fonction d’onde totale sur la sphère atomique : 


uala) Yk(a + 0) = ugla) Ykla — 0) 


entraîne alors une discontinuité de la fonction x. 
La continuité en a de la dérivée de la fonction 
d'onde totale impose la même discontinuité pour 
les dérivées de dx de part et d’autre de la surface 
de la sphère : 


uala) ES cit cales 


car : 


du1 dug 
Es de (GE). os 

Par conséquent, la dérivée logarithmique de dx 
est continue en a. Dans ces conditions, les dépha- 
sages créés par le puits de potentiel dans les compo- 
santes sphériques de 4 sont les mêmes que pour 
une fonction d’onde d’électron libre de même éner- 
gie satisfaisant aux conditions habituelles de conti- 
nuité pour r — a. La fonction d£ est simplement 
ua(a) 
ug(a) 
teur étant très voisin de l’unité, nous ne tiendrons 
pas compte de cette petite correction. 

La détermination précise de l’énergie de bas de 
bande £A1 d’un métal est très difficile ; quant au 
calcul direct de Æ3, il semble à peu près impossible. 
Cependant, on peut déterminer 4,. de la façon 
simple suivante. En effet, on sait que si atome B 
dissous dans le métal monovalent À a la valence 
Z +1, les électrons de conductibilité de l’alliage 


ainsi formé se polarisent autour de B de façon à 
faire écran à sa charge positive excédentaire Z. 


multipliée par le facteur pour r < a. Ce fac- 


1 
Quand le puits de potentiel a la profondeur 5 kè 


convenable, les déphasages n, qu'il introduit dans 
les composantes sphériques de la fonction d’onde 
d’un électron au niveau de Fermi satisfont à la 
condition de Friedel [2] : 


(27) 


ain 


RE D (27 + 1) kr). 
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Cette condition assure que la charge d’écran 
déplacée dans le gaz d’électrons par la pertur- 
bation est bien égale à Z, comme l'exige l’équilibre 
électrostatique du métal. En calculant les dépha- 
sages au niveau de Fermi à l’aide des formules (14), 
il est donc facile de déterminer k, par approxi- 
mations successives. 

Lorsque les métaux À et B ont des volumes ato- 
miques différents, on peut en tenir compte suivant 
le procédé utilisé avec succès par Blatt pour cal- 
culer la résistivité résiduelle du cuivre ou de 
l'argent contenant de petites quantités d’atomes 
étrangers dissous [13]. La méthode de Blatt aboutit 
à remplacer dans la formule (27) la charge d'écran Z 
par une charge effective Z’ plus grande ou plus 
petite que Z selon que l’atome B est plus petit ou 
plus gros qu’un atome À. 


C) CAS D’UNE IMPURETÉ PONCTUELLE EN POSI- 
TION INTERSTITIELLE. — Dans le cas où l’on étudie 
la densité électronique sur une particule ponctuelle 
de charge positive (positon ou muon par exemple), 
introduite dans un métal, il semble raisonnable de 
prendre un potentiel coulombien avec écran de la 
forme : 

Z 


Vp = — e—Àr 


; (28) 


si Z est la charge de la particule. En effet, au voisi- 
nage de l’origine, où l'effet d’écran devient négli- 
geable, ce potentiel se comporte correctement 
en Z/r. Aux distance plus grandes, il a la forme 
requise par l’approximation de Thomas-Fermi. 
Cependant, le paramètre d’écran À doit être déter- 
miné de façon que les déphasages au niveau de 
Fermi, calculés en résolvant l'équation de Schrü- 
dinger pour ce potentiel, vérifient la relation de 
Friedel (27). | 

Si on prend des fonctions d’ondes non perturbées 
de la forme (18), les fonctions d’ondes perturbées 
correspondantes s’écriront dans ce cas : 


FR, 7) = u4 xl), 


l'équation d’onde pour les fonctions dx étant : 


ua [Au + (+ Lean qu] + 2qua.vuk = 0. (29) 

Chaque fonction x peut être considérée comme 
la superposition d’une onde plane incidente et d’une 
ondelette diffusée, conformément à la formule (8), 
si on néglige le terme wua y4f dans le métal pur. 
Dans ces conditions, le terme Aua Adg qui appa- 
rait dans l’équation (29) est aussi négligeable. En 
effet, l'amplitude de l’ondelette diffusée 49 n’est 
importante que pour r inférieur à 4 /X, c’est-à-dire 
de l’ordre de l'unité atomique, donc dans une région 
où vu est pratiquement nul. La raison en est que 
les particules chargées positivement sont repous- 
sées par les noyaux ; elle est donc tendance à se 
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localiser dans les régions interstielles où les fonc- 
tions d’onde des électrons de conductibilité sont 
assez proches de celles d’électrons libres. Dans ces 
conditions, l'équation (29) pour x se réduit à celle 
de la diffusion d’électrons libres pour un potentiel 
coulombien avec écran : 


Aÿk + (ee de 2 ex) Uk = 0. 


6. Conclusion. — L'étude théorique qui précède 
nous à donné les résultats suivants : re 
19 La perturbation d’un gaz d’électrons libres 


par un potentiel central à portée, même très courte, : 


produit dans le gaz des variations de densité d’allure 
oscillante en fonction de la distance au centre dif- 
fuseur. Pour les électrons d’un intervalle étroit 
d'énergie, l'amplitude de ces variations décroît seu- 


lement comme l’inverse du carré de la distance. Par 


conséquent, la perturbation des électrons de con- 
ductibilité d’un métal pourra se manifester à grande 
distance par des effets sensibles aux variations 
locales de densité. 


20 La charge d’écran est presque toute concen- 


trée dans le volume où le potentiel est important. 
À l'extérieur de celui-ci, l'intégrale des variations 
de charge est très petite ; mais les changements 
locaux de densité électronique peuvent ne pas être 
de moyenne nulle si on les considère seulement en 


des points particuliers, tels que les noyaux des. 


atomes dans un cristal, comme c’est le cas en réso- 
nance nucléaire. 

30 On pourra, avec une bonne approximation, 
traiter les électrons de conductibilité d’un métal 
comme des électrons libres en prenant un potentiel 


diffuseur de la forme (28) si la perturbation est pro- | 


duite par une impureté ponctuelle chargée positi- 
vement : positon, muon. S'il s’agit d’une impureté 
de substitution telle qu’un atome dissous, formant 
un alliage infiniment dilué, il sera meilleur d’intro- 
duire un puits de potentiel de rayon égal au rayon 
de l’atome dissous et de profondeur telle que la 
condition d’écran soit satisfaite. 

Dans la deuxième partie de cet article et dans. 
l’article suivant nous appliquons ces résultats au 
traitement de problèmes physiques réels. 


DEUXIÈME PARTIE 


DÉPIACEMENT DE KNIGHT 
DANS LES ALLIAGES 


1. Introduetion. — Dans un métal pur, la réso- 
nance magnétique nucléaire d’un isotope est dé- 
placée vers les fréquences élevées par rapport à 
celle qu'on observe pour le même métal dans un 
sel [5] (1): Le déplacement relatif de fréquence : 


1) Récemment, Rowland a observé un effet Knight 
négatif sur le platine [14]. 


(0) 


on dE TT EN A OS TT UE er PI EME 0 OT PU OT 
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K = AyJvest de l’ordre de 0,1 à 4 %,. Il est dû au 
paramagnétisme de Pauli des électrons de conduc- 
tibilité [15]. Au niveau de Fermi du métal, ces 
électrons s’orientent dans le champ magnétique 
statique À appliqué à l’échantillon ; l'interaction 
de contact de leurs moments magnétiques avec ceux 
des noyaux équivaut à un accroissement local AH 
du champ magnétique tel que : 
AH 87% 


RE y PO >. 


H 53 (31) 


Dans cette formule, x désigne la susceptibilité 
paramagnétique de Pauli par unité de volume du 
métal ; on sait qu’elle est proportionnelle à la den- 
sité g(Er) des états au niveau de Fermi, par unité 
de volume et unité d’énergie [4]; < |[T#(0)? > 
est la valeur moyenne, par état, de la densité élec- 
ironique au noyau pour les états d’un petit inter- 
valle d'énergie à la surface de Fermi, les fonctions 
d'onde W} étant normalisées dans le volume V du 
métal. | 

Dans un alliage AB, la densité des états au 
niveau de Fermi varie avec la concentration. De 
même, les fonctions d’onde Ÿ; se modifient tant 
aux noyaux À qu'aux noyaux B. À eux seuls, ces 
deux effets doivent suffire à modifier les dépla- 


. cements de Knight respectifs des noyaux À et B. 


Poutant, les premières expériences faites sur des 
alliages semblaient indiquer que le déplacement de 
Knight de la matrice était peu sensible à la pré- 
sence du soluté jusqu’à des concentrations pouvant 
dépasser 10 % [5], [16], [17]. Ces expériences 
n'étaient peut être pas assez précises pour révéler- 
de petites variations systématiques du déplacement 
de Knight ; ces variations devaient être masquées 
par de forts élargissements quadrupolaires dans 
des alliages tels que Cu Zn, Cu Al, AI Cu, Ai Mg 
et AI Zn dont les constituants portent des spins 
nucléaires 3/2 ou 5/2. De son côté, la résonance du 
soluté semblait aussi peu sensible à la concen- 
tration bien qu’elle eût lieu à une fréquence diffé- 
rente de celle du métal B pur. | 
L'absence de variation du déplacement de 


Knight de la matrice pouvait s'expliquer dans, 


l’approximation de Thomas-Fermi étendue au cas 


d’une concentration non nulle d’atomes dissous [2], : 


[18]. On peut d’ailleurs s’en rendre compte au 
moyen de la formule (6). : 
Cependant, les mesures systématiques de Drain 
sur les alliages argent-cedmium révèlent une dimi- 
nution du déplacement de Knight de l'argent, 
linéaire en fonction de la concentration en cad- 
mium [19]. De plus, le déplacement de Knight du 
cadmium dissous diminue lui aussi, suivant la même 


loi. La linéarité se conserve pratiquement dans tout. 


le domaine de stabilité de la phase & (C. F. C.) à la 
température ordinaire, soit jusqu'à plus de 40 A 
de cadmium dissous. Ces résultats sont confirmés et 
étendus, .en ce qui concerne l’argent, par les mesures 


. 
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de Rowland sur les alliages Ag Au, Ag Cd, Ag In, 
Bg Sn et. Ag Sb [20]. Enfin, Bloembergen et Rimaï 
dans des alliages alcalins liquides [21], Oriani et 
Webb [22] dans Na Hg et Na Au observent éga- 
lement une variation du déplacement de Knight 
du solvant proportionnelle à la concentration tant 
que celle-ci reste assez faible. 

La raie de résonance du solvant s’élargit consi- 
dérablement dans les alliages, mais l’élargissement 
croit nettement moins vite que la concentration. 
Dans la série des alliages d’argent, la variation du 
déplacement de Knight et l'élargissement de la raie 


de l’argent augmentent quand on passe de l’or à 


Pantimoine, c’est-à-dire à mesure que la valence du 
soluté croît. L'intérêt de l’argent vient de ce que 
ses isotopes magnétiques ont un spin nucléaire 1/2, 
donc née présentent pas de moment quadrupolaire. 

Ces effets ne s'expliquent pas dans l’approxi- 
mation de Thomas-Fermi. Ils sont beaucoup trop : 
forts pour être dus simplement à la distorsion de la 
matrice par des atomes dissous de volume différent. 
(Cf. annexe IT.) Ils varient de façon trop régu- 
lière avec la concentration pour être attribués à des 
changements de caractère s-p des fonctions d’onde 
au voisinage d’une limite de zone de Brillouin [23]. 
Nous en proposons au contraire l’explication sui- 
vante : la variation du déplacement de Knight 
avec la concentration et l’élargissement de la réso- 
nance sont des conséquences des fluctuations de 
densité électroniques dues à la polarisation des élec- 
trons de conductibilité de l’alliage par les atomes 
dissous [24], [251. 

Dans la première partie, nous avons montré 
l'existence de ces fluctuations dans un gaz d’élec- 
trons libres perturbé. Nous allons maintenant déve- 
lopper le calcul qui doit justifier notre interpré- 
tation. Nous commencerons par le cas d’un seul 
atome dissous dans un métal monovalent, puis nous 
passerons au cas réel d’une concentration faible 
mais non nulle de soluté. Après avoir calculé la 
variation du déplacement de Knight et l’élargis- 
sement des raies pour une solution solide, nous 
étendrons la théorie aux alliages liquides. 


2. Effet d’un atome dissous sur la résonance 
nucléaire de la matrice. — Nous allons montrer 
que la présence d’une impureté polyvalente B dans 
un métal monovalent À entraîne un étalement de 
la fréquence de résonance nucléaire de la matrice. 
Pour simplifier les calculs, tout en dégageant 
l'essentiel du phénomène, nous supposons que les 
électrons du niveau de Fermi se comportent dans 
le métal À pur comme des électrons libres dont la 
fonction d’onde serait modulée par celle du bas 
de la bande de conductibilité. Nous avons vu que, 
dans ces conditions, la présence de l’atome B en 
position de substitution peut se représenter par un 
puits de potentiel dont nous avons déterminé les 
caractéristiques au paragraphe 5 de la première 
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partie. Nous conservons les mêmes notations, à ceci 
près que, puisque seuls les électrons au niveau de 
Fermi de l’alliage nous intéressent pour l'étude du 
déplacement de Knight, nous omettons provisoi- 
rement l'indice F qui les caractérisait ; nous dési- 
gnons simplement par k le rayon de la sphère de 
Fermi. 

À la fonction d’onde d’un électron de vecteur 
d’onde k dans le métal À pur, 

Palr) = ualr) g£(r), (32) 
où 
— eikr, 


1 
0 — 
%k VT 


nous avons fait correspondre, dans l’alliage, la fonc- 
tion d’onde de l’électron diffusé par l’atome B : 


Pr) = ul) dx(r). (22) 


Soit & la densité des états par unité de volume 
dans l’espace des k, au voisinage de la sphère de 
Fermi. Au niveau de Fermi du métal À pur, la 
densité électronique par unité de volume et unité 
d’énergie est : 


at 
Ha kS u* u1 YU* 9 dO, 


— guË(r) ua(r), 


e(r) — 


g — 4n k © étant le nombre d’électrons par unité 
d'énergie et unité de volume au niveau de Fermi, 
et dQ l’élément d’angle solide dans l’espace réci- 
proque. 

Dans Palliage, au point correspondant de la 
matrice, la densité électronique est devenue : 


4T 
al + Bob) = 7 kan ua gt ma 0 


, guï U A ét * 
= se QE Uk AO. 


L’aceroissement relatif de densité au point r vaut 
donc : 
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Puisque ua est normalisée dans l’unité de volume, 
les fonctions dx et di de l'expression (33) sont iden- 
tiques à celles de la formule (9). Par conséquent 
nous obtenons pour Sp/p exactement le même 
résultat que dans la relation (11) : 


9 


p(r) 10) 


(21 + 1) 


ÉUne(kr) — je(kr)] sin? ne — j(kr) mu(kr) sin 2m}. (34) 

Soit R; le rayon vecteur du jè% noyau À à 
partir du noyau B. Dans le métal À pur, la den- 
sité o(R,) a la même valeur p pour tous les noyaux. 
Dans l’alliage, cette densité devient : 


e(R;) — p + p(R;) 


et varie d’un noyau à l’autre suivant la relation (33). 
D’après la formule (31), le déplacement de 
Knight Æ est proportionnel à la densité électro- 
nique au noyau: p(R,) — < [F(R;)? >. Les 
variations de p(R,) dans l’alliage font que chaque 
noyau À y a sa fréquence de résonance particulière, 
donc son propre déplacement de Knight : 
K, — K +35K,;. La variation relative ÔX;/K du 
déplacement de Knight du 7% noyau À est égale 
à celle de la densité électronique au même noyau : : 

SK;[K — Se(R;)]e. (53) 

Par conséquent, au lieu d’une fréquence de réso- 
nance unique (en supposant pour simplifier que le 
métal À ne comporte qu’un seul isotope), on 
obtient, en principe, dans l’alliage, tout un spectre 
de fréquences de résonance voisines. 

Nous avons calculé l'effet sur la résonance de 
l’argent des impuretés suivantes : Cd, In, Sn et Sb. 
Nous avons utilisé pour cela les déphasages que 
Blatt avait calculés pour l’étude de la résistivité 
résiduelle des alliages correspondants [13]. Nous 
reproduisons ces déphasages dans le tableau I, avec 
les charges effectives Z' qui tiennent compte de 


Sel) =f" (US Uk — y2* 49) dO. (33) leffet de taille, en limitant leurs valeurs aux déci- 
p(r) 47 males que nous avons effectivement utilisées. 
- TABLEAU I 


CHARGES EFFECTIVES Z’ ET DÉPHASAGES DANS L'ARGENT (D'APRÈS [13]) 


| 


ÂTOMES DISSOUS 


Z 71 TES 2% 
PÉPRASAGES NO 76 Sn et NT Re De 3.67 
0,532 0,830 1,00 111 
. 0241 0.558 0.908 1215 
“ 0.028 0,063 0109 0178 
s 0,002 0 00% 0006 0,010 


Pratiquement, comme le montre ce tableau, ils 
deviennent négligeables au delà de ! = 2 ou 3, ce 
qui limite la somme (34). 


La figure 8 illustre les résultats. On remarque 
que la polarisation des électrons par l’atome B 
affecte fortement la fréquence de résonance des 


ee 
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noyaux À qui sont ses premiers et seconds voisins. 
On voit qu’elle crée dans les couches d’atomes plus 
éloignées des fluctuations d'autant plus grandes 
que la charge Zest plus élevée : l'effet va en s’accen- 


tuant quand on passe du cadmium (Z —1) à 


Pantimoine (Z = 4). De part et d'autre d’une raie 
centrale intense, et quelque peu élargie, corres- 


ES =! ) 1 
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pondant aux noyaux À éloignés, on devrait, en 
principe, observer des raies satellites dues aux fré- 
quences de’ résonances différentes des couches 
d’atomes À voisines de l’impureté. Naturellement, 
même aux concentrations les plus faibles qui don- 
nent un signal observable en pratique, les raies 
sont élargies du fait qu’on a alors un grand nombre 


SK % dt Ms EX % 


Fi. 8. — Premières raies satellites de la résonance de l’argent dans Ag Cd, Ag In, Ag Sn et Ag Sb infiniment dilués. 


d’atomes B dissous répartis de façon aléatoire dans 
l’alliage. Dans ces conditions, chaque atome À est le 
voisin éloigné de plusieurs atomes B. On peut 
cependant espérer distinguer les raies satellites 
des premières couches voisines des atomes B. 
Les observations de Rowland semblent confir- 
mer ce point [20]. 


3. Calcul du déplacement de Knight moyen. — 
Nous considérons maintenant un alliage de substi- 
tution AB, parfaitement désordonné, dans lequel la 
concentration c des atomes B est faible, mais non 
plus infiniment petite. Nous admettons de plus 
que les diffusions multiples des électrons au niveau 
de Fermi de l’alliage par les atomes B sont négli- 


geables. Nous indiquons, comme au paragraphe 


précédent, par À; la position du j% noyau À à 
partir d’un noyau B pris comme origine et nous 
affectons de même un indice & aux atomes B. 

A la fonction de Bloch (32) du métal À pur corres- 
pond maintenant dans la matrice de l’alliage, la 
fonction : 


‘ 


Pr) = ualr) Qk(r) 
Où : 
dlr) = kr) + À eikeRi[qx(r — R;) — Qg(r —R:;)] (36) 
i=1 


est la fonction d’onde d’un électron libre d’impul- 
sion k diffusé par les r puits de potentiel corres- 
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pondant aux n atomes B. Un calcul tout à fait 
analogue à celui du paragraphe précédent montre 
que la variation relative d0/(R;)/0 de la densité 
électronique au jè% noyau À des électrons du 
niveau de Fermi est maintenant : 


Sp’(R;) _ 1 


CRT ET 


ar 
[ [UE CR;) PK(R;) — DE(R;) W(R;)] dQ. (87) 


Ce qui détermine le déplacement relatif SX /Æ du 
centre de la raie de résonance, c’est la valeur 
moyenne Sp /o de la variation relative de densité 
électronique sur les N noyaux À de l’alliage : 


SE _o_1$ do (R;) (38) 
ARR A TE ER 
Posons, pour simplifier l'écriture : 
dk(r) — dr) = Fk(r). (39) 


Dans un alliage parfaitement désordonné, tous 
les atomes B sont statistiquement équivalents. 
Dans ces conditions, l’équation (38) devient, après 
sommation sur tous les atomes B : 

Son. 
or Re) 
e N;=e(R; 

EL e—ik.R; ‘ 
Fins SJ AIR —R)FUR)aQ UWo) 


j=1à 


L’alliage étant constitué de N atomes À et de n 
atomes B répartis de façon tout à fait aléatoire, un 


nœud du réseau a la probabilité d’être oc- 


N 
N +n 


cupé par un noyau À et la probabilité ce — 


N +n 
de l’être par un noyau PB ; cest simplement la con- 
centration de l’alliage en atomes B. En trans- 
formant l’équation (40) de façon à y mettre Ja con- 
centration c en évidence, on obtient : 
So Se(R) 


7 AD R#o P(R) 


Dans cette expression, les sommes sont main- 
tenant étendues à tous les nœuds du réseau autres 


Sp(R) . : 
À (R) se calcule à 
parür de la formule (34) où l’on remplace r par les 
vecteurs R du réseau. D’après la définition (39), la 
fonction Fy(R) représente l'amplitude, au nœud de 
rayon vecteur R, de la partie de l’onde d£ qui a été 
diffusée par l’atome B situé à l’origine. Dans 
l’expression (41), le coefficient de c et celui de «2 
doivent être du même ordre de grandeur. Comme 
les hypothèses de notre calcul supposent en parti- 


que l’origine ; la somme 
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+ [2° Fk(R)|° 40. (41) 
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culier que la concentration c est faible, nous abou- 
tissons au résultat suivant : 


Ke > >) 


= (42) 
K R#0 p(R) 


valable si e € 1. Autrement dit : aux concentra- 
tions faibles du soluté, le déplacement de Knigth 
moyen de la matrice varie proportionnellement à 
cette concentration. , 
Sp(R) 


En remplaçant, dans la relation (42), 50) Pr 
sa valeur déduite de la formule (34), on obtient : 


13K 


cK— > (oy sin? m + Brsin 2%). (43) 
Dans cette expression, les coefficients : 
ai = (21 + 1) 2, DER) — UAR)] . 


ts (2141) 2 KR) nkR) 


ne dépendent pratiquement que de la matrice À... 
En effet, pour des concentrations faibles d’atomes 
dissous, les distorsions du réseau modifient peu les 
distances entre noyaux ; on ne commet pas d’erreur 
sensible en conservant aux noyaux À de l’alliage 
leurs positions dans la matrice non déformée. 
L'essentiel de l’effet de taille est inclus dans les 
déphasages lorsqu'on les calcule en remplaçant la 
charge Z par la charge effective Z’ de Blatt dans la 
relation de Friedel (27). Dans la mesure où les élec- 
trons de conductibilité se comportent comme des 
électrons libres, le produit kr, du rayon de la 
sphère de Fermi par celui de la sphère atomique de 
la matrice est constant : kr, — (7 
métaux monovalents. Par conséquent, les quan- 
tités ÆR ne dépendent que du réseau cristallin du 
métal À. Quant aux déphasages », ils dépendent 
essentiellement du soluté, par sa valence Z +1, 
et un peu du solvant, par l'effet de taille. 


1/3 
| pour les 


4. Comparaison avec les résultats expéri- 
mentaux.* — Les principaux résultats actuels 
portent sur des alliages d’argent [19], [20] et 
quelques alliages de cuivre [26]. La structure de ces 
métaux étant cubique à faces centrées, nous avons 
calculé les coefficients et &; pour un métal mono- 
valent ayant cette structure. 

Les fonctions m(kR) — j(kR) et j{(kR) n{KR) 
prennent des valeurs oscillantes dont l’amplitude 
ne décroît que comme l'inverse du carré de la dis- 
tance À à l’origine. Par conséquent, leurs sommes 
sur les nœuds éloignés dans la matrice ne s’annulent 
qu’en moyenne. Pour évaluer numériquement les 
coefficients «et f;, nous avons effectué les sommes 
sur les N’ premiers noyaux voisins de l’atome B 
pris pour origine, c’est-à-dire, pratiquement, sur 


Le RC, À 
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tous les nœuds du réseau de la matrice situés à l’inté- 
rieur d’une sphère de rayon R, = V{(N' + Tr, 
sauf l’origine. À l’extérieur de cette sphère, nous 
avons remplacé les sommes sur les nœuds du réseau 


. par des mtégrales de volume. Celles-ci ont respec- 
_tivement pour valeurs : 


ñ (KR 
4 1) B 


k3 V4 


» _ BilkR:;) 
k3 y 
où a, et 8, sont les fonctions données par les for- 
Ar 
mules (13) et v, SS 
métal À. 

En prenant N' — 140, ce qui correspond à 
kR, = 10 = 5,2 kry, les &œ et À; ne contribuent 
que pour 10 à 20 % aux « et B, et l'erreur sur 
ceux-c1 ne doit pas dépasser quelques centièmes. 
Nous avons obtenu par ce procédé les valeurs 
numériques du tableau suivant : 


rà est le volume atomique du 


TABLEAU II 


COEFFICIENTS DE LA FORMULE (43) 
POUR UN MÉTAL MONOVALENT À GC. F. C. 


CADRE 652=70,2 
a = — 2,4 Bi — — 1,2 
a = —1,9 = 2, 
ou, = 4,9 Ps = 3,5 


Le calcul de « peut se faire avec précision par 
une autre méthode, qui donne «; = 1,052, mais ne 
s'applique pas au calcul numérique des autres coef- 
ficients [25]. | 

Dans le tableau III, nous comparons aux valeurs 
expérimentales de Drain [19], Rowland [20] et 
Sugawara -[26] les variations relatives de dépla- 
cement de Knight, par unité de concentration, que 
nous avons calculées suivant la formule (43) en 
utilisant les déphasages de Blatt. 

L'accord entre ces résultats est satisfaisant dans 
l’ensemble, surtout lorsque l'effet de valence est 
assez fort. [l faut cependant faire quelques réserves 


* pour Cu Zn et Cu Sn, à cause du moment quadru- 
polaire du cuivre. Notre résultat trop faible pour 


Ag Cd est à rapprocher du fait que Rowland 
observe ainsi une diminution du déplacement de 
Knight de l’argent dans les alliages Ag Au, propor- 
tionnelle à la concentration en or, avec 

1 0K 

Ar 0,25. 
Dans ces alliages, où il n’y a ni effet de valence, ni 
effet de taille important, les approximations que 
nous avons faites ne conduisent à aucune variation 
du déplacement de Knight: En négligeant les 
termes en yu.y% dans les équations (20 et (25) 
nous n’avons pas tenu compte du changement de 
caractère s-p des fonctions d’onde avec le remplis- 


sage de la bande de conductibilité. En fait, les 
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TABLEAU III 


COMPARAISON DES RÉSULTATS THÉORIQUES 
ET EXMÉRIMENTAUX 


ALLIAGE 
Ag Cd 0,3 
Ag Cd 0,4 
Ag In 0,7 
Ag Sn 454 
Ag Sb 452 
Cu Zn 0,1 
Cu Sn (92 


mesures de chaleur spécifique électronique mon- 
trent que les densités d’états au niveau de Fermi 
sont un peu différentes dans l'or et dans 
l’argent [27] ; ceci indique que les modifications de 
caractère s-p des fonctions d’onde quand on passe 
du bas de la bande de conductibilité au niveau de 
Fermi ne sont pas les mêmes dans ces deux métaux. 
La dissolution de l’or dans l’argent doit donc créer 
des perturbations qui se traduisent par des dépha- 
sages petits mais non nuls pour les fonctions d’onde 
dans l’argent [28], [29]. Des corrections analogues 
seraient sans doute à faire dans la série Ag Cd, In, 
Sn, Sb, mais, l’effet de valence l’emportant de 
beaucoup, elles ne sont sensibles que pour Ag Cd. 
La discussion de l’annexe I montre que ces correc- 
tions interviendraient dans le sens favorable au 
moins en ce qui concerne Ag Cd. 


5. Élargissement des raies. — Il est évident sur 
la figure 8, que les fluctuations de densité élec- 
tronique sur les atomes de la matrices doivent 
entraîner un élargissement considérable de la raie 
de résonance nucléaire. Comme le calcul de la 
forme de la raie n’est pas possible pratiquement, 
nous calculerons seulement son second moment, en 
supposant que cette quantité représente assez bien 
la largeur de la raie. 

Soit AH le déplacement de Knight absolu et 2, 
la demi-largeur de raie du métal À pur dans un 
champ Æ. Sa demi-largeur dans l’alliage doit être, 
pour le même champ : 

h=VR +R 
où 
h SK\° (x) - / pes. 
HN Lo ei er 

Avec les hypothèses du désordre parfait et de 
concentration c très faible que nous avons faites 
précédemment, on trouve, par un calcul tout à fait 
analogue à celui qui aboutit à la formule (42): 

So\2 de?) 
(©) Æ Mets p°{R) 
Le nombre D ne dépend que de la nature de 


(45) 


=, DD 
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l’alliage tant que la concentration c demeure faible. 
D'autre part, la formule (41) montre que le 


terme (+) est du second ordre en c. Par consé- 
£ 


quent : 


h, == Vhg + cD'AIAH)? pour c<1. (46) 


Le déplacement de Knight absolu AH étant pro- 
portionnel au champ magnétique A appliqué, on 
voit que si la raie de résonance du métal À pur 
était parfaitement fine, sa largeur dans l’alliage 
serait proportionnelle au champ À appliqué. 

Pour largent, 5 æ 0,1 gauss?; le terme 
cD'(AH)? a sensiblement la même valeur dans le 
cas de l'alliage AgCd pour c — 0,01 et 
H = 10 000 gauss ; il devient beaucoup plus grand 
quand on passe à des atomes dissous de valence 
plus élevée (In, Sn, Sb). Il en résulte que la largeur 
de raie doit croître sensiblement comme c aux con- 
centrations de quelques centièmes pour lesquelles 
la formule (42) est encore valable. Nous avons 


Caloule les valeurs de fe) pour Ag Cd, Ag In, 
Ag Sn et Ag Sb et les largeurs k, qu’on en déduit 


0H gauss 
pour H|- 10000 gauss 


© Ag Cd [19] 
æ Ag Cd 
e Ag lin 
o Ag Sn 
æ Àg Sb 


[20] 


0 0,08 c 


F1G. 9. — Comparaison des largeurs de raie de l’argent 
observées par Drain [19] et Rowland [20] aux seconds 
moments calculés suivant la formule (46). 


d’après la formule (46). La figure 9 permet de com- 
parer les résultats expérimentaux de Drain et de 
Rowland aux courbes déduites de ce calcul. 
L'accord est assez bon dans l’ensemble, surtout si 
l’on tient compte des approximations faites dans les 
caleuls et de l'incertitude de plusieurs des valeurs 
expérimentales. 

Le fait de remplacer certains atomes d’argent par 
des atomes étrangers modifie la partie de la largeur 
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de raie due aux interactions dipolaire et d'échange 
indirect entre noyaux. Les courbes de la figure 9 
indiquent que ces modifications sont petites devant 
l'élargissement dû aux fluctuations de densité élec- 
tronique, tout au moins avec les champs magné- 
tiques de cinq à dix mille gauss utilisés dans les 
expériences. C’est ce que justifie un calcul déve- 
loppé dans l’annexe IT. 


6. Alliages liquides. — Le calcul que nous avons 
fait du déplacementde Knight moyen de la matrice 
peut être modifié de façon qu’il s'applique au cas 
d’alliages liquides près de leur point de fusion. Près 
du point de fusion, en effet, il subsiste un certain 
ordre local autour de chaque atome dans le liquide. 
L’arrangement des atomes voisins correspond géné- 
relement à peu près à la structure cristalline du 
solide juste au-dessous du point de fusion. D'autre 
part, à cause de l’absence d’ordre à grande distance, 
la surface de Fermi du métal liquide doit être plus 
voisine d’une sphère que dans le solide cristallisé : 
c’est ainsi que les semi-métaux, peu conducteurs à 
l’état solide, deviennent bons conducteurs par 
fusion, la propagation des électrons n’étant plus 
gênée par la structure des zones de Brillouin. 

Pour appliquer la formule (43) aux alliages 
liquides, nous calculons les coefficients 6; et « de la : 
façon suivante : nous effectuons les sommes des 
fonctions n£ — 7? et 7, n, sur les positions des MW, 
noyaux premiers voisins d’un atome dissous, en 
admettant que la structure locale est restée celle 
de la matrice à l’état solide ; la répartition des 
atomes plus éloignés étant aléatoire, nous évaluons 
leur contribution aux & et 6; en intégrant les fonc- 
tions nf — 7} et nr j, à l’extérieur d’une sphère de 
volume v3 + NV, v, ; v, est le volume atomique du 
solvant À et v, celui du métal dissous B. Pour déter- 
miner la distance d’un noyau B à ses premiers 
voisins À, en tenant compte de l'effet de taille, nous 
avons considéré un cube de volume v3 + v, dans 
le cas d’un environnement cubique centré autour 
de atome B, et un cube de volume v3 + 3 v, dans 
l'hypothèse d’un arrangement cubique à faces 
centrées. En principe, il faudrait prendre pour v, le 
volume molaire partiel du métal B dans le métal À 
à de très faible concentration. 

Nous avons appliqué ce procédé au calcul du 
déplacement de Knight de la matrice dans les 
alliages Na K, Na Rb, Rb Na, Na Au et Na Hg 
sur lesquels des mesures ont été faites [21], [22]. 


- Ne connaissant pas le volume molaire partiel du 


soluté dans ces alliages, nous admettons que 
l’atome dissous y conserve son volume normal. 
À l’exception de Na Hg, les deux constituants de 
ces alliages sont monovalents, mais les électrons de 
conductibilité y ont des masses effectives un peu 
différentes [4]. D’après la discussion du para- 
graphe 1, nous ne devons pas nous attendre ici à un 
très bon accord avec l'expérience car seul l’effet de 
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taille intervient dans notre calcul des déphasages 
suivant la méthode de Blatt. Le tableau ÎV repro- 
duit les déphasages que nous avons calculés pour 
ces alliages : 

TABLEAU IV 


CHARGES EFFECTIVES Z’ ET DÉPHASAGES 
DANS LES ALLIAGES LIQUIDES 


ALLIAGE HE Mo M M2 
Na K° 0,92 2,4 1,43+ 0,541 — 0,227 — 0,044 
NARD 42 02,57 1.85% 075. 0,344 10.082 
Rb Na 0,58 1,44 1,29 0,59 0,10 

Na Au 0,57 1,44 1,29 0,59 0,10 

Na Hg 1,41 1,61 1,82 1,02 0,36 0,02 


+ Pour Z’ < 0, le potentiel est répulsif : Vp — (ko a)?[2 


Les volumes atomiques du potassium et du rubi- 
dium étant de l’ordre du double de celui du sodium, 
nous envisageons deux possibilité pour l’environ- 
nement d’un de ces atomes dissous dans le sodium : 
soit huits premiers voisins, comme dans la struc- 
ture cubique centrée du sodium, soit en environ- 
nement plus compact de douze premiers voisins, 
comme dans une structure cubique à faces centrées. 
Nous comparons dans le tableau V les résultats de 
notre calcul aux mesure de Bloembergen et 
Rimaï [21] d’une part, d’Oriani et Webb [22] 
d’autre part. 


TABLEAU V 
COMPARAISON 
DES RÉSULTATS THÉORIQUES ET EXPÉRIMENTAUX (1) 
1 K 
CHER 
ALLIAGE MESURÉ CALCULÉ 
cac: CI C: 
Na K 0,35 [21] 0,1 0,25 
Na Rb 055121] 0,15 0,4 
Rb Na — 0,25 [21]  — 0,23 
Na Au Ce Eee 
Na Hg 2 ee OA 


On voit que notre calcul rend compte non seu- 
lement du signe, mais de l’ordre de grandeur des 
effets observés. Dans Na K et Na Rb, où le soluté 
est plus gros que le solvant, nous obtenons le 
meilieur accord avec un environnement compact, 

tandis que l’empilement cubique centré doit 
demeurer pour les autres alliages, dont le soluté est 
plus petit que le solvant. Les mesures faites aux 
rayons X par Gingrich et Henderson [30] sur Na K 
seniblent légèrement favorables à nos conclusions 
pour cet alliage. 

(’) Note ajoutée à la correction. Depuis la rédaction du 

manuscrit, Oriani et Webb nous ont communiqué les ré- 


sultats suivants de leurs mesures récentes sur d’autres 
alliages de métaux alcalins : 


10K 
ce K 
ALLIAGE MESURÉ  CALCULÉ 
Na Li — 0,5 — 0,57 4 
Na Cs 0,7 0,67 (avec 12 premiers voisins) 
Cs Na —0,43 —0,3 
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Afin de permettre une comparaison meilleure des 
résultats théoriques et expérimentaux, compte tenu 
des incertitudes considérables qui entachent les 


ÊK 


F1G. 10. — Variations du déplacement de Knight du 
sodium dans Na Au et Na Hg. Comparaison des points 
expérimentaux [22] aux droites calculées. 


mesures, nous reproduisons sur la figure 10 les 
résultats d’Oriani et Webb pour Na Au et Na Hg 
et les droites correspondantes calculées selon notre 
théorie. [l est très probable que les fluctuations de 
densité électronique telles que nous les avons 
décrites ne soient pas la cause principale de la varia- 
tion du déplacement de Knight observé dans 
Na Hg. En effet, le mercure est assez fortement 
électronégatif vis-à-vis du sodium et la dissolution 
du mercure dans le sodium s'accompagne d’une 
contraction considérable du liquide [31]. Le volume 
molaire partiel du mercure est très petit aux con- 
centrations faibles et, d'autre part, le mercure et le 
sodium forment plusieurs composés définis, ce qui 
indique que l’écran électronique du mercure doit 
affecter fortement les atomes du sodium voisins. 
Notre calcul n’est probablement plus valable dans 
ce Cas. der 
réi 

7. Conclusion. — Dans cette partie, nous avons 
expliqué les caractères essentiels des variations de 
déplacement de Knight observées dans les alliages 
en les reliant aux fluctuations de densité électro- 
nique créées par les atomes dissous. Nous avons 
prouvé ainsi qu’elles ne sont pas dues à des inter- 
actions entre atomes premiers voisins seulement. 

Dans le cas d’alliages solides, nous avons rendu 
compte d’une façon quantitativement satisfaisante 
de la proportionnalité des variations aux concen- 
trations d’atomes dissous. Nous avons montré 
pourquoi cet effet est d’autant plus important que 
la différence de valence entre le soluté et la matrice 
est plus grande. 

Nous avons également pu étendre notre théorie 
aux alliages liquides et trouvé un bon accord pour 


les métaux aicalins. SE 
Manuscrit reçu le 17 avril 1959. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


’ 


TOME 20, OCTOBRE 1959, PAGE 784. 


SUR CERTAINS TYPES D'ÉQUATIONS D’ONDES NON LINÉAIRES 
GÉNÉRALISANT LES ÉQUATIONS DE LA THÉORIE DES PARTICULES A SPIN 


Par GÉrARD PETIAU, 


Institut Henri-Poincaré, Paris. 


Résumé. — Formalisme général d’une théorie des particules à spin dans laquelle celles-ci sont 
décrites par des spineurs déterminés par des équations d’ondes, systèmes d'équations aux dérivées 
partielles du premier ordre non linéaires, dont les solutions ondes planes sont des fonctions ellip- 
tiques ou multiplement périodiques de type donné. Exemples d'équations d'ondes non linéaires 
décrivant des particules de spin 0, %/2 ou #. Exemples et solutions complètes d'équations non 
linéaires décrivant des champs unifiés de spin 0 et % ou de spin 0 et #/2. 


Abstract. — The general formalism of a theory of particles with spin.is given where these are 
specified by spinorial functions which are solutions of systems of non-linear partial differential 


wave-equations of the first order. 
hyperelliptical functions of a given type. 


The plane wave solutions of these equations are elliptical or 
Examples and complete solutions are given of non- 


linear equations for unified fields of spin 0 and % or 0 and %/2. 


1. Introduction. — Caractérisation de fonctions 
d’ondes généralisées. 

Dans plusieurs publications antérieures j’ai déjà 
étudié divers types de généralisations de l’équation 
d’ondes de Klein-Gordon [3], [4]. 

Je me propose ici de généraliser la méthode 
utilisée pour obtenir des équations d’ondes non 
linéaires extensions de l’équation de Dirac ou plus 
généralement des équations de la théorie linéaire 
des particules à spin de Dirac-Fierz [1], [2]. Ces 
équations seront des systèmes d’équations aux 
dérivées partielles du premier ordre, non linéaires 
susceptibles de décrire des particules à spin plus 
complexes que celles décrites par l’équation de 
Dirac ordinaire ou ses extensions. 

Je vais d’abord rappeler brièvement les idées 
que j’ai été conduit à prendre pour base d’une 
théorie généralisée des particules élémentaires. 

Dans la théorie linéaire ordinaire le corpuseule 
n’est caractérisé que par une seule constante intrin- 
sèque : sa masse propre réduite y — Mo C/h. La 
fonction d’ondes est alors, tout au moins dans 
certains systèmes de références privilégiés que lon 
peut considérer comme repères propres du corpus- 
cule, une combinaison d’ondes planes, fonctions 
trigonométriques de la variable 


u= ot Kcot— (Ka) — K,av (1) 


avec 
K?—K° 


= = 


(p = %K, W = xk). 


Suivant l'expression de M. Louis de Broglie dans 
son système propre le corpuscule élémentaire est 
essentiellement une horloge associée au corpuscule 
par l’intermédiaire de la constante w. 

On peut penser que cette limitation à une seule 
constante intrinsèque des éléments caractéristiques 
du corpuscule est peut-être insuffisante. 


En effet la plupart des corpuscules élémentaires 
sont instables et possèdent une « vie moyenne » 
caractéristique. D'autre part il semble qu’une 
théorie satisfaisante des corpuscules élémentaires 
doit pouvoir unifier leur description en présentant 
les différents corpuscules comme états ou dégéné- 
rescence d’une représentation générale qui ne 


semble pas pouvoir être obtenue par les théories. 


linéaires considérées jusqu'ici. 
J’ai alors été amené à formuler l'hypothèse 


qu’une représentation plus adéquate des cor- 


puscules élémentaires serait peut-être obtenue en 
considérant que dans le système propre la fonction 
d’ondes, onde plane élémentaire la plus simple, 
serait non plus une fonction simplement périodique 
mais une fonction multiplement périodique. 

Dans le cas le plus immédiat où la fonction 
d’ondes reste uniforme, celle-ci au lieu d’être une 
fonction trigonométrique dé u serait une fonction 


elliptique et si l’on se restreint ce qui semble physi- 


quement nécessaire aux fonctions régulières d’am- 
plitudes restant 
tique de Jacobi. 

Dans les cas plus généraux, les fonctions à consi- 
dérer, multiplement périodiques cessent obliga- 


partout finies, une fonction ellip- : 


toirement d’être uniformes. La ciasse des fonctions 


les plus simples que nous rencontrons est celle.des 


fonctions hyÿper-elliptiques rentrant dans la caté-. 


gorie des fonctions abeliennes. 
Si l’on considère l’équation de Klein-Gordon 


1,0 d? d? d? 
ES (È dE dE SE =) on 7] Fiæ;t).-= oo) 


les solutions ondes planes les plus simples Ÿ{u) 
sont déterminées par l'équation différentielle 
dy 
du? 


+ Vu) = 0 
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dy ———— ALU 
— = VC — w? Î —— 
0 ) w Ho y, VC2 > 


F = Cosin (u — us + Eo), 


PSE sine ECHO Wa) 

Fest obtenu ici par inversion de l'intégrale 
abelienne la plus simple. 

Dans le cas général où l’équation de Klein- 
Gordon est remplacée par une équation telle que 


D + F(Y) = 0 (3) 


les ondes planes W(u) sont déterminées par une 
équation différentielle telle que 


dY ; 
er AT) = 0 
ou 
dw\2 u dy 
l 2j FF) du — Cu 
ou encore 
“he = VG(Y), 
U 
h 
Re le (a) 
To VG(Y) 


Si la fonction F(Y) est un polynôme de degré 
impair soit 2n — 1 de la fonction Y, la fonc- 
tion Wu) obtenue par inversion de la dernière 
intégrale sera une fonction hyperelliptique possé- 
dant 2n — 2 périodes. 

Réciproquement, on peut considérer F(u) donnée 
a priori comme fonction hyperelliptique d’un type 
déterminé caractérisé par la localisation dans le 
plan complexe de la variable W d’un certain nombre 
de points singuliers. W(u) sera alors complètement 
définie et sera une fonction hyperelliptique dont 
les périodes seront déduites à partir des racines 
de G(T) — 0 par des intégrales prises le long de 
contours convenables. 

Si cette méthode permet de construire des fonc- 
tions d’ondes de types très généraux et d'obtenir 
les équations d’ondes correspondantes, la généra- 
lisation des équations de Dirac est beaucoup moins 
directe. 

On admet généralement que la symétrie relati- 
viste et le rôle attribué dans le schéma général de la 
mécanique ondulatoire à la fonction d’ondes exigent 
que l’équation d’ondes ne contienne que des déri- 
vées premières des fonctions d'ondes celles-ci figu- 
rant linéairement. 

Je conserverai cette condition et les généra- 
lisations des équations d’ondes des particules à spin 
que je vais envisager seront encore linéaires par 
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rapport aux dérivées partielles du premier ordre 
des fonctions d’ondes. 


2. Le formalisme linéaire et l’équation de Dirac. 
— Les équations de Dirac sont généralement for- 
mulées suivant divers algorithmes. 

L’algorithme de Dirac utilise les variables 
d’espace-temps, x, y, z, ctet des matrices @, «, her- 
mitiennes et anticommutantes. L’algorithme de 
von Neuman utilise des variables d’univers, x, y, 2 
etx, — ict avec des matrices y, en introduisant une 
symétrie complète entre espace et temps. Ceci 
entraîne notamment la perte du caractère her- 
mitien des observables et aussi masque certaines 
propriétés associées au caractère topologique de 
l’espace. 

L’algorithme des spineurs de van der Waerden 
utilisé dans la théorie des particules à spin, notam- 
ment dans les travaux fondamentaux de Dirac et 
de M. Fierz nous a semblé plus adapté aux généra- 
lisations que nous avons en vue. 

Le terme spineur ayant été depuis les travaux 
initiaux de van der Waerden dévié de son sens 
original tant par les physiciens anglo-saxons que 
par les mathématiciens je rappellerai brièvement 
comment se présentent les équations de Dirac 
suivant cet algorithme. 

Dans le formalisme de van der Waerden l’élec- 
tron n’est pas décrit par un être mathématique Y, 
à quatre composantes, mais par l’association de 
deux êtres géométriques complexes chacun à deux 
composantes soient met y; (l = 1,2). Leurs com- 
plexes conjugués sont jet y; Les variables x, y, 
z, tsont remplacées par des variables spinorielles x;, 
qui sont des spineurs du second ordre. Les x;, sont 
liés à æ, y, z, t par les relations 


e L3 
Fe ee 22 un — I = 21 
Li =t+z—= x", Tia = TZ + Yy — LE, 


TR LR, CRT Re 


(il y a toujours commutation entre indices pointés 
et non pointés). Les indices spinoriels sont élevés 
et abaissés en utilisant le spineur 


Qi 0 4 Un 
cm = [1 0 L an le, o |? El jp = dy) 
Éa= cap 6, 4 = cPA êr, 

La relation 
LRO) 
entraîne notamment 
AE 


Spineurs du second ordre &, et vecteurs sont 
associés biunivoquement par la correspondance 


: Tr (2 
Eim = AU guim = Au 8 im? 


1 um £ . gl 
= — , ONE RENE 
AU — CE im ; Pim ? 


A 
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avec la relation 
Dre —— Û o = È è gutm Evim — 28 
Sum — j 0 ; 52m Fa 0 , 
nous donne 
1 0 1 0 | 
Sam | 0 art So:im — | 0 1 Due K, D - en se 
im 0 Œs 1 (22 im = 0 1 : 
LM (NN ER PA EE En | Pour des spineurs ondes planes ,(4) nous aurons 
. = A, 0 d 1 0 ; te, ù x d 
SD à PR ETS | 1 | dim — du 8m — Au Bim 3g — Him (7) 
i Si 0 M Te OR 0 | et le système (5) devient 
510 — 0 1 o 529 — 0 i ; : 
î 1 DAREEEST ee Pmlu) = %1 Xi (u), (8) 
Érne | 1 0 L 80m | 1 0 | 


: ù : 
Au vecteur gradient Sn correspond le spineur 


ue 
gradient 
dim —= dy Em im — dy gui, de — dy sr" 
soit 
pre dt + de dit (de 3 1dy) 
Hess Far (dx FR S 1dy) dt Li dx 
rene en 
ds AE dy dt ETS) de 


Les équations de Dirac s'expriment par la corres- 
pondance différentielle linéaire entre , et y : 
D om = 1%; À, Xi = X3 Pme (5) 


x, et x, désignant deux constantes complexes. 
Paritérationona 


sm tn 
dd) Pm = %1 dÙ, Xi — X1 X2 Po 


Or 


d’où 
E1 Pm + X1X2 Pm = 0, C'Xi + %1%2 Xi = 0. 
On pose ordinairement 


= en =) 
L 


avec x réel et positif. 

Les solutions ondes planes du système (5) jouent 
un rôle particulier. On les obtient en cherchant des 
spineurs Yi, Pm SOlutiOns de (5) qui soient fonctions 
de la seule variable 


u — Kct — (Kx) = K, 2. 
Si aux vecteurs X, et 14 on fait correspondre des 
spineurs 


et aim — av gun, 


0 
en) (u) = #2 ontu). 


Pour résoudre ce système nous poserons 
Pm(u) = Am filu), Xi (u) = bi filu), (9) 
Am, bi désignant des spineurs constants, f(x) dési- 
gnant une même fonction pour p, et p», f(u) une 
même fonction pour y. et y. Nous avons alors | 


: d 
KP @m 5 fi = a bi fe, (10) 


: ù 
HieEy Su /2 — X2 Am Ji. 
Les fonctions f,(u), f{u) d’une part, les spineurs 


constants 4m, 0, d'autre part sont déterminés main- 
tenant par 


[im GE nr, (11) 
Œn 
’ 1 x 
fra Et, bi = — am. 
de 


Le premier de ces systèmes donnant 


fi = fe = auf: 
nous montre que les fonctions f, et f, ne seront 
bornées que si nous avons &, > << 0 
Nous pourrons sans perte de généralité poser 


4, — 1,4; = — 1 et nous aurons 
EN tm de (12) 

d’où à 
fh = 2sm(ut+Eé),, f; = xcos (u +5). (13) 


Le second système (11) donne alors 
ki KŸ° am = 9 Ki bi Æ — KX1 Xo Gr 
ou 
KE ST am = %1 Xe QU. 
Nous voyons apparaître une nouvelle condition 


K$ = x, x, réel et positif, 
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Um désignant un spineur arbitraire, pour x, % 
réels et positifs, 


Em(u) = À Viry um Sin (u + Ë) (14) 


DE À Vs K7° Um COS (u + Ë). 


La marche suivie ci-dessus nous montre comment 
nous devons procéder pour introduire des fonctions 
d’ondes doublement ou multiplement périodiques. 

Nous avons introduit les fonctions circulaires 
sin w et cos u par des fonctions f.(u), f.(u) caracté- 
risant chacune un spineur @n ou y;. C'est le système 
du premier ordre 


hi = Gifs Î2 = fi 


avec 
di Gore 0 


qui a déterminé le type fonctionnel des solutions. 


3. Définition des fonctions elliptiques de Jacobi 
par un système différentiel du premier ordre. — La 
recherche des fonctions d’ondes et des équations 
d’ondes généralisées nous conduit donc à préciser 
d’abord la forme possible des systèmes différentiels 
généralisant le système ci-dessus. 

Nous avons déjà rappelé que les fonctions pério- 
diques les plus simples après les fonctions circu- 
laires sont les fonctions elliptiques de Jacobi. Or si 
nous considérons des fonctions réelles, la théorie de 
Jacobi introduit trois types de fonctions 


en (u, k), sn (u, k), dn (u, k) 


avec 


dan w 
Pour k—0, enu—cosu, snu—>sinu, dnu—1. 
R—=1 : t} dn u— L 
PRET SN U—>ENn U Se 
Pour k—1, cn + : He 


Nous pouvons alors considérer deux types de 
généralisation du système (12) qui déterminait les 
fonctions sin w et cos u : 

19 La généralisation de (12) où l’on considère 


5 T ë T 
cos w — sin (u +3), SIN NU — COS u Hg): 


Deux types d'équations sont alors acceptables : 
a) À partir des relations 


sn’ u — dn?u cd u — (1 — k? sn? u) cd u, 


sn w 


ee — — (1 — k? cd? u) sn u, 
n?u 


cd’ u = — k’? 
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On voit que Y, —snu, y: — cd u vérifient le 
système 
EU 2 2 2.2 
Yi Ya yiys Ye = — y + k'y y. 


a2) À partir des relations 


nu = —snudnu = — (k? + k?sn?u) sd u, 
en uw 
sd’u = — Fe (4 + k? sd? u) en u, 


on voit que 


DR = CNE Ye = Sd w 


vérifient le système 


p= —kK'y—kRyiy, ya = y + YU 


2° La généralisation de (12) où l’on considère 
que sn uw, en w généralisent sin w et cos u, dn u 
généralisant 1. 

Les relations différentielles 


SD CN U ANA, en’ u = — sn uw dn w, 


Oneu = trASN UIenNu, 


conduisent à considérer trois fonctions principales 


Yi Sd EX) eys=ren (u EL y = "dns. 
telles que 
V1 = YoYs V2 = — Yiÿs VS = — Kyiye 


Plus généralement nous en concluons que le sys- 
tème ci-dessus a pour solution particulière 


or (u + Éo; k), Yz — CN (uw + Eos k), 
Ya = dn (u + 60, À) 


si les valeurs initiales y%, y8, y$ sont liées par les 
relations 


(Gi) + (= 1, {y} + (yo) = 1, 
k(y9)7 + (y2)° = (y3)° 
avec 
k?=1—R 


Ceci nous conduit à discuter la forme des solu- 
tions du système différentiel généralisant le pré- 
cédent soit 


Het Balise fs = — asfifs (15) 


œ, %», «4 étant réels et positifs. 
On trouve alors les solutions 


02 0\2\ 1/2 
filu) = AN er sn (pu + E, k), 


0\2\ 1/2 
folu) = LE = æ(2) ) cn\(o u + E, k),. (6) 
FORCE EE) 


SI 
2 
(ee) 


p° Œ. (2 3 7 (3)? a) L> 
(12)? CRE (13)? a) L3 


Re 
(1)? 3 + (15)? a) Lo 


Ces expressions se simplifient lorsque nous 
posons 


à = Lo) + PDA, À = CO) ce + (8) 1 ]/2 
puis 
M AVE a p = À Vo, 
NE 


Ji 


À 
— = — Va Xo La, 
ho Lo e 


alors nous écrivons 


— À RERO 
AU = NX V'a TS (e u + £, Fa Ce 2), (18) 

A L 7, 
folu) = AV cn (e CT a V Yi X2 Ka]; 

E UE 
al) = dn (e UE, Virus |: 
V & Ê 
On voit que si 
&3 = , k = (DJ, fs = (OLC = Ê 
Va do 


Le système différentiel devient 


dont la solution est 


h=A1Vasnmout+e), fe = AV cos (pu + Ed). 


Si &1, &, «4 sont différents de zéro nous pouvons 
nous ramener à un système différentiel analogue 
arlt5havecrienmi= a." /Pour-cekeuous 
poserons 

1 1 î 


M 1 fa — P2 Re 
Vo G3 oc : Var Go 


P3- 


Les + sont alors déterminés par le système diffé- 
rentiel 


Q1 = Pa Pa, P2 = — Pr Pas P3 —= — Ps Pre (19) 


Les fonctions @,, +, +3 sont liées par les rela- 
tions 


(p1) — (pi) = — (pa)? + (p2) = — (3)? + (9$)2. 
(20) 
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Les solutions s’écrivent 


À 
1 — À Sn ( u + Ex) (21) 
4 À 
2 — À Cn Le u + Go; à), 
M 
®e — À9 An (uu LS és). 
avec 
x = (99)? + (92), x = (o1)? + (63), 
2 (o9)? QE (3)? Ni CAE tee À Bal (o2)° 
16? —= 7 02 2 …. f 0 2 == 12? in Co = 22° (22) 
(o1)° + (o3) 12 (0) 
Nous avons supposé 11 
(eo)? < (02)? < (08). (23) 


4. Introduction d’un système différentiel du 
premier ordre définissant des fonctions multi- 
plement périodiques. — Nous voyons facilement 
que le procédé suivi ici pour introduire des fonc- 
tions doublement périodiques se généralise pour 
introduire des fonctions multiplement périodiques. 

Nous sommes conduits à considérer le système 
de n équations différentielles du premier ordre non 
linéaires : 


fi = €E1@fo -.. fn 
fear ne ER (24) 
1 — En En J1 - ++ ]n—1. 
1 Lo, -.., On SOnt des constantes > 0, €;, 
EN ere M POSAIE 
Îp = Xp Pr 
avec 
Xp = (op) 0/2 (ati do. am) V4, 


les fonctions @,(u) sont déterminées par le système 


Pp — Ep O1 --- Pp—1 Pp+1 - -- Pm (25) 


ou 


Pp Pr — En P1 Pa... Pp .-. Pre 


L’équivalence entre les systèmes (24) et (25) 
n’existe évidemment que si tous les &, sont difté- 
rents de zéro. 

S1 l’un des & soit « devient nul nous sommes 
ramenés au système différentiel de rang (n — 1), 
les à, (p —1,2, ..., n —1)étant remplacés par 
les valeurs (x, f}) en désignant par f? la valeur 
constante que prend alors la fonction f,. Ce résultat 
montre comment par confluence les «, tendant 
successivement vers 0 notre système différentiel se 
réduirait au système le plus simple f, — «fs, 
le = — fr ; 


484 
1 
0 


MEUTEE € 
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Le système (25) s'intègre immédiatement en 
. donnant 


L 


E1 Pi P1 — E2 Pa P2 = ... — En On Pn 


d’où 
PRES + Dr d Le 2 
1 1 Bi E2 P2 + Bo = . — €p Pr + Pr. (26) 


D'autre part 


2 2 

CPRENO INT pra Do+i at. où 

ou 
2 2 
op = IT egl(ey op + By — Po) (27) 
20 
d’où 
d 
du = —— = (28) 


Fi eal Ep ep a — Bo) |"* 


Pour que les fonctions @,(u) soient acceptables 
du point de vue de la mécanique ondulatoire il nous 
. faut montrer qu’il existe au moins un choix des 
. constantes B, et des signes e, pour lequel toutes les 
fonctions œ,(u) sont bornées. 
Le choix des signes 


Et = > ls £o == 
et la condition 


=) en 1020) 


COR TS, (30) 


permettent de satisfaire à cette condition. On voit 
alors que 


CG = 1 (Be — 81 — 
sera positif sous la seule condition qi < PB; — f.. 
(42)? = (8: Hi on çà) Le (cè UE Poe B) 


sera positif sous la seule condition q5 < B,; — f, ct 
DL 
[IT 


(o)? = (Bo no o%) [ — (Bo ÉCS Fa] 


2<Q1< 


ñ 
IT (oë des: Bu) 
d=?+1 


avec la condition 
Bp — Br < D < es Cr 


(B, — B2est la plus grande des différences PB, — B,) 
Les conditions ci-dessus sont réalisables si l’on 
prend des valeurs initiales convenables. 
On prendra 
Bi = — (9), ..., Bo = + (on) 
avec 


lo) (on) << (ou). 
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Pis Pa; -.., ©, Seront alors astreints aux conditions 
Con) < (on + (9), (en)? < (90)? + (02)2, 
Co — (9) < (on)? < (op)? + (2). 

On pourra prendre sans perte de généralité 
Ba = (1) = 0. 


La fonction @, sera obtenue en inversant l’in- 

tégrale 
"Mes 1e RE LD A A 2 nee 
Jo (ea — 9?) (Ce) — op) ….. (2 — o2)17? 
(31) 


et les autres fonctions à partir de celle-ci par les 
relations 


CRE (02)? ae oÀ. 


9. Formalisme général d’une théorie non linéaire 
des particules à spin. — Les considérations précé- 
dentes me conduisent à admettre que dans la 
théorie non linéaire introduisant les fonctions mul- 
tiplement périodiques les plus simples c’est-à-dire 
les fonctions elliptiques il faudra nécessairement 
introduire trois spineurs [5]. 

Nous écrirons ceux-ci 


ÉAB, ncD, GEF, 


À, BP, C, D, E, F désignant des groupes d’indices 
spinoriels que nous préciserons. 

Lorsque nous considérerons des ondes planes, 
ces spineurs seront des fonctions de 


DR RUE =; Re 
et j’admettrai qu'à chaque type de spineur est 
attaché une fonction de w soient f,(u), fa(u), f3(u), 
indépendante du spin (c’est-à-dire la même quelque 
soient les indices du spineur). Les ondes planes 
associées aux trois spineurs ci-dessus ont donc les 
expressions 


Cer = cprfslu) (32) 


fu fo, fa sont liés par les équations qui caractérisent 
le type fonctionnel, pour des fonctions elliptiques 
nous écrirons le système (15), «;, &, «4 réels et 
positifs et lorsque ces coefficients seront différents 
de zéro, en les incorporant dans les constantes spi- 
norielles a, b, c, 


A di ile 


Pour les ondes planes les équations d’ondes sont 
ramenées à des équations différentielles spino- 


rielles. 
Le premier membre de ces équations de la forme 


Ca — aaBfilu), vob = bob folu), 


(33) 


L, Em 
don C4) 
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((A) désignant un groupe d’indice pointés ou non 
pointés) s’écrira pour l’onde plane 


PS AA . 
4) A 4) f1lu) 


K? ar fi(u) = Ki, am lier (34) 

Le second membre devra donc, afin de ne 
dépendre de w que par f,(u) fA(u) contenir obliga- 
toirement les spineurs n et © chacun une seule fois 
et selon une combinaison présentant la variance 
totale ({(A)). I] s'écrit donc 


: . \l 
%1 F(rob, CEr)(4), 


x, est une constante. F sera un spineur contracté, 
symétrisé ou antisymétrisé des deux spineurs cr, 
Cp qui ici sont de la forme 


cb — bob folu), Car = cr falu). 


L’équation précédente s’écrit alors 


1 Er 
Ke 0n filu) = %a F(bob, ra) fa fe. (35) 


Les spineurs constants a, b, c s’obtiendront à 
partir du système 


IE et 
Ka, = %1 F(bcb, cEr)(a), (36) 


et des équations analogues pour les autres compo- 
santes. 

Les trois spineurs définissant le champ corpus- 
culaire seront donc déterminés par un système de 
trois équations que nous écrirons 


1 ï l 
d AA ru Fin, C)(4); 
. eue & 
din TB) = %2 Pal, É)(B), (37) 
D ne : 
Des (ex 3 PE, n)(C); 


F;, F3, F4 sont des combinaisons de deux spineurs. 
(4), (B), (C) des groupes d’indices, »,, x,, #3 trois 
constantes réelles. 
Pour les ondes planes, nous écrirons 
Lo = Co Alu), 1%, = dE fu), Ce, = cf, falu), 
(38) 


et les fonctions fu /» fa étant déterminées par les 
équations différentielles précédentes les spineurs 


constants 4°, D", "(0 devront être solutions des 
équations 


.1 e 
A! au, — 4 HD, c)ay, 
-l # = 
KE Em = x, F(e, a)!B), (39) 


m  (B) 
L'Z 1 
KG = %3 Fsla, bio. 


Ce n’est que si ce système admet des solutions 
que le système (37) sera acceptable. 
e même nous introduirons des ondes planes 
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multipériodiques eñ considérant 7 spineurs que 
j'écrirai 
a) 
Ë (an 
p —=1,...,n;les (Ap) sont des groupes d’indices 
spinoriels pointés et non pointés. 
Ces n spineurs seront déterminés par nr groupes 
d'équations qui s’écriront 
AU (D (D 


(2+1) Er 
dm © (4?) = MES ARNO 


ÉJ{an. (40) 


Les x, sont » constantes réelles fixées à l’avance. 
Les F, sont des combinaisons par contraction, 
symétrisation ou antisymétrisation contenant 
chacune les n — 1 spineurs autres que & et chacun 
seulement une fois. 

Les ondes planes s’obtiendront en écrivant 


(D) » (D) ,» 
Émp = 4 (4p fou) (41) 
avec 
fo = Ep Up fafo +: fp—afp+a - +: fn. (42) 


Le système sera admissible du point de vue des 


ondes planes si le système déterminant les cons- | 
(D). 

l m c 

tantes 44,» SOIb 


() (P—1) (2+1) (2) 


1 (D) 1 
w a J(am, (43) 


. x 
Km &(An — ep? Fy(a, nc (3 OR EL PEU 
p 


admet un système de solutions autres que O0 pour 
(D) 
tous les a. 


6. Équations d’ondes non linéaires représentant 
des corpuseules de spin 0. — Nous allons préciser 
ces constructions sur quelques exemples qui nous 
serviront de guide dans des constructions ulté- 
rieures. 

Considérons les trois spineurs 


l l L 
ê;, Np» (é 


L’équation la plus simple du type indiqué ci- 

dessus serait de la forme 
1 À 1 ,° 
d END ne Er, (44) 

Mais bien que cette équation possède une 
variance spinorielle fixée, il n’en est pas de même 
de sa variance tensorielle. 

Nous fixerons cette dernière en considérant des 
combinaisons symétriques ou antisymétriques des 
indices let p soient (lp) ou (lp). 

Considérons le cas antisymétrique, nous pren- 
drons alors comme système fondamental les trois 
équations 


Ve Less (DID 
dm EP — = (ré C2 — nt) 
ne 3 NC DE 
dm NP — (££ Eap — ae ca) (45) 


a . #3 (,1 . DES 
Din Co = (5 ri — EG nt). 
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Les solutions ondes planes seront de la forme 


Emp = am j(u), nm? — bn f,(u), Cp — cv f;(u) (46) 
avec 
# 4 | € 
fi = fofs, fe = —fifs Î8 = —fifs 
Les spineurs 
amp, bmp, cmp 
sont déterminés par les équations 
4 . a es DS 
Km am =. — (CAE) 
l . Le ss De 
Km bmP — 2. (cc ag — cy aü) (47) 
l . 4 ne D . 
Ky emD — _—. (a baD — ag ba). 
Pour résoudre ce système nous poserons 
amp — c, KP, br — ©, KP, cp — c, Kw. (48) 


Par identification nous voyons immédiatement 
apparaître les conditions 


He OO) Hands M) No e ARU (49) 
Si nous prenons alors 
a > HE 0) LEO; 
les constantes c,, ©. c,; sont données par 
1 1 1 
CC CS + pol 
a [ke %3l V bei %3l V ka %ol 


Par suite, nous obtenons la solution générale du 
système (45) au moyen de trois constantes arbi- 
traires : 


247 — ni Km» sn Lu UT A à) 
VIe %3l 2 
: ke . À - 
nm? — Km» cn ( u + Es) (51) 
Vlr Xs| d. 
: Fe : RS 
MP — Km» dn (au + ES): 
V'|x1 Xl À 


La valeur de X, est ici encore arbitraire. Mais si 
l’on remarque que 


u — Ky 26 = Koeuam = Ko 


e, désignant un vecteur unité, nous voyons qu’en 
posant 


ue een 
à 
Ê = de NV 
la solution générale ci-dessus s’écrit 

Em — 2 K'm» sn (v + 6, à Vhel) 
Vlrel 

na? — P K'i en (v + €, À Vhel) (53) 
Vlil 

: mp En 
Cmp — dan (v + Éo: à Vel). 


1 Xe 
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Ici 
DRE Re im 
= ms ar —= 2 im (A HAL 
Pour 
Xg = 0, Kmp — Cite — se : 
Ve X»| 


Le système étudié ici est susceptible d’une écri- 
ture tensorielle simple. Il correspond en effet aux 
équations vectorielles 


du AU = x1 Bu CV, du Bu = x: Cu AU, dy CH = x3 Au Bu 


(54) 
liant trois vecteurs AÀ#, Br, C. 
Posant u — K, x" 
Abu) — av fi(u), Bu(u) — bu fu), Cu(u) = cv f;(u) (55) 


les constantes at, b, c! sont déterminées par 

Kua= bu, Hybt=xen at, Kuct = x3 au bte (56) 

sous les conditions (49) 

et avec le choix x, > 0, x < 0, x3 << 0 on trouve 
Kv Kv Kv 


me ——————— , ) re RE = i}) cu FE , 
Vh Xal Vha X3l V le %o| 


ab 


(57) 


et par suite pour À, BY, CY la solution précédente. 

Le système (45) ou (54) décrit un corpuscule de 
spin 0 mais ce système n’est pas la généralisation 
du système usuel qui s’écrit 


| . 
ou en spineurs 
à, 47) = ko en di, ne = Elo. (59) 


Nous généraliserons ce système en considérant 
trois spineurs du second rang 


Én MCP o te 


ou un vecteurs J4 et deux invariants J, et J.. 
L’équation d'ondes du corpuscule considéré « de 
spin Ovest alors 


, U 

, Emr = 1 ne Car 
5: _ m?D ia -? 
D nf = x, Cf ee (60) 

: ia -P 

de ce = 2% n° (4 

ou en notation tensorielle 

du Ju Ii; du Li = #2 Le Ju, du Le = 43 Li Ju. (61) 


Intégrons ce système dans le cas des ondes planes. 
Nous avons encore ici avec u — K, 2" 


13 = aifslu), Ju = ju filu) (62) 


I, = ai fau), 


792 


Les constantes jt, a,, a, sont déterminées par les 
équations 
Ky JU = 1 Qu do, Hu 1 = — %2 Ge lu Eu Ge = —%3 Gi]u 


(63) 


Nous sommes encore conduits à prendre 


K12- 00; Xo = — |x|, Ha = — |X3| 
et nous obtenons 
K; KG Ky 
1 nn or  — = ———. (64) 
Vi %3l V1 %| V'hte %3l 


La solution générale « onde plane » du système 
est alors 


À À 

Ju(u) = — Ky sn Le un Es) 
Vlk Xl \2 
1 € À 

FAO RO 6. NES =) (65) 
V ki %3l Xe 
À À 

Tu Due an (. u + Ëo, ©): 
Vhu %| ko 

Si nous posons alors 
2 Ko D He ho Ky = K; 
DU RC K an (66) 
x Re” 
k Nue À Vhl, 


cette onde plane s’écrit encore 


À ur 
Jatv)e= ae SOUDE A Ve 
2 
1, (v) = 7 He en (DES 1 Vhel), (67) 
1 
Lo) = 0 dn (o + Eo, à Vpal). 
V [1 %ol 


Sur cette dernière forme on voit comment s’effec- 
tue la réduction au cas linéaire. En effet si x3 — 0, 


EN / bel 0. 
I, se réduit à la constante 


K 


15 = 68 
FAT [1 Xe) fo 
et 
Cr | 
Ju(v) = -—=sin (v + 6) (69) 
Vlr 
À K° 
Tv) = — cos (v + 65) 
Wal 
avec 


DA RACE 


Le système étudié ci-dessus et ses solutions se 
généralisent facilement au cas multipériodique. 
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Nous considérons en adoptant la représentation 
vectorielle un vecteur J# et n invariants, 1;, 
JR 1, nes parles net 1 équations aux dérivées 
partielles du premier ordre (p = 1, ... n) 


du Ju(at) = %o Li Lo .….. 17 
de le = lie US 
du Ip — Xp Li 


(70) 
. Te To+i muette 1F Ju. 


Les solutions ondes planes s’obtiennent en écri- 
vant que J4, 1,, ... 1, sont des fonctions de la 


variable 
NT En 
alors 
en 
Ceci nous conduit à poser 
Julu) = jupolu), Zo(u) = ap Pplu). (71) 
Jus Ets ce. Up, +. An SON des constantes. 


La méthode générale précédente nous conduit à 


déterminer les fonctions p(u), ..., @p(u), ... par 
le système 
Qolu) = Pi Pr +. On 
Eylu) = —— Po Ps + Pp—1 Po+i - ++ Pr (72) 
VD 0 500 6) 


Les constantes j,, 41, .. 


., 4 S’obtiennent alors à 
partir des relations : 


TE Can Gt Se 5 Ce 


Ky Ap = — Xp 3 . .: An Ju (73) 


La recherche des solutions de ce système conduit 
aux conditions 


Gp—1 Ap+1 


0 7 0, Se Pépl, (p = 1, SES n). (74) 
On obtient alors 
: 13 K$ 1 Fr 
E Ne . . Hu, Ge 


= : 1 
Xl ( Ho ne 
Xo \fk1%9 ... Xl ; 


7. Equations d’ondes non linéaires représentant 
des corpuseules de spin %. — Nous pouvons égale- 
ment construire dans notre théorie des équations 
généralisant le système ordinaire décrivant les par- 
ticules de spin À. 

Ces équations s’écrivent dans le cas linéaire 

a) En notation vectorielle 


du Tv — dv Ju = %1 Fos, 


Gp —= 


dy FU Lo Ju (76) 


en introduisant un vecteur J4 et un tenseur anti- 
symétrique du second ordre F4”, 


| 


. 
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b) En notation spinorielle 
da 14, + De J = 2x1 Fim, 
di Fam = x, Ji, 


en introduisant un spineur du second rang symé- 
tique Fm et un spineur du second rang J". 
Dans notre schéma ces systèmes sont susceptibles 
de deux généralisations : 
19 Une généralisation conservant le type spino- 
riel ou tensoriel. 
Celle-ci s’écrit en notation vectorielle en intro- 
duisant à côté du vecteur J# deux tenseurs anti- 


symétriques du second ordre F,, et G,, liés par les 
équations 


du US Fa dy Ju —= x1( Fu Gov = Fe Gou); 


dy Fuv = x, GX, Je, (78) 


dy Guy — X3 EE, J®, 


Sous forme spinorielle le spineur J® sera lié à 
deux spineurs symétriques du second rang F,, et 
Gym par les équations 


Dia 19, + Da J8 — 21 (Fig im + Fma 4) 


Di FE = %3 Grp JD (79) 


OCR Fiv JD, 


La recherche des solutions ondes planes de ces 
systèmes s'effectue suivant notre méthode générale. 
Considérons la représentation vectorielle, nous 
posons 


Ju = Ju. fi(u), Fyv = fuv fou), (80) 


Cuv — £uv fa(u), 


les fonctions f,, f2, f, sont solutions de notre sys- 
tème général. 

Les vecteurs et tenseurs constants 4. CU RS 
sont déterminés par le système 


Eu jv — Ky ju = (fur Lov — fw° ul, 


KE JE = — x 20e Je, (81) 


BL : 
Ky ge = — x, f eo JP. 


La condition de Lorentz à, J" — 0 du système 
classique (76) est ici remplacée par les conditions 


Fu J,) = 0 et Guy ds) = 0 


qui sont identiquement satisfaites. 

20 Une seconde généralisation semble beaucoup 
plus intéressante. Elie permet de considérer les 
équations linéaires des corpuscules de spin 0 et de 
spin % comme deux dégénérescences d’un même 
système non linéaire. 

Nous introduisons un champ général caractérisé 


TYPES D'ÉQUATIONS D’ONDES NON LINÉAIRES 


793 


du point de vue tensoriel par un invariant 1, 
un vecteur JY, un tenseur antisymétrique du 
second ordre F4, Ces trois grandeurs sont liées par 
les équations 


du JS = dy Jy —= X1 TFyv, 


2, Fuv = x, IJu, (82) 
dl IN 
On a encore ici 
du JU = 0 (83) 


car en multipliant la 2€ équation par à, on est 
conduit à 


H9 I dy JU — 0. 


Pour x3 —= 0 on a 1 —Cte — J,le système (82) 
dégénère et se ramène au système des équations du 
corpuscule de spin * (76). > 

Pour x: —=0et F,, = Cie —},,, posant 


Sun 
le système se réduit à 
du 1 = %3 Sy 
du SM — %3(fUUNT juv) Z (84) 
ou encore 
CG 1 = x 43 JEUN] ftuv JL. (85) 


On retrouve les équations du corpuscule de spin 0 
pour %1 X3 J fran < 0. 

Nous obtenons facilement les solutions ondes 
planes de (82) en écrivant 


Ju = Julu) = ju folu) (86) 
Fuy — Fuvlu) = fus fafu), 
fa fo, fa étant liés par les relations 
= Erfofs, fe = tofifs fs = esfifr, (82) 
avec 
nn NC Ca Nu 3) 


Les constantes L,, j,, f, sont alors déterminées 
par les équations 


Ky jy — K, lu = a X1lo fav 
K, JU = es x2 lo jb, (88) 
Ky Lo = €3%3 uv J" 
et 
Ky ju = 0. 
On trouve immédiatement la condition 
Enes Hyae Ù 
et 


KE — [1 %ol 18. 
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Lo dernière équation donne immédiatement 
€1 X1 15 — Es X3 Ju JM. 


De même les deux dernières équations donnent 
la relation 


Pr tag 
KŸ = — € €3 40 H3 Ju JV. 


Si le vecteur 7, est du genre temps soit 7, j > 0, 
nous obtenons le système de conditions 


CE Hs Has Ùn  Extes ka Rat Des es rie US 00) 
Si nous considérons par exemple le choix 
€, —= +1,e, ——1,e, — — 1, les constantes x 


seront assujetties aux conditions 


PCT NS 20 Fe Ier LA (90) 


qui pourront être satisfaites soit par x, > 0, 
D = Ù, 22 UNOU par QU, 2 40%, 70: 
On a alors 


KG = lxs noel 16, 


fes 1 2 Pr 
LE ne Eee, 
LE re TE 


(91) 
Supposant les conditions précédentes réalisées 
nous nous donnerons un vecteur du genre temps 
constant 7, tel que 
72 


EAN) TU — 0ù, 
u / ? Ju | ko Xl 


Nous en déduisons 


To = KolV Pr %ol, (92) 
Paca Edit UV 
LV %1 Éo (Au IN = K, ju). 
Avec le choix er => 0,6 = 1,e, = 1, 
%1 > 0, x2 > 0, x3 — — |x;] nous obtenons les 


fonctions ondes planes 


; À 
Ju = À Ju Sn Le ue 


ho 
UK À 
ee ON Ge a), 93 
se x, (93) 
A1 VA Lo : ; À 
Fu — LEE Fes (Æy Jv LER AS ju) cn Le U, + 


Mais X, est jusqu'ici arbitraire. Nous pouvons 
alors redéfinir w en écrivant 


Aou = À, Kyat = Kjat = v (94) 
en posant 
PATES 
De même nous redéfinirons un vecteur constant 7, 
par js = V42 li j,, arbitraire sous les conditions 


KL IU—= 0, pu ARE, (95) 
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avec 


Alors la solution onde plane complète considérée 
s'écrit 


ie. an (v, à Vlx3l) 
Vr Lo 
ee Je sn (», AV) (96) 
ve 
1 et 7 Es 
Fus = KA, Ve (AG 5 — 5 jé) en (v, à VPesl). 


À et K, sont ici arbitraires ainsi que 7, sous les 
conditions (95). 

Le système non linéaire considéré ici peut être 
interprété comme réunissant en un champ unique 
le champ de spin % électromagnétique ou mésique 
vectoriel et le champ mésique de spin 0. Le choix 
particulier des signes ci-dessus définit une solution 
qui peut être considérée comme voisine de celle du 
champ de spin #, celle-ci subsistant seule dans la 
dégénérescence À V/}x,| —0 quiréduit le système au 
cas linéaire. 

Le choix des signes €, —e; = —1,e, — +1, 
conduirait à une solution voisine du spin 0 et se 
réduisant à celui-ci pour x, = 0. 


8. Équations d’ondes non linéaires représentant 
les corpuscules de spin #/2. — Nous allons main- 
tenant examiner de plus près la possibilité de géné- 
ralisation des particules de spin #%/2. Nous allons 
donc chercher une généralisation de la théorie de 
Dirac, c’est-à-dire chercher à construire une théorie 
faisant intervenir au moins deux spineurs de rang 1 
de la forme Et. 

Nous voyons immédiatement qu’il ne peut y 
avoir de systèmes de trois spineurs Et, #}, €! liés 
par des équations du type précédent. 

On ne peut en effet construire de système tels 
que 


dé, Em — (ne. 


On voit alors immédiatement que, si l’on ne veut 
utiliser que des spineurs de van der Waerden du 
type &' on ne peut avoir de solutions ondes planes 
représentées par des fonctions multipériodiques que 
si l’on considère au moins quatre spineurs. 

Soient 


D @) E) «) 
61, se a A 


2? , 


ces spineurs. Nous désignerons par (a), (b), (c), (d) 
une permutation de 1, 2, 3, 4 


? ? ie 
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Le système le plus simple que nous pouvons 
écrire est alors 


L (a) . (b) (c) (D)EA(C) (Ad) 
Du les, Leo eu 
(€) (d) (c). (d) () (4) (b) (d) (b) (c) 
( EH UEr EUL (re Eo a 64 4) ad) 


(97) 


Ce système possédera des ondes planes multi- 
plement périodiques (ici triplement périodiques) 
(a) 


si nous considérons les spineurs Ë! fonctions de uw 
de la forme : 


(a) (a) 
Eu) = Clfalu) 


la — Ea fo fe Ja. 


; (98) 
avec 
(99) 


(a) 
Les C! seront alors déterminés par l’équation 
entre constantes 


l (a). @)  (c) (D). (c) (à) 
K; Cm— Xq [CG a Ca 009 CG ) C1 


(b) (c) 
+[ JOt+[.]Ct] (oo 


9. Équations d’ondes non linéaires associant les 
champs de spin %/2, 0 et %. — Nous avons fait 
l'hypothèse dans l’exemple précédent d’un seul 
type spinoriel £. Mais nous pouvons également 
considérer le cas où deux spineurs El, 1” seraient 
associés à un ou plusieurs spineurs de rang supé- 
rieur par exemple de l’un des types 

ARR 


FA 


Le cas le plus simple est celui d’un champ du 
type C7 correspondant tensoriellement à un inva- 
riant ou à un pseudo-invariant. Notre formalisme 
général nous conduit aux équations 


1 & u 
de Em Ai de fa: 


de, nn = 3 00 En, (101) 


dE CU = net no + EP ni). 


Les solutions périodiques ondes planes si elles 
existent sont de la forme 


Em(u) = am filu), nu) — bm f(u), 
avec 


fi = E1 fo fs; fe — Eof1 fs fs — €3 f1 12 


les constantes a”, b”, c"? doivent satisfaire aux 
relations 


Ki. am — EX de ba 
Kb, — c3%5 € ds (103) 
Ki CP = eyvglal b? + a? bi). 


0 
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Par combinaison nous obtenons 
KG ai = €, e3 x al +) (as bgbt. 
Multipliant par b{ on en déduit 


al by = 0 
d’où 
at D; = 0 


et par suite a! = 0 à moins que l’on ait x, x3 = 0. 
De même la seconde équation donnerait b! — 0 
à moins que X1 X3 = 0. 

Ce n’est donc que si x; — 0 c’est-à-dire 7 — Cie 
que le système considéré admet des solutions ondes 
planes et dans ce cas il se ramène à l’équation de 
Dirac linéaire usuelle. 

L'association de deux spineurs du premier 
rang El et d’un spineur antisymétrique du 
second rang Cl? ne conduit donc pas à un système 
possédant des solutions doublement périodiques. 

Dans le cas de deux spineurs &!, nlet d’un spi- 
neur symétrique €?’ nous écrivons les équations 

d, Em — x, Ca no ;! d, nm — x; CE (10%) 


d, Cm —= x3(E1 nn? = ED n). 


On montre par un calcul semblable à celui du cas 
précédent qu’il n’y a pas de solution du type onde 
plane pour ce système à moins que l’on ait x3 = 0. 

Il nous reste à considérer les systèmes associant 


aux spineurs &, nlun vecteur ou pseudo-vecteur €. 
Deux types d'équations sont possibles 


a) d,Em= mt ne, 


D CT — Don CP = x 3(E7 No — En M). (105) 


On peut encore montrer par un calcul moins 
immédiat que ce système n’admet pas de solution 
ondes planes sauf dans le cas où x3 = 0. 


b) D, Em = v1 CE mp, 
de n% = Ko œ Er, 
Din CP + Don CP = Ma(E No + ÉD M). (106) 


Nous allons obtenir des expressions complètes 
pour les solutions ondes planes de ce système. Nous 
écrirons 
Cu) = cd falu), 


(107) 


Emlu) = amfilu), mu) — bm fo(u), 


les fonctions f.(u), f{u), fa(u) seront déterminées 
par le système 


A — E1f2fs 


Ei —= Ej— — 1, 


f = Co f1 f3 fs — Eg fi fe 
ER = ol (2, k = 1,2, 3). 
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Les constantes a”, b", c? sont alors déterminées 


par le système 


Et, am = €; x d, br, . 
KT, dm — €, x, ch av, (108) 
Em y + Kom ue — €3 X3(t bp + ay bi). 


Ce système se résoud sans difficulté. Supposant 
X1, X2 réels, on est conduit à la condition 


E1 €o X1 Ko > 0. 


Ceci nous conduit à deux types d’équations et 
par suite à deux types de solutions. 


10 

Ey = + 1, EEE | 

#1 = fx, 2 = — |X»| (ou #1 %2 < 0) 
90 

Ey = — 1, €Eo = — 1, Es — 1 

X1 = Pril, DEN 


Dans ces deux cas a" et Ds exprmeront au 
moyen d’un spineur arbitraire soit uw" et l’on aura 


am — um pbm — 


VE Val 


um. (109) 


Pour c on obtiendra avec Es = — À pour le 
premier cas, €; = + 1 dans le second cas 


cm = 5 tx à 
Te #2] KG 5 


Les ondes planes complètes s’exprimeront au 
moyen des solutions (21) pour €, = + 1, 
EE Es — 

Nous serons encore conduits à redéfinir la masse 
propre X, et la vecteur X, (ou le spineur X£) en 
posant 


(410) 


Ko — X° K5; Ur Fe = 


FANS 
= m 


0 


Ne 10 


Nous poserons À,/h, — À V’hel et le spineur un 
étant.arbitraire nous écrirons 


La 
À Um —= Um: 


Nous obtenons alors les solutions générales ondes 
planes du système (106) sous les formes suivantes : 


Dans le premier cas €, = +1,e, —=e; =—1: 
pm Un (n + Es, A Vhal), 
ve : Pesl 
v en (0 + Eo, À Viral), (111) 
: 1 K'MQ uQ ur h 
om. — == CET RPANES 
Rorieders 
X dn (v ai Eos à Vhril), 
et dans le second cas €, —=e, —=—1,e,) — +1" 
en a” (o + 6, À Vhal), 
(112) 


#1 


PA = —ûn (o + &o À VPesi), 


= L (ex K'ma UQ AE 
< Vel LE : 


ren (DEA VEN 


Ce système peut être interprété comme rassem- 
blant en une entité physique unique les ondes spi- 
norielles &!, n! décrivant ordinairement l’électron et 


lPonde C7 (vecteur ou pseudovecteur) décrivant 
ordinairement le champ mésique. 

La mécanique ondulatoire non linéaire permet 
donc de construire des modèles unifiant les concepts 
habituellement séparés de particules et de champs 
d'interaction. 


Manuscrit reçu le 24 février 1959. 
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APPLICATION DES ÉMULSIONS DE FAIBLE SENSIBILITÉ 
À LA DISCRIMINATION DES IONS LOURDS 


Par CH. GEGAUFF et J. P. LONCHAMP (1), 


Laboratoire de Physique Corpusculaire, Strasbourg. 


Résumé. — On expose l’ensemble des faits expérimentaux démontrant la dualité de structure 


des traces d’ions lourds. 


Le calcul de la perte d’énergie spécifique « primaire » a suggéré une méthode de discrimination 
entre ions lourds (Z = 6, 7, 8, 10 et 18), enregistrés dans des émulsions à faible sensibilité (IlfordÆ,, 
K — 1, À — 2). Différents critères sont examinés : longueur ! de la partie continue, comptage 
de lacunes en fin de parcours. L'identification est possible à une unité de charge près. 


Abstract. — This paper presents a review of all the experimental facts showing that heavy ion 
tracks in nuclear emulsions are formed of two distinct parts which we call ‘‘ primary ”” and ‘‘secon- 


dary ”. 


The determination of the specific energy loss of the ‘‘ primary part ” of the track suggests a 
method of identification of heavy ions (Z = 6, 7, 8, 10 and 18) recorded in low sensitivity emul- 


sions (Ilford Æ,, À — 1, K — 2). 


Several means of identification are discussed : length without gaps at the end of the track, 
number of gaps counted from the end of the track. IL is possible to identify the ions with an accuracy 


of one unit of charge. 


Structure des traces d’ions lourds dans les émul- 
sions nucléaires sensibles. — L’explication du phé- 
nomène d’amincissement des traces des ions 
de Z > 2 vers la fin de leur parcours nous avait 
amené à rejeter l'hypothèse de Freier [1] attri- 
buant ce phénomène à la diminution du pouvoir 
ionisant de la particule à la suite de captures élec- 
troniques lorsque l’ion était suffisamment ralenti. 
Nous proposions alors une explication basée sur la 
densité spatiale des rayons Ô de faible énergie, 
éjectés par le passage de l’ion. Ce point de vue 
devait logiquement nous amener à postuler une 
dualité de structure des traces [2]. Celles-ci com- 
portent deux parties que nous appelons «primaire » 
et «secondaire ». La partie primaire est constituée 
par les grains effectivement traversés par l’ion inci- 
dent. La partie secondaire est constituée par les 
grains situés au voisinage immédiat de la trace pri- 
maire et rendus développables par des électrons 
d'énergie supérieure à un certain seuil, dont nous 
discuterons ultérieurement la valeur. La partie 
secondaire peut devenir de plus en plus impor- 
tante au fur et à mesure que l’on s'éloigne de 
l'extrémité de la trace, parce que l’énergie maxi- 
mum des électrons projetés croît assez rapidement 
avec la vitesse de la particule. Cette dualité de 
structure peut être confirmée par l’examen des 
traces enregistrées dans des émulsions nucléaires de 
sensibilités différentes. 

La caractéristique essentielle des émulsions 
nucléaires est leur sensibilité aux particules et spé- 
cialement aux électrons. L’émulsion du type 


(1) Actuellement Faculté des Sciences de Lyon. 


Ilford G; enregistre les électrons de toute énergie, 
sous forme de traces visibles. Le diamètre moyen 
de ses grains vierges est de 0,27 1. Il en est de même 
pour l’émulsion Ilford Z,. Le diamètre moyen de 
ses grains vierges est de 0.13 p. 

L’émulsion du type Ilford C, est moins sensible, 
elle n’enregistre sous forme de traces visibles que 
des électrons d’énergie inférieure à 30 keV ; le dia- 
mètre moyen de ses grains vierges est de 0.16 p. 
La comparaison des figures À montre très clai- 
rement la différence d’aspect d’un même ion lourd, 
ion 59Ne de 200 MeV d'énergie et de 130 u de 
parcours, enregistré respectivement dans une émul- 
sion G;, Li et C2 On remarquera d’abord que 
l'épaisseur moyenne est plus faible dans la €, et 
dans la L, que dans la G;, ce qui s’explique en par- 
tie par la différence de la taille des grains de ces 
émulsions. 

Dans le tableau I, nous avons représenté l’aug- 


lonser sG 10Ne 18À 

NOR 0,48 ue 0,84 be 1,60 be 
PROS NT PE 
PV Ro O18u 033 072 


2 


mentation d'épaisseur moyenne Ae:59:10, C’est-à- 
dire la différence entre l’épaisseur moyenne à 100 
de parcours résiduel et l'épaisseur moyenne à 10 4, 
pour différents types d’émulsions et différents ions. 
Les épaisseurs moyennes sont évaluées à partir de 
20 traces. 
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Rappelons que Ae est influencé par le degré de 
développement et parfois par des fluctuations de 
sensibilité d'un même type d’émulsion (spécia- 
lement pour la G;). 

Néanmoins, les variations qui en résultent sont 
relativement faibles, ainsi Ae passe de 0.80 & 


Fic. 1. — Ÿ0Ne E — 200 MeV, R = 130 y. 


pour le Ne dans la G; développés 20 minutes à 
Pamidol, à 1.11 u pour le Ne développé 120 mi- 
nutes. 

Dans tous les cas, Ae est significativement infé- 
rieur quand on passe d’une émulsion G; à une 
émulsion ZL,, ou à une émulsion C,, pour un même 
type d’ion. Nous attribuons cet effet moins à la 
différence de dimensions des grains qu’à une diffé- 
rence de portée pour des électrons de même énergie 
dans les émulsions L, et G; pour lesquelles le dia- 
mètre des grains passe du simple au double. En 
effet, un modèle de répartition des grains dans 
l’émulsion, l’évaluation statistique de la portée et 
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le parcours des électrons chiffré an nombre de 
grains rendus développables [3], [4], nous ont 
permis de montrer qu’il fallait un nombre plus 
grand d’électrons d’une même énergie, dans le cas 
de la L, que dans le cas de la G; pour atteindre leur 
portée maximum. ae 
Les courbes de la figure 2 représentent l’épais- 
seur moyenne en fonction de 8 = V/c pour eC 
enregistré dans les émulsions G;, Laet C2. À vitesse 
égale, la différence d’épaisseur entre les traces 
enregistrées dans la G; et la Z, croît lorsque la 
vitesse augmente, ce qui confirme nos prévisions 
quant à la portée. Le calcul du spectre à pour 4C 
montre que pour 8 > 0.09, les électrons de 9 keV 
et plus apparaissent, mais que leur nombre reste 
petit de l’ordre de 3 à 4 par p de parcours. Leur 


épaisseur en Jl 


d 42 
—C$— 


Fr 


efficacité géométrique n’est plus la même pour les 
deux émulsions. Ce phénomène est plus marqué 
pour des ions de Z = 6 pour lesquels le nombre 
de à est plus important. 

Dans le cas de l’émulsion C,, il se surajoute à ce 
facteur géométrique une moindre sensibilité des 
grains. La perte d’énergie spécifique initiale des 
électrons de 9 keV et plus, est faible, de l’ordre 
de 5 keV/u et de ce fait, leur efficacité photo- 
graphique dans l’émulsion C, est faible elle aussi 
comme il résulte des données de Brown [5]. Leur 
contribution à l’épaisseur de la trace est quasi 
nulle, contrairement à ce qui se passe dans des 
émulsions à plus grande sensibilité comme la G; 
et la L, dans laquelle leur efficacité photographique 
est importante. Par contre, le nombre d'électrons 
d'énergie inférieure à 9 keV décroît lentement, 
lorsque la vitesse croît à partir de B — 0,09. Ces 


remarques expliquent qualitativement l’allure des 
courbes de la figure 2. 
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Figure 3, nous avons reproduit côte à côte, les 
microphotographies de deux traces d’ions 0 dans 


Fre 8: 


une émulsion G; et D,. Pour la G; Ae —0.79 pu 
et pour la D, Ae — 0.06 pu. Cette dernière valeur 
de Ae, s'explique par la très faible sensibilité intrin- 
sèque aux électrons des grains de l’émulsion D. 
La différence entre partie primaire et secondaire 
peut encore être démontrée par le sous-dévelop- 
pement des traces ou par le développement à l’aide 
d’un révélateur discriminant. Il paraît raisonnable 
d'admettre que la période d’induction du dévelop- 
pement de l’image latente des grains traversés par 
la particule incidente, donc des grains de la partie 


APPLICATION DES ÉMULSIONS DE FAIBLE SENSIBITÉ 


799 


primaire, soit plus courte que celle des grains tra- 
versés par les seuls à. Un sous-développement par 
dilution et par abaissement de température devrait 
permettre de saisir ce phénomène, en arrêtant le 
développement avant l’apparition de la partie 
secondaire. 

Le sous-développement à l’aide de l’amidol ne 
permet pas d’attemdre ce stade dans des condi- 


RS ESER 


Fi1G. 4. — ?0Ne, E — 200 MeV, R = 130 y. 


tions reproductibles. Nous avons choisi un révé- 
lateur au glycin pour ces expériences. La figure 4(a) 
montre une trace de 2%Ne dans une G; sous-déve- 
loppée au glycin par dilution et abaissement de 
température. La variation de la dimension trans- 
versale est très faible, la partie lacunaire provient 
essentiellement de la partie primaire de la trace. 
Par contre, un révélateur à l’oxalate ferreux à 
potentiel redox élevé (— 225 mV) qui développe 
sélectivement les germes d’image latente en fonc- 
tion de leur taille, nous a donné un résultat encore 
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plus frappant (voir fig. 4b) ; la partie primaire de la 
trace est nettement prépondérante. 


Structure des traces d’ions lourds dans les émul- 
sions nueléaires de faible sensibilité. — Les résul- 
tats précédents nous ont incités à examiner l’aspect 
des traces d’ions lourds dans des émulsions prati- 
quement insensibles aux électrons, du type X5, 
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Fic. 5. — Ÿ$Ne, E — 200 MeV, R — 130 v. 


K—1, K— 72. Le diamètre moyen des grains 
vierges de ces émulsions est de 0.2 p. Leur sensi- 
bilité, indiquée par la perte d’énergie nécessaire 
pour rendre un grain développable, va en diminuant 
de 5 à K — 2. Si nos hypothèses sont exactes, 
on devrait s'attendre à la disparition totale de la 
partie secondaire de la trace, due aux électrons à. 
La figure 5 montre un exemple des résultats 
obtenus. Il s’agit d’ions Ne enregistrés respec- 
tivement dans de l’émulsion K,, K — 1 et K — 2, 
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et développées à l’amidol. On constate que la trace 
est lacunaire à son entrée dans l’émulsion, que la 
densité des grains croît régulièrement et que les 
grains finissent par devenir jointifs, la partie termi- 
nale de la trace étant continue. Cet aspect rappelle 
celui des traces dues aux particules légères. L’allure 
de la variation de la granularité est certainement 
liée à la variation de la perte d’énergie spéci- 
fique dE dx. 


a) ÉVALUATION DES PERTES D'ÉNERGIE {SPÉCI- 
FIQUE. — Le dÆ/dx total peut être calculé, d’une 
part à partir des relations parcours-énergie expéri- 
mentales par dérivation graphique des courbes, 
opération assez délicate, et d’autre part à partir 
du d£Æ/dx d’un proton de même vitesse en multi- 
pliant cette valeur par le carré de la charge efficace 
de l’ion à la vitesse considérée. Nous utilisons les 
valeurs de Z efficace proposées par Papineau [6] 
pour Z < 10. 

Cette perte d'énergie ne contribue pas intégra- 
lement à l’ionisation dans la partie primaire de la 
trace. En effet, une partie de l’énergie dissipée est 
communiquée à des électrons qui acquièrent une 
énergie suffisante pour provoquer de l’ionisation en 
dehors du domaine de la trace primaire. Phénomé- 
nologiquement, on peut dire que les électrons rece- 
vant une faible énergie, c’est-à-dire la grande majo- 
rité, ne provoquent de l’ionisation que dans la 
trace primaire, alors que les électrons recevant une 
grande énergie (E > 5 keV) ne cèdent qu’une partie 
très faible de leur énergie dans la partie primaire. 
L’ionisation due aux électrons d’énergie inter- 
médiaire se partage entre la partie primaire et la 
partie secondaire de la trace. Ce phénomène de 
partage étant continu, il est arbitraire de fixer des 
limites à ces énergies électroniques. Il semble assez 
raisonnable d’admettre au vu des relations par- 
cours-énergie et énergie-grains développés [5] que 
des électrons d'énergie Æ < 2 keV ne contribuent 
qu’à former la trace primaire. En adoptant cette 
limite, on sous-estime légèrement la perte d'énergie 
spécifique que nous appelons primaire. Si nous 
adoptons la limite Æ < 5 keV, on risque de la 
surestimer. Désignons par n, la limite choisie pour 
le transfert d'énergie. La probabilité d’avoir un 
transfert d'énergie compris entre Q et Q + dQ a 
une section efficace différentielle de 

de — 2% #40 

mr O 

le calcul de cette section est légitime dans la 
mesure Où n est grand par rapport aux énergies de 
liaison des électrons atomiques ; ceux-ci peuvent 
alors être considérés comme libres. Par ailleurs, 
n ne doit pas être trop grand par rapport aux 
dimensions atomiques, de sorte que la particule 


chargée incidente puisse être considérée comme 
ponctuelle. 
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À No 10 


Par intégration de do sur tous les transferts 
possibles, on obtient la perte d'énergie par chocs 


Q>n: 
des à 
dE|dr0>n — asus 


mu? 


2 mv? 


NZ représente le nombre d’électrons par em° du 
es traversé ; z la charge efficace de l’ion inci- 
. dent. 

La différence entre dE/dzx totalet dE/dxo>, donne 
d£EJdx primaire. Dans le cas de l’émulsion nu- 
cléaire Ilford : 


dE z? 
dro5n —10 001 ga 10810 


2 mv? 


dE Jdx est exprimé en keV/u. Dans le cas du ZNe, 
on a représenté les pertes d’énergie en keV/u en 
fonction de l’énergie en MeV (fig. 6). 


FreNC: 


courbe (1) d£Æ dx total, courbe (2) dE/dto>orey 
courbe (3) dE ]dXo>5ske ; 


courbe (4) dÆ/Jdx primaire avec la limite 
 — 2 keW'; 

courbe (5) dÆ/dæ primaire avec la limite 
HD keV: 


On observe sur ces courbes que le d£/dx pri- 
maire diffère peu du d£ /dx total et que l’influence 

du choix de la limite n est peu important. 

Nous avons utilisé ce calcul du dE dx primaire 
pour examiner les possibilités de discrimination 
entre ions lourds de types différents dans les émul- 
sions à faible sensibilité. Différentes grandeurs sont 
liées à cette perte d'énergie : la densité des grains, 
des amas de grains, ou des lacunes ; la distribution 
des longueurs des amas, ou des lacunes. Chacune de 


« . 


ces grandeurs peut théoriquement servir à iden- 
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tifier la nature de la particule incidente. Il convient 


de rechercher celle qui apparaît comme la plus 
satisfaisante. 


b) APPLICATION A LA DISCRIMINATION. — Dans 
une phase préliminaire, nous nous sommes atta- 
chés à rechercher un paramètre facile à évaluer. 
La longueur / de la partie terminale continue nous 
a semblé pouvoir constituer ce paramètre. 

On remarque sur les courbes figure 7. que la 
variation de dE /dx primaire est très rapide, ce qui 
laisse prévoir une localisation assez nette du début 
de cette partie continue. 

En fait, l'examen d’un très grand nombre de 
traces, à l’aide du microscope nous a montré qu’il 
était difficile de trouver un critère reproductible 
pour fixer le début de cette partie continue. La 


HE] 
dx | primaire 


20 40 60 80 100 120 
ioe PE 


longueur des amas croît assez régulièrement, et il 
paraît préférable de choisir comme critère l’appa- 
rition de la première lacune de longueur supé- 
rieure à 0,5 u. Ce critère est sensible au degré de 
développement et au degré d’effacement spontané 
(fading) qui peut devenir important pour ce type 
d’émulsion. Les mesures effectuées sur deux lots de 
plaques exposées et traitées séparément ont con- 
firmé ces hypothèses. 

La première série de ces plaques a donné les 
résultats suivants : 

Ne À — 38,3 p + 6 u représentent la moyenne 
de 150 mesures avec une déviation stan- 
dard 6 — 6 uw. En reportant cette valeur sur la 
courbe (1) de la figure 7, on constate que le début 
de la longueur continue correspond à un dE /dx pri- 
maire de 1 360 keV/u. La base de l’histogramme est 
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très grande, indiquant de grosses fluctuations. Elle 
s'étend de { — 25 [dÆ/dx prim = 1 720 keV/u] à 
| = 55u[d£E/dx prim = 1 050 keV/u] 

Ces valeurs reportées sur les courbes de la figure 7 
permettent les prévisions suivantes : 

Pour ,/N l =14u. 

Pour $C ! — 0 avec des fluctuations de / allant 

jusqu’à ! = 11 p. 

Les mesures expérimentales ont donné les résul- 
tats suivants : 

Pour ,N l —11,4u + 3u: 

Pour $Q ! —0en moyenne. 

La deuxième série de plaques a donné des valeurs 
différentes : 

Ne L — 28,7 up + 5 u ce qui correspond à 
d£Æ [dx primaire de 1 640 keV/u. 

La base de l’histogramme très large elle aussi 
s'étend de { — 19 u [dÆ/dx prim = 1960 keV/u] 
à = 44 u[dEÆE/dx prim = 1 200 keV/u] 

Ces valeurs permettent les prévisions suivantes : 


Pour $0 ! = 14. 
Pour. Nat ==8p; 
Pour $C L —0 avec des fluctuations allant 


IUSqU a 9 Le 
Les mesures expérimentales ont donné les résul- 
tats suivants : 


Pour S0 = M40, 700 
PORN OUEST 
Pour $Q LL —0 avec des fluctuations allant 


JUS A OU 

Pour 44, le dE /dx primaire à 320 MeV, c’est-à- 
dire 80 y de parcours résiduel est de : 2 376 keV/u. 
Quelques traces d’argon trouvé dans ces émulsions 
nous ont permis de vérifier la prévision / > 80 u. 

Ces valeurs sont à nouveau en concordance avec 
les prévisions théoriques, mais différent des pre- 
mières. Ceci peut s'explique par une différence de 
développement et par des variations de sensibilité 
de la plaque d’une part, et d’autre part, par le 
choix du critère qui n’est pas suffisamment objectif. 

Les émulsions À — 1 et À — 2 ont donné des 
résultats peu reproductibles. 

Ces résultats nous ont amené à choisir un para- 
mètre moins difficile à définir : le comptage du 
nombre de lacunes. Nous comptons toutes les 
lacunes, quelles que soient leurs dimensions par 
cellules de différentes longueurs dont l’origine coïn- 
cide avec la fin de la trace dans l’émulsion. Les 
résultats relatifs à l’émulsion X, sont présentés 
sous formes de courbes. 

La courbe (fig. 8) représente le nombre moyen de 
lacunes en fonction de la longueur de la cellule 
pour 4C, 40 et 4 Ne. Les points relatifs à chaque Z 
s’alignent sur 3 droites, qui sont pratiquement 
parallèles. On voit d’après ces courbes qu’il est 
inutile d'augmenter la taille des cellules, la discri- 
mination ne s'améliore pas. Les deux lots de 
plaques que nous avons mentionnés précédemment 
ont donné des résultats comparables. 
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Sur la courbe (fig. 9) nous avons représenté le 
nombre de lacunes en fonction de Z pour des cel- 
lules de taille croissante. La déviation standard 6 
de la distribution du nombre de lacunes est repré- 
sentée sur la courbe. 


nombre 
de 12 
lacunes 


longueurs des 
celluks en 


EMULSION Ko 


nombre de lacunes 


A0 20 30 0 50 60 70 


Fic. 9. 


En extrapolant la droite correspondant à la cel- 
lule de 784, nous voyons que n — 0 pour Z > 16. 
Un certain nombre de traces d’argon Z —18, de 95, 
à 100 u de parcours dans cette émulsion Ko, nous 
ont permis de vérifier cette extrapolation.: 

Cet ensemble de courbes permet à l’expéri- 
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mentateur de discriminer les ions dans ce domaine 
de charge avec une incertitude d’une unité. On voit 
que 23 u de parcours constituent une longueur 
suflisante pour discriminer des ions de Z < 12. Si 
l'expérimentateur a des longueurs de traces supé- 
rieures à sa disposition, il peut identifier des Z de 
valeur plus élevée. 

Des comptages sur des cellules de longueur crois- 
santes par exemple : 23, 39 et 78 4, permettent par 
recoupement de réduire l’incertitude sur Z. 

Les résultats obtenus par la même méthode à 
l’aide d’émulsion À — 2 sont moins satisfaisants. 
La séparation entre les différents Z est moins 
bonne. La diminution de sensibilité de la À — 2 
par rapport à la K, ne se traduit pas forcément par 
une augmentation du nombre de lacunes, mais pro- 
bablement par une augmentation de la longueur 
des lacunes. Ce serait alors ce paramètre qui devrait 
être examiné. 

Quelques comptages de lacunes sur l’émulsion 
K — 1 ont donné des résultats fort peu différents 
de ceux obtenus sur X,. 

Il reste évidemment à évaluer de façon quanti- 
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tative, l'influence du degré de développement sur 
ces critères de discrimination. 


Conclusion. — L'utilisation des émulsions à 
grande sensibilité en vue de la discrimination des 
ions lourds a fait l’objet de nombreux travaux [2], 
[7], [8]. 

Ce travail montre l'intérêt de l’utilisation des 
émulsions de faible sensibilité. Les domaines inté- 
ressants de Z ne sont pas les mêmes dans les deux 
cas. Ce point fera l’objet d’une publication ultérieure. 

Ce travail a été exécuté au Laboratoire de Phy- 
sique Corpusculaire de Strasbourg. Nous tenons à 
remercier vivement M. le P' Cüer pour l'intérêt 
qu’il a bien voulu témoigner à ce travail. Nous 
exprimons notre très vive gratitude à M. le 
P' Barkas de l’Université de Berkeley et à ses 
collaborateurs MM. H. Heckmann et F. Smith qui 
ont assuré l’irradiation des émulsions au Hilac. 
Nous remercions le Commissariat à l'Energie Ato- 
mique pour son aide matérielle. 
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MÉCANIQUE STATISTIQUE D'UN PLASMA LORENTZIEN INHOMOGÈNE ET ANISOTROPE : 


ÉTUDE DE LA DISTRIBUTION ÉLECTRONIQUE 


Par Raymonn JANCEL et THéo KAHAN 
Institut Henri-Poincaré, CG. N. R. $., Paris. 


Résumé. — En s'appuyant sur les résultats de la 1'e partie (1), on donne explicitement le sys- 
tème d’équations permettant de calculer les termes d’ordre zéro et d’ordre un du développement 
de la fonction de distribution. Ces résultats sont étendus au cas d’un champ électrique extérieur 
oscillant, en cherchant des solutions stationnaires par un développement en série de Fourier, super- 
posé au développement en fonctions sphériques. Enfin, on indique les problèmes physiques relevant 


de ces méthodes de calcul. 


Abstract. — Based on the results obtained in the first paper (1), the authors ‘give explicitly 
the system of equations allowing calculation of the terms of zero order, as well as the terms of 
order one, of the expansion of the distribution function. For an external oscillating electric field, 
a stationnary solution is obtained by a Fourier analysis superimposed upon the expansion in 
spherical harmonics. The earlier results are generalized to this case. Physical problems which 
can be dealt with by these methods are indicated. 
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6. Équations générales du premier ordre. 
Nous allons appliquer le développement précédent 
au caleul des termes d’ordre zéro et un du dévelop- 
pement de f et nous utiliserons à cet effet l’équa- 
tion générale (4.1) en nous réservant de montrer 
par la suite que les résultats que nous obtiendrons 
sont également valables pour un champ magné- 
tique d’orientation quelconque. Faisons succes- 
sivement dans (4.1) ! —0, m —=0 ; puis L — 1, 
avec m = + {et 0. Il vient pour! = m —0: 


da 1 : Le dax” 


, : doi : daŸ 
Eur, Len) par(, 08 | Dr, pau, or 2% 
e de 


+ HA, 1) D BK, 0) a] 5 (0). 
(6.1) 


D’après les formules (9) dei l’appendice, on a 
d'autre part : 


Dao 1 2 ee =) ; da da 
y + 5/5 AA if dy )| 


(1) Voir la 1re partie : J. Physique Rad., 1959, 20, 35. 


1 daf JAUNE ES da M: 
PRO NAS EE D a 
de 1/2 rs Vinlr(é . 
: da de) VE doŸ 
— 11, = + — + = PET, ==. 
: 1 "Ve ALLIE 


Eee | à VERTE 2 NS Dre) 


= Jj 1. (20): (6.2) 


fl 
5; pl V3 De a 


On voit apparaître dans les parenthèses les 
expressions figurant dans (5.18), de sorte que l’on 
obtient, par simplification, l’équation correspon- 
dant =D 


de 
dt 


Ÿ re 
= + LP = Jin. (/) 


Te dir D 
+ 3 divr fe FPS 


(6.3. 


qui, mise sous une forme condensée, s’écrit : 


D 
se a es divrfe + NA TIZ (VE Te") = Jrr, (9) 
(6.3’) 
De même, les équations relatives à L — Fe 


m = + 1, 0 s’écrivent : 


dat | do .doÿ 
+ +5l—C0,0 ee KE 2%) | 


(Te —1l;) e C(0, 0) & 


4 
A 


2 Fe 


+ io oi + vod = 0 (6.4a) 
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b) 1—1,m——1 
dit 1 da ous 
7 +32£E(0, 0) ‘ae y 
+ LU + ir) Et, 0 2% 
: 2 LA x y , A 
— à op AT + val — 0 (6.46) 
CEE MS (72 —00 
doQ da$ 1 daû 
ne 410, 0) Su 7; T2 AU, 0) S7- + vlaî— 0. (6.Ac) 


Pour faire apparaître les combinaisons des « 
de (5.18), on remplace les équations (6.4a) et (6.4b) 
par celles obtenues en les ajoutant et en les retran- 
chant de manière convenable. Il vient ainsi : 
SERRES 2 I da0 

x CAN ES 


— i wB(ar + «1) + vi(art— xl) — 0 (6.5a) 


D + on) , ; /2000, ; /2 Ty doù 
dt 3 dy V 
— à œp(art — ai) + vi(ar ? + al) = 0. 


En tenant compte de (5.18), les équations (6.5a), 


(6.5b) 


. (6.58) et (6.4c) prennent la forme suivante : 


RTE SEM ro AU ) (L) k 

J … | 7. NA E wp Jey + Vifez — 0 (6.64) 

à ME) T0 /() @) 

D nn 9 46 
«) (0) (0 

nn (6.6c) 


Ro oz ANA 


Ces trois équations peuvent se mettre sous la 
forme vectorielle condensée : 


r 2/0 


(0) | À 
*& V. 5. 


D 
ni r /€ 


(1) 
PE V = {) 
dt | 1e 


(6.7) 


— (or As) + 


Cette équation vectorielle ayant été obtenue dans 
le cas particulier w$ parallèle à Oz, il nous reste à 
vérifier qu’elle reste valable dans le cas général 
d’une orientation quelconque de wg. 

Il suffit pour cela de remplacer dans les équations 
précédentes les termes en im wg a par l’expres- 
sion (3.38) ; on trouve 


GE) EVE 
0 


(tr op, ©B,) Ver OB, oi (6.8a) 
REMEOM EE 
 — és 
— (op, — ti &8,) TE — L OB, 1 (6.8) 


ON VERRE 
—1 


1 
j œ1 = : A1 ,/; ï 
( op, + 08, Re (op, ui ©B,) RIVE: (6.8c) 
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. On remplace, comme précédemment, les équa- 
tions (6.8a) et (6.8b) par leur somme et leur difté- 
rence et l’on fait apparaître les composantes de f{® 
(d’après (5.18)) pour trouver en définitive : 

on, D — un, AD 
O8, Drop, JDE (on A) cat d.) 
oBy 19 — 8, (1) (6.9) 


Nous aboutissons ainsi au système suivant de 
& équations 


RTE 
EVE a Nr 
4 d 


3 V£dr, PET IE Jr) (6.10) 


de + y, 1 
dt ke 


T {© 
LL Fe ES Af®) un V1 10) — 10} 
e e 


Ce système permet de calculer en principe les 
2 premières approximations de la fonction de dis- 
tribution électronique du plasma considéré. 


7. Résolution de l’équation de Boltzmann en 
présence d’un champ électrique alternatif et d’un 
champ magnétique constant. — Dans le cas où l’on 
a affaire à un champ électrique oscillant que nous 
écrirons F cos ot, l’équation de Boltzmann revêt 
la forme suivante : 


dfe V. T 17% o — J 
Se + Ve Vrfe + [FE cos ot + (wg À We)].Vre fe = (fe) 
(RL) 


où Fest indépendant du temps. Pour résoudre 
cette équation on utilise le développement en har- 
moniques sphériques (2.1) ; dans ce développement 
les «?* sont fonctions du temps et il est naturel, 
étant donné la variation sinusoïdale du champ élec- 
trique, de développer ces a} en série de Fourier ; 
nous poserons donc : 
+Co 


UE 


k=—00 


dr. Ve) FE (7.2) 
où, d’après la relation (2.4) et la condition de réalité 
des coefficients de la série de Fourier, on a : 


Eh (— 1)" an. (A) 


D’après (7.2), la fonction f, revêt la forme sui- 
vante : 


+oo +oo m= +1 


DAS D RAT VOREFe, 0, o) eikot, 


k=—00 1=0 m=—1 


fer, Ve, t) = 


(7.4) 


En portant le développement (7.4) dans (7.1) et 
en identifiant terme à terme ainsi que nous l'avons 
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fait précédemment, on obtient une équation géné- 
rale qui permet de calculer les fonctions 
inconnues «%. Nous allons établir cette équation en 
nous appuyant sur-les résultats obtenus dans le 
calcul des «”. 

D'une manière générale on pourra remplacer 
simplement « par «% toutes les fois où l’opérateur- 
temps n’agit pas sur f; c’est le cas pour les 


— — 
termes V,.Vife, (wB A Ve) Ve fe et J(f). On a 
done immédiatement pour ces termes, d’après 


(2.23), (39), (3.34) et (3.38) : 
COM ECO ET 
HAE Al ) 


DD RTC IE 
k=—0c0 1=0 Mm—=—1 

qui peut s’écrire en vertu de la linéarité de l’opéra- 

teur J : 


J(fe) = 5 J{ao Yo) efrot (7.5) 
k=—0co 
+oo +co +1 : 
M D Sr Jimi 


k=—0c0 l=1 m=—1 


+00 + +co +oo +1 l 
D rien lire OO VONT CRE 


k=—0co k=—00 1=1 m=—1 
puis : 


L# Ÿ D ÿmei 
e-Vrfe — Limit CFNESX 


1 
| 2 | C(1—1,m 1) (ES 


m—1 
: cs) 


é dx dy 
+ D +1, m—1) Vi doRrri — ee) 
dx dy 
Jen m+1 
+ E( 1, m 1) É. Là dARI—1 
dx dy 


7 m+1 


Pl Ad, mA) ps (£e fre Dar 
dx dy 


# 22 


m m 
RATE 7) Le + Ve BI +1, m)  _ 
De même, si w% est parallèle à Oz : 
(os À Ve). Ve fe = ion 2 mod Yi ét 
M 


et si l’orientation de w3 est quelconque 


(op A Ve). Vvefe = 2 Yr citot 
kl,m 


m—1 
À lion, + om) EL VTT à 1) 


Tir Ear 0e. 


TT EE RM ER, 
sn (7.8) 
2 
+iop, m CA 


Reste les deux termes où figure l’opérateur 


ù ; A 
temps ; pour ele résultat est immédiat et on a : 
dfe (nt m M ikot . 
To D kapy Yy erot ; (7.9) 
dt L'ART 
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D'autre part le terme en T' cos wt s'écrit 


Re : 
T cos ot Vre (Z ou YT citut) 2 (7.10) 
M 


Comme l'opérateur T.V est indépendant du 
temps, (7.10) peut s’écrire : 


eiat + CSSS 


co | T'Vr Doi y) È (7.11) 
2 % l,m 


Le crochet de (7.11) a été calculé et figure 
dans (3.22), (3.26) et (3.27) ; on a donc : 


TVre (So Yt)= Wet We (412) 


m ml 
1 (ee Fe A 
ins rl ) 


1 N 
+ V2 daki+1 D(L + 1,m— 1)) 
dVe 


Do 
dVe 


IP 
= a DREEE FUI + 1, m + 1) 
Ve 


RTS r,) ( Ei—1,m+i) (7.13) 


dar dax. 1 
= He AU 1, m) + EEE pq + 4 mi) | 
ù dVe 


e 


wi, = Z (+3) ry ÊE _ Ly p{ + dm — 1) oh 
als il 
2 


(7.14) 

Il reste à effectuer la sommation par rapport à 
Pindice k, ce qui peut s’écrire : 

eik+lot L eilk—1)ot 

DIRE RP ER ET 


ù - (Wx + Wr) ou 


RW: 
D eikot ps de (7.45) 


Win + Ws 
2 ? £ 
Établissons maintenant l'équation générale 
pour w7 parallèle à Oz ; on a, d’après (7.5), (7.6), 


(720,1(7.9);1(74S ) ET EREr (745), l'équation sui- 
vante : 


m—1 m—1 
ko ay + ; - | CES AE ___ _: 


m1 1 
+ D(I +1, m—1) Va pe dt) 


dx 7 
doi ee 
ELEMENT E ( Ê É 
ee po M 1 dy 
m+i m+1 
EE EUR EALArS pese 1 jai 


AE DaR—1 2 dAkU+1 
nr 0) Beuvrus + Va B(I +1, LÉ 
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; Vente 
k ane ct n x 
" ee é dAx+ — 

7. NA 


mM—1 mM—1 
dAR—AU+1 | dARHI,UHI 
L 
dVe dPe 


+ DU+1; m—1) a( 


] 


«a Dai 
dVe 


] 


(7.16) 


m+1 
DAT ,1—1 


Dont, | 
res eut ( 57, 


) 


m+1 
dARHI,1+1 
de 


. Il 


g D +1,m—1) (of ny+1 + ofriur1) 


m+1 
dUR—1+1 


= F4 14, m +1) ra = 
e 


+T [40 — 1, m) Re ee dar 1,11 
| dPe dVe 
DOk—1.4+1 daRe11+1 
dVe | dVe 


+ B(1 +1, m) ra( 
2143 


+7 [ 
ÉDer 
D 


VF( + A, m + 1) (af ua + afiysi) 


+ T, B(I + 1, m) (af is + dau) | 


Jr.xlo%o) si 1 = 0 


+ im op af — 
ou pour 


AT le": 
8. Cas général du développement de f. sous 
… forme tensorielle. — En opérant comme au para- 
graphe 5, on mettra jf, sous la forme tensorielle : 


= 
: YO = Y HU : VO, 
1=0 


fe = 5 aw (8.1) 


Développons maintenant les tenseurs sphé- 
riques Œ et 40 en série de Fourier ; il vient : 


Œ) A eirot 


k=—0c0 


+0co 
eu) > (F$ coskot + GX sinkot). 
k= 0 


Il 


(8.3) 


On définit ainsi un ensemble de tenseurs sphé- 
riques, de degré (21 + 1), À®, 5%, G® associés à 
chaque valeur de k. Le développement (8.1) se 


mettra dès lors sous la forme : 


CO +0co 
fR=Y XYZ eirt(Ap : y) (8.4) 
1=0 k=—0co 
CO CO 1 1 
= YO (59 coskot + GP sinkot) : VW, 
1=0 %=0 


Le tenseur &® étant défini par l’ensemble des 
coefficients 42%, les tenseurs réels F® et GP s’expri- 
meront par des combinaisons linéaires de ces coef- 
ficients. En utilisant la relation (7.3) satisfaite par 
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les «4, on vérifie que le terme général du déve- 
loppement (8.1) s'écrit (avec m Æ 0) : 
aa eikot + ape ikot) Yy 

+ (oi ettot + etat) y," 


Var 


0) —Mm* 
cos kot cos mo(ogy + (— 1)" ay 


+ (tm ont + ot) + 
sin Æct sin mo((— 1)" 7" + (—1)m 077 — of — at) 
— : coskot sin mo(at + (— 1) Œ 


mn CU 0) 


—m*# 


L'isinkot cos mo(oty — (— 1)" ayy 


== + \. 

+ (— 1) Et — ox ) À. (8.5) 

Cette relation détermine 4 combinaisons linéaires 

des CR ay", Cr qui correspondent aux coef- 

ficients réels du développement. D’après (5.10) le 
tenseur V® s'écrit 


pr fe im Var y! m 
VU — (a. AE PE Mean TE Fer, COS Mo, 
2 VIT 
ÉCART Ce 
L 5: L = 1® SIN M0 
21 2+ 


de sorte que #Ÿ et Œ} ont pour composantes, dans 
cette représentation, 


0 0* 
AU À A, » 


FD À (on + (1) TT + (1m of + or), 
ot + (— 1) ET — (1) ani — or). 
| (8.6) 
o%x — opt); 
GP À (tm om + (1) ami — ax ), 
C1) it + (A) om — or — de. 
(8.7) 


9, Équations générales du 1° ordre. — Comme 
précédemment, nous allons maintenant tirer de 
l'équation (7.16), les deux premières approxi- 
mations pour la fonction de distribution électro- 
nique. En s’en tenant aux harmoniques sphériques 
du 47 ordre, la fonction de distribution prend la 


forme (5.12), soit : : 


PV fi, (9.1) 
OO 
+ D (F9 cos k ot + GŸ) sin k oi) 
k=1 
S d 
+ VIFS + D (FŸ cos k ot + GŸ sin k wi) 


k=1 


où les coefficients réels du développement en série 
de Fourier sont mis en évidence, En nous bornant 
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dans la suite à 4 < 1, le développement (9.1) se 
réduit à : 

fe = FU + FŸ-cos st + G9 sin ot (9.2) 
LP FO pe (FT con or GŸ sin ot) 


= f6) + f0 etot + JO ent 
+ Ve + Ve (fi eiot + fi et), 
On obtient par identification : 
AP DS OU Enr, 
À En mn ON à OU D En 
De) = 60, 


où les grandeurs réelles F, G correspondent aux 
tenseurs #® et G® définis en (8.6) et (8.7) à un 
coefficient numérique près défini à partir des O7 
(voir appendice, formules 5). 

On obtient en suivant la methode exposée pré- 
cédemment au paragraphe 5, 


Vr/D, 


(9.3) 


CE — VTT 
(9.4) 


Ro OC n ie à 


et de même 


4 Er * 
SV Ga L agi) 
1 La EE 
fo — se (c0,1 «5 n) (9.5) 
TE 
AVECE 
1” var 
NC 
(D AR NE Fe 
fire SA/ 55 (ed + ai) (9.6) 
CUS 
AVE 
NES 1% 
| Va — 11) 
+ A a t Er + Fe) 
He | 9 V/ 52 GR ua (9.7) 
Auw11E os 
\ VE 


En vertu des relations précédentes, les compo- 
santes du tenseur (8.6) s’écrivent pour 4=— 1,1 —1 : 


&,1 — di,1 Mi + ee ner NE D + 1) 
i— oi + ii + ati — ai") 
NE 
ce 
ali + aji — 2 VE + 19). (9.10) 
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De même les composantes de (8.7) ont pour 
valeur : 


27 
EN à. * (1) (111000 
a Lai ot ee 2 VE (1 — A )e 


A Le 27 4901 « 
— ie + oùt + ali + Li) = 2 V7 PERTE 


T 
dat = 2 EP 0 


Ceci étant, nous allons déduire de l’équation 


générale les équations du 1 ordre, en suivant une 
méthode en tous points analogue à celle du para- 
graphe 6. A partir de l'équation générale (7.16) on 
trouve pour ! = k —=m —=0 


, [pen Ve ee _. (9.12) 
__ FU,1) Và de + aa] + V8 B(1, 0) sea L 
ET 


da 1,1 
dre 


1 
à dVe 


REA — 
. Pee D(1,— 1) (XI: + art) 
IR 


+ I, V8 B(1, 0! ( 


2 


jan 
— 2 F(1, 1) (ol1,1 + o,1) 


+ F,B(1, 0) (ali + 2 | = J11,(08,0), 


soit en utilisant les formules (9.10), (9.11) et (9.3) : 


Leo 


si 57 VAT) = JF). 


(9.13) 


V2 . 
vue 


(°2 


De même pour # = + 1, 1 —0, m —0 on 
trouve deux équations qui s’écrivent : 


ES | 


2 ral ee 1 de ” 


ei | RE 
dVe so 

2 3 
se À qe [Te (œo:i 


nl 
VE F, ai = J1.1 (09.0) 


— 01) — iTy(aoi + aou)] 


(9.41). 


4 


*.5°. (94500 


No 10 


ef 


2 ue : 
— L Dee + Ie va = sr de 


— TG + db1)] + A5 Ta oda = Jr ne no). 
945} 
On ajoute et on soustrait membre à membre 
(9.14) et (9.15), ce qui nous donne compte tenu 
de (9.3) : 


? 


1e 
© ED SV, GT — Ji x (GS) 
ë * 


Fes 
© GP + VrET + 


4 d 


(9.16) 
ANA 


F (PS LUF) = Ji (FO). 


Pour # —0 et ! — 1, on a 3 équations corres- 
pondant à m = 0 et m — + 1 ; on trouve donc : 
AU di lent 0 


7e 0 0 
d%0,0 A(0, 0) / dx—1,0 | 2 : 
A(0, 0) D LT, Ve ( TA le Die 
| (9.17a) 


DR UT ma 


__ C(0, 0) ee Le ee) 
2 dx dy 


À 


da2 1,0 , dXio je 
(GE) Ct0N0) ” (9.476) 


np. Si 
ACTE tem LUE 
DNERURESS RRES | 
E(0, 0) focÿo . De 
2 ee ST 
> Le (Le + il) E(0, 0) er | , i &p 0.1 


= y ao. (9.17c) 


En combinant ces trois équations comme précé- 
demment, on obtient 


va (D () 
NN le ENS de 9.18 
V, F0 LEA Sr. (wB AFo ) V1 F0 (9.18) 
Voyons maintenant le cas où 4 = + Letl —1; 


on à 2 jeux de 3 équations correspondant aux trois 
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Valeurs possibles de m; m —0 et m—+1.Il 
vient ainsi, toujours d’après (7.16) : 
Cl Em UE 
iwo®,1 + + A(0, 0) _ 
1,1 3 , dz 
HET, 0) Do _ = ÿ 591: (9194) 
PA N ole D ME AE 
DPI 
LE) _ . 
er 2 dx dy 
C(0, 0) ; 250,0 : 1 1 
1 P,) —> t ŒI,1 = — Vi U, 
4 Ve (Lx — 11%) di OB 41,1 1 X1,1 
(9.196) 
CRE 
0 EL UN OMLoS con 
LO &1,1 9 ( S et | 
E(0, 0 Oo Tee —1, 
+ : . Tr + il 5  — LB 11 = — Vi di; 
(9.19c) 
et l’on a pour le groupe # = — 1 : 
aa LS 
D UUREET 
OO ES 1 ler 9 Fee 
JP daÿ 0 0 
AO 0j ee 00e) 
orne GE 
DIF TEA 
RS C{(0, 0) ET do 
ere ÿ dx 7 
Ce 1e IT) US RAR CN mn Ci 
(9.206) 
CORSA 
, ue EAU O = 0e OT ) 
ee LONX t 1 D dx Ï dy 
0 
_. (lex +11) eu OBS E  T E 71) 
e 
(9.20c) 


ce qui conduit par des calculs semblables aux pré- 
cédents, aux deux équations suivantes : 


SE Vi Oro GUY GT, (9.21) 

ee T oF0 (1) 

PCM TA UNE — A UE Fi 
Cr 


En réunissant les équations (9.15), (9.16), (9.18) 
et (9.21), on obtient le système d'équations diffé- 
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rentielles suivant, dont la solution fournit Ia pre- 
mière approximation à la fonction de distribution fe 


v2 (nn ù : «) (0) 
LE ypoe A VA F — J Fo), 
NU ir pape Ve PA SEA 
URL 1(1) loire CNE 0) 
o Gi k— V..F Ar x (VE TK) = Ji.1 173 
0 3 7 = de 
Vé Y 
ü) pi 75 VrG LT Te Ji.r Ge 
(0) Poe (1) (1) 
NO (Os AFS) ——v,Fo, (9.22) 
e e 
on ED = V: Ge on N GE y, AD), 
(0) 
NET RTE v A + r em UT RARE OUR 
MÉMOE 


10. Développement en présence d’une force indé- 
pendante du temps. — Nous allons étendre cette 
méthode au cas où le plasma est soumis, en plus 
. d’une force électrique alternative, à une force élec- 
trique K(r) constante dans le temps. Il figurera dès 
lors dans le premier membre de l’équation de 
Boltzmann, un terme supplémentaire additif 
K(r).V ve Je. 

Pour résoudre l’équation de Boltzmann avec ce 
terme, appliquons la méthode précédente du déve- 
loppement en harmoniques sphériques et en série 
de Fourier de f.. La contribution de la force K(r) 
s’écrira avec des notations identiques : 


K(r).Vr,f= 2 Yieitot x 
k,l,m 


1 | Ke —1Ky ( DIR ; 
Le ET C(l—1, m—1) 


AR 
va “OÙ MONT 1) 
e 


RS LUiR, (or à 
EE = Ell—1,m +1) 


, Mm + 1) 


Te DORA A(l—1, m) + pad p( 1 
Î FA dVe rs ; E dVe [l , m 


e 


dar 
V2 2— Fi] +41 
e d7e { Î 


(10.1) 


+ (21 + 3) 


Rae 
PE ot + 4m — 1) 


RE 
— © F( +1, m+1) ain 


A 


e 


+ K, B(I +1, m) us] 


En introduisant la contribution du premier ordre 
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de (10.1) dans (9.22), on obtient le système suivant 
de 12 équations : 


Æ (1) 1 0 pr r® 
3 Vr-fo rer 57. V8 Tri) 


ce T _ (V3 Kr).F6) — J(F0)), 
GE - V,.F® + TA à. (V3 T.F) 
+ sr Kb. EP) = J(RP), 
F9 7 v,G£ _— 5 V8 K(r).GŸ) —— J(G), 
© FT — Vr GP + ns + os A GT = GT, 
(0) 
6 AV AE F- SE 
.. D — AT = — VF? (10.2) 


Les systèmes d’équation (9.22) et (10.2) per- 
mettent d'aborder l’étude théorique de divers 
phénomènes physiques tels que : diffusion, effet 
Hall, conductivité et propagation des ondes électro- 
dynamiques dans un plasma inhomogène et ani- 
sotrope. 

Remarquons enfin que nos formules (mutatis 
mutandis) peuvent trouver des applications dans 
la théorie de la diffusion des neutrons lents dans 
des milieux inhomogènes (réacteurs nucléaires). 


APPENDICE 


Nous utilisons dans le texte les fonctions sphé- 
riques sous la forme : 


PSS ON EE 


&n(l + m) | 
dr 
sin #0 cos Or P; (cos 0)eimp (1) 
avec la relation 
VE = (im. (2) 


Le choix des phases que nous faisons ainsi 
conduit à l'apparition du facteur —1 seulement 
pour les valeurs positives impaires de m. Dans la 
théorie usuelle des fonctions sphériques on choisit 
les fonctions de normalisation en mettant le fac- 
teur (— 1)* dans les facteurs de normalisation. 
Avec les notations de Condon et Shortley [31, nos 
YF s’écrivent 

Yi — VI Of On, (3) 


N°10 
avec 


@" — Of" (C. S.) = (m > 0) 


[am SE FAIT mi d 
(©) 1)" Qu PS ENS RT 5 7) Sess 
1 1) 2 QG misit m0 DE 


P, (cos 60) (4) 


= 2 + A1 (1— m)l . dm 
O7" — que En no — 
L “ ne ere ES Gin 


\ (P, cos 0). 


OÙ — (— 116; ”. 


Donnons à titre d'exemple quelques valeurs 
explicites des O7 


eB— 4/5: ®? — NEC Gt = + 4 /Ésin 0 
0 5 
(©) = 4/5 8 006 0 — sine 9; 


41 f15 
O2 — F/T 0 0 sin 0 ; 


7 
O$ — . (2 cos? 0 — 3 cos 0 sin20). (5) 
En outre 


eimop, (6) 


Dm — V2r 


avec la condition de normalisation 


2T 
Dim) Dim 


à 1) Ap — Ovumt. (7) 


Avec ces notations les formules de récurrence 
pour les fonctions ®* prennent la forme: 


A mu on 
cos 6.0; — Ne F1) 1 ere 3) 


m Um me 
(8) 
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1) (1+ m + 2) 
1) (21 + 3) 


m+1 (L— m) (1 — m—1) 
nie V Gi 1) @i + 1) 


Norte re 
GET V DIRES) 


_ er NN ERNER EE 
rs 21 —1) (2 +1) 


sin 60" __@"tt /U+mi 
ï Oi+1 Ÿ TE 


Nous posons dans le texte pour simplifier l’écri- 
ture: 
e {+1 m) (L + 1 + m) 
ARS V Uri +S 


Pre No U+ m), 


(21—1) (2 +1) 
EDUEm+, 


_ rm 
= RER 


} (21 + 3) 


= UV = mA) (9) 
Den = Een 
}(=m +2) 


_ En SR ; 
LEE ü Qi +1) +83] 


C+rmtrm—i) 
CYRSTEN CIRE 


HO 


Indiquons enfin qu’à l'aide de l’équation diffé- 
rentielle de définition des 7, on peut établir la 
relation suivante : 


1 Le 
i Of" =; VU—m) Um +1) 0"! 
vue ee no ee 000 
m PRESSE AD 
m cotg 0.0 = —-<V{I— m){l + m +1) @ 
1 
2V(I+ m) 1— m +1) ®& (11) 


Manuscrit reçu le 21 mars 1958. 
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DIFFRACTION DE LA LUMIÈRE PAR LES ONDES ULTRASONORES 


Par R. R. AGGARWAL, J. S. O. 


Defence Science Organisation, Ministry of Defence, New Delhi. 


Résumé. — Par une méthode d’approximations successives [1], on a poursuivi le calcul de l’inten- 
sité des franges de diffraction, en supposant que la figure de difiraction ne contient qu’un nombre 
limité de franges. Les résultats sont identiques à ceux que donnent d’autre méthodes plus pénibles. 
Celle-ci a l’avantage de pouvoir être aisément généralisée, la précision s’accroissant par substitution 
en retour. On voit que, lorsque l’angle d’incidence varie l'intensité d’une frange de diffraction 
particulière atteint son maximum maximorum pour l’angle de réflexion de Bragg, et ses maxi- 
mums secondaires pour les angles de Bragg correspondant aux autres ordres. 


Abstract. — The procedure of obtaining closed intensity expressions for [1] the diffraction lines 
on the assumption that the pattern exhibits only a limited number of diffraction lines has been 
elucidated further by applying the method of successive approximation. The identity of most 
of the results obtained by others through laborious methods is established. The method has the 
advantage that it can be generalised easily, the accuracy being increased by ‘ back substitution ?. 
It is seen that as the angle of incidence is varied the intensity of a particular diffraction line attains 
a principal maximum at the corresponding Bragg reflection angle and secondary maxima at the. 
Bragg reflection angles corresponding to the other orders. 
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Introduction. — Lucas et Biquard [2] ont pro- 
posé une théorie d’optique géométrique pour la 
diffraction de la lumière par un milieu agité par des 
ondes ultrasonores. Cette théorie a été développée 
récemment par Nomoto [3]. Brillouin [4] a consi- 
déré le phénomène comme dû à la réflexion de la 
lumière par des couches successives de liquide 
d'indice convenable ; plus récemment [5], il Pa 
considéré comme dû à la propagation de la lumière 
à travers un milieu dont l'indice présente des varia- 
tions périodiques. Cela donnait une équation 
d’onde où l'indice était variable, qu’il transformait 
en une équation du type étudié par Mathieu. Les 
solutions étaient alors données par Brillouin, 
Rytov [6] et, dernièrement, Mertens [7] sous la 
forme de fonctions de Mathieu, tandis qu’Exter- 
mann et Wannier [8] utilisaient le déterminant de 
Hill pour résoudre cette même équation. 

Raman et Nath [9] considéraient le phénomène 
de diffraction comme dû à la transmission à travers 
un milieu stratifié, et exprimaient l’intensité à 
l’aide de fonctions de Bessel. Comme cette approxi- 
mation élémentaire ne pouvait expliquer toutes les 
particularités du phénomène, ils développèrent 
l’idée de l’équation d’onde de Brillouin et arri- 
vérent à une équation aux différences mêlées [10]. 
Nath, plus tard, résolut cette équation [11] et 
obtint une solution sous forme d’un développement 
en une série en général non convergente. 

À côté de la théorie optique géométrique de 
Lucas et Biquard, et de la théorie dynamique 
d’Extermann, certaines autres utilisent la même 
équation fondamentale, et l’auteur [1] a pu, par 
conséquent, établir l'identité de différentes expres- 


sions obtenues pour les intensités des franges de 
diffraction par divers auteurs. Pour cela, 1l a résolu 
l'équation aux différences finies mentionnée précé- 
demment, en admettant que la figure de diffraction 
ne comprenait qu’un nombre limité de franges. 

Le calcul par approximations successives a été 
poursuivi dans sa thèse [13]. 

Par une méthode analogue à celle de Nath [14], 
l’auteur a obtenu les expressions exactes de l’inten- 
sité quand la figure de diffraction ne comprend que 
les franges d’ordre zéro et un, obtenues par ré- 
flexion. Ces résultats ont été depuis publiés par 
Pariseau [15]. La méthode d’approximations 
successives à été étendue au calcul approché des 
intensités des franges du 3€ ordre sous incidence 
oblique, ainsi que de celles des franges de 1re, 2e et 
3e ordre pour des angles d'incidence correspondant 
aux réflexions de Bragg des mêmes ordres. 

La méthode de substitution en retour (back 
substitution) indique les conditions pour lesquelles 
les expressions approchées des intensités sont vala- 
bles. La très grande similitude entre ces résultats 
et ceux de Mertens [7] pour l'incidence rasante 
a été également discutée. 

Le présent article esquisse les résultats qui se 
rattachent à la thèse de l’auteur, mais n'ont pas 
encore été publiés. Les notations utilisées sont 
explicitées dans l’appendice 1. 


2. Expressions très approchées de l’intensité. — 
1er Cas : franges de diffraction du premier ordre. — 
En supposant que pour l’ordre zéro, la lumière tra- 
verse la région des ondes ultrasonores sans aucune 
perte d'intensité et que la figure de diffraction ne 


Cr se 


N° 10 


montre aucune frange d’ordre supérieur à 1, 


_ l'expression, de l'intensité Z,, pour les franges 


. dudie ordre est : 


Sin? [u +00) | 


CI 


Cette même expression a été obtenue par Bril- 
louin [4], Rytov [6], Erwin David [16] et, récem- 
ment, Bathia et Noble [17]. Ceux-ci tirent leur 
résultat de l'intégrale de la théorie de la diffusion. 
Connaïssant l'intensité pour les franges du 1erordre, 
on peut corriger par récurrence celle de l’ordre zéro 
dans l’équation 


dY 
Crea Yes + Vs = 0. 


111 (1) 


et on obtient 


4 Sin? [a — 20) | 4 Sin? [a + 2a) | 


# & A 
L = 1 — p?(1 EX 2)? p2(1 Te 2x)? 
; pa [u—2%] sin?[u + 2) | 
Ta (1 — 2) 1 (1 + 20) 


Le TL Se [u 2x) . 


Sin [u + 2@) Se] Cos se 


2€ 
+ ou — 45 
Sin [a —_ 9) | Sin [a + 29) *] | 
Done ee. 
Peu 


Pour l'incidence rasante, cette expression tend 
vers celle de l'intensité exacte obtenue par 
Nath [14] seulement si p est grand et E petit, condi- 
tions de validité de l’expression (1). De même, 
pour « = 1/2, l'expression (1) se réduit à 

a 
Ii 

L'expression approchée (1) de l'intensité est donc 

fausse pour £ > 2, puisque /,, ne peut dépasser 1. 


Cas n° 2 : Franges du second ordre. — À partir 
des expressions des intensités d’ordre zéro et 1, on 
obtient pour le second ordre [1] 


1 __ Sin? [(1 — x) pé] 
Tr: = tot) 20) k(1 — «) 


ee [u—20 €] sir je 1] | e) 


(1 —20) (3 — 24) 
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identique à l'expression obtenue par Erwin 
David [16] et Noble et Bhatia [17]. Pour l’inci- 
dence rasante, il vient 


rs O0 

< Sin? |[—— 
SH Sin? (b£) | q s (E al 
= 3 ê z - Sin ne) > asile (3) 


Lé maximum pour ,,est 4/9p4, petit seulement 
si p est grand. L'expression approchée (1) n’est 
donc valable que quant p est grand, sinon il] faut la 
corriger par substitution en retour. 

D’autre part, J,, possède un maximum prin- 
cipal pour « — 1, et des maximums secondaires 
pour & — 1/2 et « = 3/2. Pour ces valeurs de «, 
l’expression (2) prend une forme indéterminée. 
Extermann [12] avait rencontré la même difficulté 
pour calculer le déterminant de Hill. Pour la sur- 
monter, il faut reprendre tout le calcul avec la 
valeur particulière de «. Prenons, par exemple, le 
cas où «x — 1: le système d’équations différen- 
tielles devient : 


d Ë 
2 TE +ipd=1 
d Ye 
PRIS d = 0 
Leur résolution donne 
LP 
La = Sin? [eE Ji] (a) 
et 
11 
(RER er Zpi 


RE Sin (5) __ $E Cos (5 F + 9? E2 Sin? (S)]. (5) 


L'expression (4) donnant la première approxi- 
mation de l’intensité n’est valable que si (5) est 
petit, ce qui n’est possible que pour p grand 
et E petit. Donc, les conditions de validité sont plus 
restrictives quand « = 1, ce qui est généralement 
le cas, quant « vaut un nombre quelconque de fois 
1/2. On trouve, par un procédé analogue : 


(l+ 2)@= 1/2 (l+ 2la= 3/2 


Res Ne ee (5 F 
= [12 sn pË Cos à 
Lie coin &)] 6 
+ p? &? Sin (£ (6) 
Les intensités maximums aux différents angles 
satisfont aux conditions. 
(L+o)e=1 À (L+ o)oœ= 172 (La) = 0- (7) 


Cas n° 8 : Franges du troisième ordre. — En 
appliquant la méthode adoptée dans le travail pré- 
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cédent de l’auteur [1], on obtient pour Ÿ, l’équa- 
tion supplémentaire : 


et, en utilisant la valeur déjà obtenue pour d,: 


. pË 
DRE | exp { z{1 | 
ir C0) de2u ee2 


exp | (2 — 2x) pË } 
(5 — 2) (3 — 2) (1 — a) 


exp | (9 — 6) E 
DNS TE) ts NE 20) 
+ (8) 
3(14 — 2x) (1 — x) (3 — 2) 
d’où 
{ es (1 — 22) | 
Î+3 66(3 — 2a)? | (2 — x)? (1 — 2)? 
c Sin? (— 20) | 


2 Sin Sin { (1—&) PE } 
7 (2 — à) (1 — &) (1 — 2a) (5 — 2a) 


pÈ 
(1 — 20) z 


Cos { (3 — 2) | 


3 Le 
2 Sin | (1 = 2e) ee | 


Z Sin | (3 — 22) L 
ci 32 — aj° {1 — 20) (5 — 20) 
Cos | (2 — à) PE } 
3pÈ 
F| 


2 Sin { (1 — 0) pË } Sin } (3 — 2a) 


3(5 — 2)? (1 — x) (2 — à) 
Gos | (5 26 | (9) 


à l'incidence rasante 


_ cë) 
, [sn Sin? (bË) 
A 25 


Sin ? F5) 
DATE 
—= Sin ( Sin (pË) Cos F5) 


se _ Sin (5 Sin FE) Cos (2p6) 


— 7 Sin (bË) Sin () Cos (5) | (10) 


dont la valeurs maximum est 121 /(8 100 96). Donc, 
même pour p — 1, la franges du 3€ ordre a une 


N°10 


intensité négligeable. En recommençant les calculs 
donnés ci-dessus pour des valeurs de &« multiples 
impairs de 1/2, on a : 


= ts (40 (À) 


[s Sin? (pË) + 72 Sin? (©) 


(L+ sla= 3/2 (11) 


1 
(T4 3)a=1 Too pf 


8 Sin? (5) — 36E { Sin (36É) — 9 Sin (pë) } 


+ 96? €] (12) 


Re 1 
BR TEN 


(TR [s Sin? (pË) + 72 Sin? 5 


Soin: FS) + 36€ { Sin (36%) — 9 Sin (66)} 
+ 967€], (13) 

L'expression (9) suggère l'existence d’un maxi- 
mum maximorum pour « = 3/2, et de maximums 
secondaires pour « = 1/2 et «x — 1. Pour avoir une 
idée des grandeurs relatives des intensités pour 
æ —0,x — 1/2, x = 1 et « — 3/2, considérons le 
cas où pË est grand par rapport à 1. Les valeurs 
maximums sont approximativement : 


Wien = APE 100 
(Lslu= 12 = 621576 p* 
(Tslu=1 — 67/36 p° 
(Lsla= 32 = 62/16 p°. 


Aiïnsi la frange du 3° ordre, absente pour l’inci- 
dence rasante, apparait pour une incidence égale à 
l’angle de réflection de Bragg correspondant au 
1er ordre, se renforce pour une incidence égale à 
l’angle de Bragg du 2€ ordre et devient très intense 
pour l’angle de Bragg du 3° ordre. 


3. Expressions de l’intensité données par Bril- 
louin et Mertens..— Brillouin [5], puis Mertens [7] 
utilisant les fonctions de Mathieu, ont obtenu pour 
expression de l’intensité une série en pT2?, con- 
tenant, pour la nitme ordre des termes en pr 2, 
p— +9 ,,,, etc. Si on admet que les figures de 
diffraction ne contiennent pas de franges d’ordre 
supérieur à 1, on peut négliger les termes avec des, 
coefficients pT#,p7 6 ..., Leur expression pour l’in- 
tensité du premier ordre Z,, à l’incidence rasante 
devient : 


LE 11 \ PE 
I = — 2 STE NRE 
+1 sen [(: F2 | =: 


L'expression (1) devient alors 


The du LÉ]. 
e k& 
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Aïnsi la formule de Mertens coïncide avec la for- 
mule approchée si 11/36? est négligeable par rap- 
port à 1. De même, la formule de Mertens donnant 
l'intensité pour le second ordre, devient, en négli- 
geant les termes en p—6, p—8 etc. : 


£ 8 
Q 2 
1 Sin (1 5 à = ee | 
Ii — 
Jp A 


slt)" 


1 3 


qui se réduit à (3) si p est grand. On peut prouver 
de la même façon l’équivalence des expressions de 
l’intensité pour le 3€ ordre. Les analyses de Mertens 
montrent donc que l’expression (1), pour l'intensité, 
n’est pas rigoureusement valable en présence de 
franges de diffraction du second ordre. Il faut alors 
tenir compte des termes de coefficient p—4 dans 
les séries donnant les amplitudes du premier ordre. 
Cette correction peut s'effectuer dans notre cas par 
la méthode de substitution en retour précédemment 
indiquée. On constate également que la formule 
obtenue pour l’intensité par la méthode d’appro- 
ximation, sans aucune correction par substitution, 
est exacte quand les franges de diffraction d’ordre 
supérieur sont négligeables. 
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APPENDICE I 
SYMBOLES ET ABRÉVIATIONS 


à, longueur d’onde de la lumière dans le vide. 
*, longueur d’onde du son dans le liquide. 
v, angle d'incidence mesuré dans le liquide. 
to, Indice de réfraction moyen du liquide. 


U, Variation maximum de l'indice de réfraction 
du liquide due aux ondes sonores. 


L, largeur du faisceau sonore. 


27 ui L 
A À 
A* 
Le o 
12 
RTS 


d, facteur d'amplitude de la frange de dif- 
fraction du nième ordre figurant dans l’équation 
aux différences mêlées de Raman-Nath généralisée. 


Manuscrit reçu le 15 octobre 1958. 


RÉFÉRENCES 


[1] AccarwaL (R. R.), Proc. Ind. Acad. Sc., 1950, 31, 
417. 

[2] Lucas (R.) et BiquarD (P.), Rev. d’Acoustique, 1934, 
3, 198. 


[3] Nomoro (O.), Bull. Kobayashi Inst. of Phys. Res., 
1951, 1, 42 et 189 ; Proc. Phys. Math. Soc. (Japan), 
1942, 24, 380 et 615. 

[4] Brizrouin (L.), C. R. Acad. Sc. Paris, 1914, 158, 
1331 ; Ann. de Physique, 1922, 17, 88. 

[5] Brizzouin (L.), « Actualités » Sc. et Ind., 1933, 59 
(Hermann et Cie, Paris). 

[6] Ryrov (S. M.), « Actualités » Sc. et Ind., 1938, 613 
(Hermann et Cie, Paris). 

[7] MerTens (R.), Sinon Steoin, 1949-1950, 27, 212. 


1936, 9, 520. 

Raman (C. V.) et NatTx (N. S. N.), Proc. Ind. Acad. 
Se., 1935, 2, 406, 413. 
Sc., 1936, 3, 459. 


] 
] 
] 
] 
] ExTermManx (R.), Helo. Phys. Acta, 1937, 10, 185. 
] 
] 
] 
] 
] 


[11] Narx (N. $. N.), Proc. Ind. Acad. Sc., 1936, &, 222. 

[12 

[13] AccarwaL (R. R.), Ph. D. Thesis, Delhi University, 
11955 

[14]-Nara (N. $S. N.), Proc. Ind. Acad. Sc., 1938, 8, 499. 

[15] PariseAu (P.), Proc. Ind. Acad. Sc., 1956, 44, 165. 

[16] Erwin (David), Physik Z., 1937, 38, 587. 

[17] BxaTiA (A. B.) et Nogze (W. J.), Proc. Roy. Soc. 


(London), 1953, 220 À, 356. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


LE SPECTRE D’'ÉLECTRONS DE CONVERSION ÉMIS 
AU COURS DE LA TRANSMUTATION Pa —>2#U 


Par Mlle G&orGETTE ALBOUY, 
Laboratoire de Physique Nucléaire à Orsay 


et M. Manuez VALADARES, 
Laboratoire de l’Aimant-Permanent, G N. R. $S., à Bellevue, Orsay. 


Résumé. — On a étudié par spectrographie magnétique le spectre d'électrons de conversion 
émis pendant la transition ##Pa — 2%#U, De nouvelles raies ont été observées, on a pu améliorer 
la valeur des énergies des différents rayonnements y et en déduire leur nature et leur polarité. 
Les intensités relatives de quelques-uns de ces rayonnements soulèvent des difficultés pour l’accep- 
tation du schéma de niveau actuellement admis. 


Abstract. — The conversion electron spectrum emitted during the ?%*Pa —> ?*U transition has 
been studied by magnetic spectrography. Newlines have been observed, energy values of various y 
rays more precisely determined and their multipolarity deduced. Relative intensities of some of 
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these rays do not agree with the level scheme at present accepted. 


Introduction. — L'étude du spectre d’électrons 
de conversion émis au cours de la transmutation 
233Pa _À, 238U a fait l’objet de plusieurs recherches. 
Après les premiers travaux de Haggstrom [5] et 
Lévy [8] une étude très détaillée a été faite par 
Keller et Cork [6] en employant la méthode de 
focalisation à 1800 dans un champ magnétique. Ces 
auteurs ont étudié le spectre depuis 19 keV jusqu’à 
500 keV ; ils utilisaient comme détecteur la plaque 
photographique et, par conséquent, les intensités 
des raies ont été simplement indiquées dans une 
classification allant depuis « très fort » jusqu’à 
«très faible ». D’après l’énergie des rayonnements y, 
Keller et Cork ont proposé un premier schéma de 
niveaux pour l’Uranium 233 résultant de la trans- 
mutation $ du Protactinium 233. 

Elliott et Underhill [3] ont repris cette étude à 
l’aide d’un spectrographe à lentille et, bien que 
leur appareil n’eût pas le pouvoir séparateur de 
celui de Keller et Cork, ils ont pu donner une pre- 
mière mesure quantitative des intensités des raies 
de conversion à partir de 53 keV, les raies au- 
dessous étant trop faibles pour pouvoir être mesu- 
rées. Par ailleurs, ils ont pu décomposer le spectre 
continu en trois spectres et étudier les coïncidences 
entre les électrons de conversion des différents 
rayonnements y et entre ces électrons et le fond 
continu. Les limites des spectres continus ne per- 
mettaient pas d’accepter le schéma de Keller et 
Cork et ainsi Elliott et Underhill ont proposé. un 
nouveau schéma de niveaux pour l’Uranium 233. 

L'étude du spectre continu et du spectre d’élec- 
trons de conversion a été repris par Brodie [2] utili- 
sant un spectrographe à lentille muni d’un dispo- 


sitif d'accélération des électrons après focalisation ; 
ainsi cet auteur a pu descendre dans l’analyse du 
spectre jusqu’à 6 keV et confirmer, en ce qui con- 
cerne le spectre continu, les résultats de Elliott et 
Underhill. Comme résultat de ses recherches, Brodie 
a proposé un schéma de niveaux pour l’Ura- 
nium 233 qui était celui suggéré par Keller et Cork 
mais inversé. 

Une nouvelle étude. des spectres émis au cours 
de la transmutation 2#5Pa — 238U a été effectuée 
par Ong [10] employant un spectrographe à len- 
tille. Ong a étudié le spectre depuis 10 keV jusqu’à 
600 keV, observé quelques nouvelles raies de con- 
version et effectué une nouvelle mesure soit des 
énergies maxima, soit des intensités des trois 
spectres continus précédemment déduits. En outre, 
Ong a mesuré les photons + avec un spectrographe 
à scintillations et les coïncidences y — e—, B—e— 
et 6e. —e . L'ensemble de ses résultats a conduit 
Ong a adopter le schéma de niveaux proposé par 
Brodie. 

Malgré ces nombreuses recherches un certain 
nombre de points restaient à éclaircir : a) la nature 
et la polarité de quelques rayonnements y ; b) un 
désaccord concernant l'intensité de certaines raies 
de conversion ; c) l'existence d’un rayonnement y 
de 475 keV. C’est dans le but de contribuer à 
éclaircir ces problèmes et, à la fois, d'améliorer la 
précision sur la mesure des énergies, que nous avons 
entrepris cette recherche. 


Spectrographe et résultats expérimentaux. — 
Nous avons étudié le spectre des électrons de con- 
version émis au cours de Ja transmutation 
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KeLzLer ET Cork 
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TABLEAU I 
SPECTRE D’ÉLECTRONS DE CONVERSION 2%%Pa ni CEA 
Bropte ONG Nos RÉSULTATS 
NivEeAU [ÉNERGIE y 
E 7 E I E I 
19,7 ) — — 19,29 m Lxr 40,23 
— 2 — = 292,89 m Mi 28,44 
23,4 10 29,5 = 23,03 m 15% 40,19 
—_ — _ == 23,23 f Mix 28,41 
23,8 se 23,5 ne 24,13 f Mini 28,43 
27,9 — 25,9 — 27,09 m Nr 28,49 
— — _ _ 28,16 f 0! 28,48 
_ — 3225 = = EE st 
35,0 es 34,3 nes 35,04 m Mi 40,22 
36,5 450 35,3 — 35,90 m Mirr 40,20 
— — 36,7 == _ — _— 
_ = 38,8 a 38,98 f Nr 40,25 
— = — = 39,24 Î Nix 40,28 
= ee 40,9 = 0,00 NE Or 40,23 
54,0 21 53,0 20 53,48 F Le 75,23 
04,9 — — —— 04,28 m Lis POS 22 
65,3 17,5 64,3 29 64,80 F Le 86,55 
66,2 de Lx 2 65,55 m Lx 86,49 
70,5 1 69,3 7 69,63 m Mi 75,18 
= _ —_ — 69,96 f Mn 95,14 
— — — — 41,7) Claus KI Lx r. Auger 
— — — = 72,59 t: f: KTx Lx r. Auger 
= ve 73,4 0,3 73,81 f M 75,25 
— — — — 74,79 ae Or 70,41 
— — — — 76,2 END KLx Lxyx | r. Auger 
— — — — JA ba f. KLxr Lys | r. Auger 
— — 80,8 — 80,99 m Mi 86,04 
_ es | 13,5 81,34 f Mir 86,52 
83 10,5 81,8 81,98 m ve 103,73 
eu = 82,61 f le. 103,55 
= 84,8 3 85,03 f Ni 86,47 
Æ = 2x = 86,05 t. Ot 86,37 
— — 88,1 (1) — — r. Auger 
| 0,3 
— — 9276 (1) — — r. Auger 
99,2 2 98,2 3 98,23 f Mi 103,78 
— — 100,4 0,8 (1) — — r. Auger 
4104 — 402,1 — 402,23 Î Nr 103,67 
BE es == — 103,40 tt Or 103,72 
— _ 105,8 (1) — — r. Auger 
— 108,5 | 0,2 (1) —_ — r. Auger 
— — 111,4 (1) — — r. Auger 
156 — 155,6 — — — (K) (271) 
— A 184,4 22 184,1 F K 2997 
198 400 196,1 100 196,0 (OM K AUS 
228 8 224,8 10 224,4 F K 340,0 
238 = 250,4 249,8 f Lis 270,7 
as 2) =. un 253,4 1:46) Litr 270,6 
— — 260,4 — 258,9 Î K 374,5 
Ex Ex 267,0 æ 265,5 rt Mu 270,7 
— — 2794 BE 21 F Lx 299,7 
— — — — 278,8 Î dr 299,7 
ED Le Le Le 282,3 m K 397,9 
292 25 290,6 20 289,6 tour, Ti 311,4 
—— — — —— 290,6 m Lax 314,9 
eo — 29557 = 294,0 m Mi 299,5 
— — — —— 298,0 f Nr 299,4 
— — 300,2 — 299,6 m K 415,2 
Le LE 307,4 4,5 305,9 F Mi 311,4 
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TABLEAU I (suite) 


SPECTRE D'ÉLECTRONS DE CONVERSION 


28pa À 233] 


1 


KELLER et Cork BRODIE ONG Nos RÉSULTATS 
NIvEAU |ÉNERGIE Y 
E I E Jl E 7 E 1 
31429 = — 312,0 455 310,1 m Nr 91179 
en 2,5 2 23 sr ns 311,1 f Or 311,4 
320,2 | 319 dr 319,8 1,8 318,1 F rs 339,9 
— — — —— — 9191 m Lx 340,0 
= 2= = — _ — 32275 Î Lrrx 339,7 
336,0 — __ — 334,9 — 334,3 MARS Ms 339,8 
340,1 = ce 340,0 338,6 m Ni 340,0 
ee 22 #3 2 2i ue 339,6 f Or 339,9 
354,5 de és _ 353,5 0,1 354,0 m Lir (2) 374,9 
= = on a 360,0 Le 357,8 te ne 375,0 
= E PS A 370,8 369,5 tr. Mir 374,7 
= ee £s 22 Les _ 373,2 CAE Nr 374,5 
378,1 = nes es 377,4 0,15 377,1 m Le) 398,0 
— —— — — — —- 380,7 Î Lx 397,0 
— — — — — — 392,8 Î Mix 398,0 
394,6 -— — — 394,8 0,25 394,2 m Lrx (2) 415,1 
_ = = = 2 — 396,5 tps Nu 397,8 
— — — — 398,9 — 398,0 (ET Lairr 415,2 
411,8 ee = se 411,6 0,1 410,3 f Mir 415,5 
Fe es JE 415,8 Le 414,0 (Re Nr 154,3 
===>) 


E, énergie en keV :l ’énergie des rayonnements y est exprimée aussi en keV. 
» ne 


I, intensité. 


1) Ces raies correspondent à des passages Auger XLM ou ÆLN ; le pouvoir séparateur de notre spectrographe 
montre que ces raies sont, en effet, des multiplets, mais leurs très faibles intensités ne permettent pas des mesures pré- 


cises des différentes composantes. 


2) Ces trois raies sont plus larges que les autres raies ; il est possible qu’il s’agisse de trois doublets Lr, Lrr, mais la 
faible intensité des sources ne nous a pas permis d’éclaircir ce point. | 


288Pa — 235U à l’aide de deux spectrographes 
(méthode de focalisation à 1800) dont les rayons de 
courbure maxima étaient de 15 cm et 40 cm. Les 
spectrographes étaient installés dans les entrefers 
des aimants permanents. Le champ dans le spectro- 
graphe à plus petit rayon était de — 300 œærsteds 
et pour le grand spectrographe le champ a varié 
de 35 à 72 œrsteds. 

Le spectrographe de rayon 15 cm a été employé 
surtout dans le but de montrer qu'aucune raie 
d'énergie supérieure à 415 keV n'existe dans ce 
spectre ayant une intensité dix fois plus faible que 
celle des plus faibles raies observées. Cet examen 
était important pour voir s’il existait un rayon- 
nement y de 476 keV (les raies L, M, ... auraient 
des énergies comprises entre 455 et 476 keV). Ce 
point sera discuté ultérieurement (voir schéma de 
désintégration). 

Toutes les valeurs données dans le tableau I sont 
basées sur des mesures effectuées à l’aide du grand 
spectrographe et découlant de l'analyse de 
70 spectres obtenus. 

Les mesures des énergies ont été faites en utili- 


sant comme étalons les principales raies du spectre 
du Th B. + C [12]. 

Le tableau I contient aussi les résultats obtenus 
par Keller et Cork (spectrographe magnétique à 
focalisation à 1800), par Brodie (spectrographe à 
lentille) et par Ong Ping Hok (spectrographe à len- 
lle). L'appareil de Brodie lui a permis, grâce à 
l'emploi comme détecteur des électrons d’un comp- 
teur à paroi très mince, de descendre l’analyse 
jusqu’à 7 keV et celui de Ong jusqu’à 10 keV. 
Comme nous employons du film photographique 
l'analyse devient difficile au-dessous de 15 keV ; 
nous pensons reprendre l’étude de cette partie du 
spectre (au-dessous de 19 keV) en employant la 
méthode de pré- ou post-accélération. 

Nous avons inscrit sur ce tableau les intensités 
données par les différents auteurs pour certaines 
raies, Mais nous nous sommes bornés, en ce qui 
concerne nos résultats, à indiquer l’ordre de gran- 
deur résultant de l’examen visuel. En effet, il est 
difficile lorsqu'on emploie comme détecteur le film 
photographique de comparer les intensités des raies 
dont l’énergie varie depuis 19 jusqu’à 400 keV. En 
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E 53,48 64,80 81,98 keV 


Fic. 1. — Spectre d'électrons de conversion dans la région de 53 à 86 keV. 
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E 69,63 81,98 98,23 1841 196,0 keV 


F16. 2. — Spectre d'électrons de conversion dans la région de 69 à 200 keV. 
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Fic. 3. — Spéctre d'électrons de conversion dans la région de 180 à 230 keV. 
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E 224,4 217,3 289,6 305,9 834,3 keV 


Fic. 4. — Spectre d’électrens de conversion dans la région de 220 à 340 keV. 
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tout cas, pour un certain nombre de cas intéres- 
sants nous nous sommes attachés à faire des déter- 
minations quantitatives des intensités relatives des 
raies spectrales, comme nous le verrons par la suite. 
Remarquons que les intensités de la colonne Keller 
et Cork sont les intensités mesurées par Elliott et 
Underhill. 

L'attribution des raies aux différents niveaux de 
conversion est parfois différente de celle qui avait 
été donnée auparavant ; le pouvoir séparateur du 
spectrographe et l'amélioration sur la précision de 


Li 
y 86 keV 


(AS ASS À ua 


Fi. 5. — Microphotogramme du spectre de raies de con- 
version dans la région de 53 à 66 keV (voir fig. 1). 
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F1G. 6. — Microphotogramme du spectre de raies de con- 
version dans la région de 64 à 86 keV (voir fig. 1). 
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Fic. 7. — Microphotogramme du spectre de raies de con- 
version dans la région de 180 à 230 keV (voir fig. 3). 
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F1G. 8. — Microphotogramme du spectre de raies de con- 
version dans la région de 180 à 230 keV. Cliché surexposé 
(voir fig. 3). 
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la mesure des énergies nous à permis, en effet, de 
faire cette attribution sans ambiguïté. 
Les figures 1 à 4 reproduisent différentes régions 


_ 


L, y 31 
Li 299 vY 
MY 31 keV 


Lu ÿ EL 


À My 239 Ky 415 
Ha AM | | 


Fic. 9. — Microphotogramme du spectre de raies de con- 
version dans la région de 270 à 335 keV (voir fig 4). 
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F1G. 10. — Microphotogramme du spectre de raies de con- 
version dans la région de 300 à 340 keV (voir fig. 4). 


spectrales de 50 à,340 keV. Les figures 5 à 10 repro- 
duisent des microphotogrammes soit des clichés des 
figures 1 à 4, soit des clichés pris avec des temps de 
pose différents, mais des mêmes régions spectrales. 


Nature et polarité des différents rayonnements y. 
— RAYONNEMENT y DE 28,45 KEV. — Nous avons 
observé la conversion de ce rayonnement dans les 
couches M,, My, Mn, N, et O, comme le montre le 
tableau suivant : 


22,89 + 5,55 (Mi) 
23,23 + 5,18 (My) 
24,13 + 4,380 (Nr) 
( 
( 


— 28,44 keV 
28,41 
28,43 » 
28,295) 
28,48 » 


27,05 + 1,84 (Nr) 
28,16 + 0,32 (Or) 


Moyenne — 28,45 keV 


Brodie a observé la conversion de ce rayon- 
nement dans les couches Z, et L,,; peut-être la 
raiede 11,6 keV signalée par cet auteur correspond, 
au moins en partie, à la conversion dans la 
couche Lx. Notre dispositif expérimental ne per- 
mettait pas de descendre à des énergies si basses 
par manque de sensibilité du film photographique 
employé comme détecteur. 
. Sur nos clichés la raie M, est nettement plus 
intense que Ÿ;; qui, à son tour, est légèrement plus 
faible que Mn. Les rapports de conversion dans les 
trois sous-niveaux Â] rappellent la conversion, 
dans la même couche, du rayonnement y de 
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33,12 keV émis au cours de la transmutation 
24Am — 257Np [9]. Il est très probable que le 
rayonnement y de 27,45 keV soit aussi un mélange 
MA + E2, à large prédominance de W1. 


RAYONNEMENT y DE 40,23 KEV. — Le fait que 
nous ayons pu décomposer la raie de conversion 
de 23,5 keV (Keller et Cork) en trois raies (22,89, 
23,03 et 23,23 keV) nous a permis de mettre en 
évidence la conversion dans la couche L;r du 
rayonnement y de 40,23 keV ; en effet, cette raie 
avait été considérée jusqu’à maintenant comme 
étant due exclusivement à la conversion en M, du 
rayonnement y de 28,45 keV. Le tableau suivant 
contient les différentes raies de conversion de ce 
rayonnement. 


m 19,29 + 20,94 (Zxr) — 40,23 keV 
m 23,03 + 17,16 (Lxrx) =D T0) 
m 35,04 + 5,18 (Mr) 4022) 
m 35,90 + 4,30 (Mrrr) 020) 
ï 38,98 + 1,27 (Nr) A) 
Î 39,24 + 1,04 (Nrrr) —40,280) 
LEO O0 EN O0 SN Orne 4025 0 
MOVENNOR EE EE — 40,23 keV 


Remarquons que les raies de 35,04 et 35,90 keV 
avaient été considérées jusqu’à maintenant comme 
des raies de conversion en Z,; et L;, d’un rayon- 
nement y de 57 keV mais soit la détermination plus 
précise des énergies, soit le rapport des intensités 
nous conduisent à classer ces deux raies comme 
résultant de la conversion en W,, et M,,, du rayon- 
nement y de 40 keV. Ainsi malgré qu’il existe un 
passage de 57 keV dans le schéma de niveaux de 
P'Uranium 233, il ne nous semble pas qu’on puisse 
parler d’un rayonnement y de cette énergie; 
d’ailleurs, Browne qui a étudié ces rayonnements y 
de faible énergie à l’aide d’un spectrographe à 
cristal du type Cauchois n’a pas observé non plus 
un rayonnement ayant cette énergie. Or, comme les 
raies de 35,04 et 35,90 keV ont des intensités du 
même ordre de grandeur que les intensités des raies 
de 19,29 et 23,03 keV, s’il y avait deux rayon- 
nements y (40 et 57 keV) les deux devraient être 
visibles par spectrographie cristalline. 

Brodie avait signalé une raie de conversion en L,, 
dont l'intensité serait du même ordre de grandeur 
que celle de la raie de conversion en Z,,. En effet, 
sur un de nos clichés les plus exposés, on voit une 
raie très faible qui pourrait correspondre, d’après 
l'énergie, à la conversion en Z,, mais son intensité 
est certainement vingt ou trente fois plus faible que 
celle correspondant à la conversion en Z,,. Les 
raies de conversion en L,, et Li, ont sensiblement la 
même intensité, mais une détermination précise du 
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rapport des intensités est rendue difficile du fait 
que la raie L;;, est la composante du milieu du 
triplet 22,89, 23,03 et 23,23 keV. 

Nous avons observé encore les doublets Mir et 
Nm les deux raies de chaque doublet ayant sen- 
siblement lamême intensité, nous n’avons pu détec- 
ter la conversion dans le premier sous-niveau, 
c’est-à-dire soit A,, soit MN.. 

D’après ces résultats on classerait le rayon- 
nement y de 40,23 keV comme un quadrupôle 
électrique pur. Pourtant, les expériences d’exci- 
tation coulombienne [9] ont conduit [1] à consi- 
dérer ce rayonnement comme un mélange M1 + E2 
dans les proportions, en photons de 20 à 80 ‘%.. 

La méthode expérimentale que nous employons 
(détermination des rapports en intensité L,/L;, 
Li lLirr, MM, ete.) est très sensible aux faibles 
mélanges de Æ2, mais peu sensible aux faibles 
mélanges de M1 ; en tout cas, il nous semble diffi- 
cile d'admettre que M1 puisse être mélangé dans la 
proportion, en photons de 20 %. Nous reviendrons 
sur ce problème dans la discussion du schéma de 
niveaux. 


RAYONNEMENTS y DE 75,19 KEV, 86,49 KEV ET 
103,7 KE V. — Le classement des raies de conver- 
sion de ces rayonnements est fait dans les tableaux 
suivants : 


F . 53,48 + 94,75 [L1) —:75,93 keV 
m 54,28 00,9 (Lin) 75520) 
m 69,63 + 5,55 (M1) — 75,18 » 
Î 69,96 + 5,18 (Mrr) —.75,14 » 
f 73,81 + 1,44 (Nr) — 75,25 » 
t.f. 74,79 0,32 (Or) — 7511 » 
Moyenner-re 7°". — 75,19 keV 
F 64,80 + 21,75 (Lr)  — 86,55 keV 
m 65,55 + 20,94 (Lin) — 86,49 » 
m 80,99 + 5,55 (Mr) — 86,54 » 
f 81,34 + 5,18 (Mr) — 86,52 » 
f 85,03 + 1,44 (Nr) — 86,47 » 
1. 1.086,05 En 0/321Or) — 86,37 » 
Moyennes tre — 86,49 keV 
m 81,98 + 21,75 (Lr)  — 103,73 keV 
f 82,61 + 20,94 (Zrr) — 103,55 » 
f 98,23 + 5,55 (Mr) — 103,78 » 
f 102,23 + 1,44 (Nr) — 103,67 » 
Br AUs 40 0,932 (01 = 108,720? 
MOYEnDE... "00 — 103,69 keV 
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Ces trois rayonnements y présentent les mêmes 
caractéristiques de conversion ; les deux raies de 
conversion en L, et en Z;, ont des intensités dans 
le rapport 10 : 1 ; pour les deux rayonnements y 
plus intenses (75 et 86 keV) nous avons pu observer 
la conversion en Y,, et le rapport des intensités 
est encore de l’ordre M, : M, = 10 : 1. Ce type de 
conversion est caractéristique des rayon- 
nements M1 purs ; un faible mélange en photons 
de Æ2 (1 %) aurait été parfaitement détectable. 

Rayonnement y de 271 ke V.— Le tableau suivant 
contient les raies de conversion de ce rayon- 
nement que nous avons pu observer : 


f 249,8 + 20,9 (Lrr) — 270,7 keV 
Lou te 253,4 + 11727 (Lu) = 9700 » 
t. 2 02655 Lt 59 {Mi —=27070 

Moyenne ter ere — 270,7 keV 


La conversion de ce rayonnement présente la 
caractéristique d’avoir un coefficient de conver- 
sion pour la couche XÀ inférieur à celui des sous- 
niveaux L. D'autre part, la conversion en Lx est 
supérieure à celle de Z;,.. Ces deux faits conduisent 
à penser que le rayonnement y de 271 keV est un 
rayonnement électrique de polarité 3, 4 ou 5 ou un 
mélange M + E à prédominance électrique. 


RAYONNEMENTS y DE 299,6 KE V, 311,5 KEV ET 
339,9 KEV. — Le classement des raies de conver- 
sion de ces rayonnements est fait dans les tableaux 
suivants : 


E 184,1 


| 


184,1 
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F 184,1 +115,6 (Æ) = A7. 
F 277,9 + 21,8 (Zx)  "—=299,705 
f 278,8 + 20,9 (Lir), —=299,75 
m 294,0 + 055 (Mn 099 
f 298,0 + 1,4 (Nr) —299,% » 
Moyenne .. "#7 —1299/6 ke 
196,0 + 115,6 (Æ)  — 311,6 keV 
289,6 + 21,8 (Lx) —311,4 » 
m 200 6 EM 0 ON PTT) OO 
F 305 0 SE OUT) STE 
m 310,4 + 1,4 (Ny —311,5 0)» 
f 3111 +. 03 0e TTaR 
Moyenne CORrPerCRre "311,5 MKeN 
F 224,4 +115,6 (K) — 340,0 keV 
F 318,1 + 21,8 (Lx) ‘© — 339,9 » 
m STOMIES TD OCT TT) ES 0 0) 
f 322,5 + 19,2 (Lin) — 339,7 » 
MON CRU ER DUR), 686 à 
m 338,6 +. 4,4 (Nr) - — 340,0 » 
f 339,627 0,3 (Orne — 33900) 
Moyenne A ET LE à 1000 JIKEN 
196,0 


temps de pose 5 


Î temps de pose 1 
196,0 keV 


F1G. 11. — Comparaison des intensités des raies de 184,1 et 196,0 keV. 


Les raies de conversion dans la couche X de ces 
trois rayonnements y sous les raies X les plus 
intenses du spectre. Comme le rapport des inten- 
sités de ces trois raies était assez différent d’après 
les auteurs, nous avons cru utile de faire une 
aouvelle détermination. Nous avons fait sur une 
même plaque deux poses avec des temps d’expo- 
sition différents (fig. 11) ; la comparaison, au MiCro- 
densitomètre Vassy, des raies des deux spectres 
permet de déterminer les intensités relatives 


(fig. 12). Les résultats sont inscrits sur le tableau 
suivant : 


INTENSITÉS RELATIVES 


Raïss (keV) 


ELLIOTT BRoD1IE OxG Nos 
RÉSULTATS 
184,1 10 4 22 19 
196,0 100 100 100 100 
224,4 10 8 10 910 


No 40 


Comme on voit, l’accord est entièrement satis- 
faisant pour la raie de 224,4 keV tandis que la 
valeur obtenue pour la raie de 184,1 keV confirme 
le résultat de Ong. Remarquons que le pouvoir 
séparateur de notre appareillage nous a permis de 
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F1G. 12. — Microphotogramme de comparaison des inten- 
sités des raies de 184,1 et 196,0 keV (voir fig. 11). 
La raie B (196,0 keV) correspond à un temps de pose 
unité et les raies À (184,1 keV) et C (196,0 keV) à un 
temps de pose cinq fois plus long. 


séparer entièrement les raies de 184,1 et 196,0 keV 
(voir fig. 3) tandis que les auteurs qui nous ont 
précédés étaient obligés de décomposer l’ensemble 
des deux raies. L’intensité de la raie de 184,1 keV 
(conversion ÆÀ du rayonnement y de 229,6 keV) 
joue un rôle important dans l’établissement du 
schéma de niveaux. 

En ce qui concerne la conversion dans la couche ZL 
nous avons pu observer, pour la première fois, la 
conversion dans le sous-niveau Z;, de ces trois 
rayonnements y ; pour les trois rayonnements le 
rapport des intensités L, : L,, est de l’ordre de 10. 
Les raies de conversion en L;;, des rayonnements y 
de 299 et 311 keV sont cachées par des raies intenses 
et la raie de conversion, dans le même sous-niveau, 
du rayonnement y de 339,9 keV a une intensité très 
inférieure à Z,.. 
= Ong a classé, d’après le rapport X/L, les rayon- 
nements y de 229,6 (K/L — 6,3) et 339,9 keV 
(XJL = 5,5) comme M1 purs et le rayonnement 
de 311,5 keV comme un mélange M1 + ÆE2 
(KIJL = 4,0). Le fait que nous employons comme 
détecteur le film photographique ne nous permet 
pas de déterminer avec précision le rapport K/L 
étant donné la différence des énergies de ces deux 
raies, mais nous avons pu déterminer le rapport 
entre les raies L, de ces trois rayonnements de la 
même façon que nous l’avons fait pour les raies X. 
Nous avons trouvé que le rapport des intensités 
des trois raies Z,; est le même que pour les trois 
raies As LL (2996 keV):L; (y 311,5 keV): 
L, (y 339,9 keV) = 20 : 100 : 9. Il est, par consé- 
quent, légitime de conclure que le rapport K/L est 
le même pour les trois rayonnements y. Par ailleurs, 
les caractéristiques de la conversion de ces rayon- 
nements y dans la couche Z nous conduisent à 
classer ces trois rayonnements comme dipôles 
magnétiques purs. 


RAYONNEMENTS y DE 374,5 KEV, 397,9 KEV ET 
415,2 KE V.— Le classement des raies de conversion 
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de ces rayonnements est fait dans les tableaux 
Suivants : 


m 258,9 + 115,6 (Æ) — 374,5 keV 
m 354,0 + 20,9 (Or) SRE) » 
Cu 357,8 + 19,2 (Lrrr) — 375,0 » 
Ho 369,5 5,2 (Mir) — 37457) 
LR 0 75 2 EPS EN TT) 7 OR) 
Moyenne cree — 374,7 keV 
m 282,3 + 115,6 (Æ) — 397,9 keV 
m 379,1 + 20,9 (Zrr]) — 398,0 » 
f 80 ET OL) 097290) 
Î 392,8 + 5,2 (Mr) — 398,0 » 
oies OO EN TE NEE) e =YES, à 
Moyenne rte te ter — 1997 J1KeN 
m 299,6 + 115,6 (Æ) = 415,2 keV 
m 394,2 +. 20,9 (Zrr) = 415,4 » 
LS O8 OC 2 ln) 2) 
Î ADS EN OO) 
RAT OEUVRES 3) 
— 415,3 keV 


MOVENNEREE EE EEE 


ai| 


Les raies de conversion en X de ces trois rayon- 
nements y ont des intensités dans l’ordre suivant : 
K (415 keV) > K (398 keV) > K (374 keV). Les 
trois raies de conversion en L,, ont des intensités 
disposées suivant le même ordre ; on peut dire que 
pour ces trois rayonnements y le rapport X/L est 
sensiblement le même et certainement supérieur à 2 
D'autre part, la raie de conversion en Z;, est plus 
intense que celle en L;r, le rapport étant de l’ordre 
de 3. D’après ces rapports, ces trois rayonnements y 
pourraient être des £2, mais il faudrait des mesures 
d'intensité plus précises pour être certain de cette 
attribution. 

Remarquons finalement que la raie de 357,8 keV 
que nous avons admise comme correspondant à la 
conversion dans la couche L;,,, du rayonnernent de 
374,7 keV avait été interprétée auparavant comme 
la conversion dans la couche ÆÀ d’un rayon- 
nement y de 475 keV. Nous ne croyons pas que cette 
interprétation soit à retenir, soit parce que l’exa- 
men comparatif de la conversion des trois rayon- 
nements y (374,7 keV ; 397,9 keV et 415,3 keV) 
justifie l’attribution de cette raie à la conversion 
en Lin, Soit parce que nous n’avons pu observer 
aucune autre raie de conversion dans les ni- 
veaux L, M, ... d'un rayonnement y de 475 keV, 
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malgré les recherches que nous avons faites spécia- 
lement dans ce but. 

Dans le tableau II nous avons comparé les éner- 
gies des rayonnements y obtenues par les différents 
auteurs. Toutes les déterminations ont été faites, 
comme nous l’avons dit, par spectrographie des 
électrons de conversion, à l'exception de celle de 
Browne qui a employé un spectrographe à cristal 
du type Cauchois. 


TABLEAU II 


RAYONNEMENTS ÿ DE LA TRANSMUPATION Pa —> 235U 
(énergie en keV) 


BROWwNE KELLER BroDIE ONG PRÉSENT 
ET CORK TRAVAIL 
— — 17e 4557 — 
28,67 28,9 29 27,5 28,45 
40,47 40,6 41 20,1 40,23 
ss 58,1 58 57,6 — 
75,4 SE) 76 74,8 75,19 
87,0 87,1 87 86,2 86,49 
— 104,5 104 103,7 110827 
— 272,6 D 2 22150 270,7 
— 301,5 301 300,9 299,6 
= 313,1 313 312,6 311,5 
— 342,0 342 340,8 339,9 
= 376,5 — 375,9 374,7 
— 399,9 400 400,3 397,9 
— 416,4 417 416,7 HASES 
— — — (475,6) — 


Pour les faibles énergies, les valeurs que nous 
avons obtenues sont inférieures à celles de Browne 
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d'environ 0,2 keV. Il semble pourtant que les 
mesures de Browne étaient systématiquement éle- 
vées ; ainsi, par exemple, tandis que Browne a 
obtenu pour le rayonnement y le plus intense de la 
transmutation 241Am — ?23%7Np la valeur 59,78 keV, 
les mesures beaucoup plus précises de Day ont 
conduit à la valeur 59,568 + 0,017 keV. On voit 
que la différence est du même ordre de grandeur 
que celle que nous trouvons. 

Pour les énergies supérieures à 270 keV nos 
valeurs sont systématiquement inférieures ; on doit 
remarquer que la différence provient en partie de la 
mesure des énergies des raies d’électrons de con- 
version, mais surtout du fait que nous avons classé 
d’une façon différente les raies. Ainsi, par exemple, 
tandis que Ong considérait que la conversion du 
rayonnement de 416,7 keV s’effectuait dans les 
niveaux X, L,, M, et N,, nous avons admis que 
cette conversion se réalisait en X, Ly, Lim Mu et 
Ny;ilest évident que dans ces circonstances, même 
si l'énergie des raies de conversion était la même, 
l'énergie déduite pour le photon y serait diffé- 
rente. 


Schéma de niveaux. — Comme résultat de leurs 
recherches sur le spectre d’électrons de conversion, 
Keller et Cork ont proposé pour l'Uranium 233, 
résultant de la transmutation $ du proctactinium, 
le schéma de niveaux représenté dans la figure 13a. 

Ce schéma traduisait parfaitement les relations 
numériques simples qu’on pouvait établir entre les 
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énergies des différents rayonnements y ; évidem- 
ment ce schéma ne tenait pas compte des intensités 
relatives des différents rayonnements y, qui 
n'avaient été appréciés que visuellement, ni des 
spectres continus partiaux jusqu'alors inconnus. 

Les déterminations par Elliott et Underhill con- 
cernant le spectre continu venaient rendre ina- 
ceptable ce schéma ; en effet, ces auteurs ont pu 
décomposer le spectre continu en trois spectres 
ayant leurs maxima à 565,256 et — 140 keV. Basés 
sur cette analyse et sur les intensités relatives des 
rayonnements y, Elliott et Underhill ont proposé 
un nouveau schéma de niveaux représenté dans la 
figure 13 b. Ce schéma est en partie analogue à 
celui de Keller et Cork ; en effet, les niveaux du 
schéma a placés à 0, 17, 75 et 104 keV se retrou- 
vent maintenant en haut du schéma b respecti- 
vement à 416, 399, 342 et 313 keV ; la différence 
entre les deux schémas réside dans le fait que le 
niveau de 376 keV (schéma a) occupe dans le 
schéma b la position à 301 keV. Cette différence 
provient essentiellement du fait que Elliott et 
Underhill n'avaient pas observé les rayonnements y 
de 272 et 376 keV correspondant sur le schéma & 
aux transitions 376 — 104 et 376 — 0. En effet, 
sur le schéma de Elliott et Underhill, il n’y a pas 
de place pour les transitions correspondant aux y 
de 272 et 376 keV. 

Lorsque Brodie a affectué ses recherches il s’est 
attaché particulièrement à vérifier si ces rayon- 
nements y existaient ou non dans la transmu- 
tation 2%5Pa — 288U, Il a pu confirmer l'existence 
du rayonnement y de 272 keV et ce fait l’a conduit 
à déplacer le niveau de 301 keV (schéma b) vers 
40 keV (schéma c). Comme il est aisé de le voir ce 
schéma proposé par Brodie est exactement le 
schéma de Keller et Cork, mais inversé ; le rayon- 
nements y de 272 keV et 376 keV y trouvent, par 
conséquent, leur place et le rayonnement y de 
301 keV qui, dans le schéma b, correspondait à une 
transition directe au niveau fondamental, devient 
maintenant le rayonnement émis dans la transition 
entre les niveaux 341 —- 40 keV. C’est ce schéma de 
Brodie qui a été adopté par Ong. 

Les recherches par excitation coulombienne du 
noyau d’Uranium ?#3, effectuées par Newton [9]ont 
apporté une confirmation à l'existence de ce 
niveau à 40 keV au-dessus du niveau fondamental. 
En effet, Newton a observé trois rayonnements y 
de 40,4 EE 0,2 keV, 215 Es 0,5 keV et 929 Et il keV 
correspondant au schéma de la figure 13 d. Des 
probabilités de transitions observées pour les dif- 
férents rayonnements y, on déduit que les rayon- 
nements de 51 keV et 40 keV sont des mélanges 
M1 +E2 dans la proportion, en photons, de 
20 %, M1 + 80 % E2. 

Le schéma c sur lequel nous avons marqué les 
valeurs de spin proposées par Newton [9] est en 
accord : 
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1) avec les énergies des divers rayonnements y 
observés, 

2) avec la polarité des y de 28, 75, 86, 104, 300, 
310 et 375 keV. 

Par contre, il ne permet d’expliquer ni les pola- 
rités que nous proposons pour les rayonnements y 
de 270 et 416 keV, ni la non observation du rayon- 
nement y de 58 keV. En effet, l'intensité de ce 
rayonnement d’après le schéma (c) devrait être 
environ la moitié de celle de 86 keV, or le fait qu'il 
n'ait pu être décelé sur nos clichés montre que son 
intensité est au moins 20 fois plus faible que celle 
du rayonnement de 86 keV. 

De plus, il est évident que l’acceptation du 
schéma € implique que le nombre de transitions 
de 40 — 0 ne soit pas inférieur au nombre de tran- 
sitions qui, à partir des niveaux supérieurs, abou- 
tissent au niveau de 40 keV. Ces transitions sont 
celles de 272 keV, 300 keV et 376 keV. La proba- 
bilité des transitions de 272 et 376 keV est faible 
par rapport à celle de 300 keV et, on peut, en consé- 
quence, ramener la discussion aux deux rayon- 
nements de 40 keV et 300 keV. On remarquera que 
dans les deux schémas bet c ces deux rayonnements 
sont en cascade, mais, tandis que dans le schéma b 
la cascade s'effectue dans l’ordre de 40 keV (tran- 
sition 342 — 301) 301 (transition 301 — 0), dans 
le schéma c l’ordre est inversé, le rayonnement 
de 301 est émis.d’abord (transition 341 —40) et le 
rayonnement de 40 vient ensuite. En résumé, 
tandis que dans le schéma b c’est le rayonnement 
de 301 keV qui doit avoir une intensité égale ou 
supérieure à celle de 40 keV, dans le schéma c c’est 
l'inverse qui est vrai. 

Malheureusement, il manque des mesures pré- 
cises des intensités de ces deux rayonnements. Le 
seul auteur qui a mesuré les deux rayonnements 
c’est Brodie et il a trouvé que les deux avaient la 
même intensité ; or, les mesures d’Ong, que nous 
avons confirmées, ont montré que le rayonnement 
de 301 keV est cinq fois plus intense (voir tableau) 
qu’on ne le croyait d’après les mesures de Brodie. 
Pour accepter le schéma c il faudrait, par consé- 
quent, conclure que la mesure de l'intensité du 
rayonnement de 40 keV, faite par Brodie, était 
entachée d’une erreur du même ordre de grandeur. 
Il nous a paru, de ce fait, intéressant d’essayer de 
comparer les intensités des raies de conversion des 
rayonnements de 40 et 300 keV, malgré les difi- 
cultés que comporte une telle comparaison lorsqu'on 
emploie la plaque photographique comme moyen 
de détection. 

Pour effectuer cette détermination, nous avons 
pris comme spectre de comparaison, calibré en 
intensités, le spectre du dépôt actif du thoron. 
Nous avons admis les intensités déterminées soit 
par Flammersfeld [4], soit par Krisyouk et al. [7] 
et par Zhernovoy et al. [13] ; les différences entre 
les résultats de ces auteurs sont d’ailleurs faibles. 
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Pour les raies de conversion du rayonnement de 
40 keV du Pa — U, nous avons employé comme 
étalons les raies Aa, B et Bb du dre du 
Th C — Th C” (rayonnement y de 40 keV) et pour 
les raies de conversion du rayonnement de 300 keV 
(Pa — U), les raies H et Z du spectre Th B—Th C. 
On activait une source de ?**Pa dans une atmos- 
phère de thoron et on obtenait simultanément 
pour les deux spectres toute l’étendue de la gamme 
des énergies à étudier, ce qui permettait de négliger 
l'intensité de la source, sa décroissance et le temps 
de pose. 

Nous avons trouvé que le nombre de transitions 
de 40 keV était 3 à 4 fois plus faible que le nombre 
de transitions de 300 keV. Au cas où cette con- 
clusion serait confirmée, la cascade 300,40 corres- 
pondant aux transitions 340 — 40 — 0 devien- 
drait inacceptable et on serait alors tenté 
d'admettre un niveau à 300 keV, comme sur le 
schéma d’Elliott et Underhill. 

Il ne nous semble pas que dans les conditions 
présentes d'incertitude sur le schéma de niveaux à 
adopter on puisse faire avancer la discussion sur les 
caractéristiques des différents niveaux. 


Conclusions. — L'étude du spectre d'électrons de 
conversion émis au cours de la transmutation 
233Pa — 238Ù nous à permis de : 


a) montrer que certaines raies de conversion 
sont, en effet, des doublets ou même des triplets et 
d'observer dix-neuf raies nouvelles ; 


b) augmenter la précision sur la mesure des 
énergies et, par conséquent, rendre plus précise 
l'attribution des raies aux différents niveaux de 
conversion ; ce classement nous a conduits à éli- 
miner le rayonnement y de 57 keV jusqu'à main- 
tenant admis et nous avons montré qu’il n’y a pas 
de raison expérimentale pour admettre, non plus 
un rayonnement y de 476 keV ; 


c) établir sans ambiguïté la nature et la polarité 
des rayonnements y de 28, 40, 75, 86, 104, 300 et 
340 keV ; il reste à préciser la nature et la polarité 
des rayonnements y de 271, 375, 398 et 415 keV ; 


d) montrer que nos mesures — malheureusement 
peu précises — des intensités relatives des rayon- 
nements de 40 et de 300 keV rendent difficile 
l'acceptation du schéma de niveaux actuellement 


admis pour l’Uranium 233 descendant de la trans-. 


mutation 8 du Protactinium 233. 

Nous remercions bien vivement M. Vernois et 
Mne Pluchet qui ont bien voulu préparer les nom- 
breuses sources employées dans ce travail. 


Manuscrit reçu le 23 mars 1959. 
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EFFETS DES CHAMPS ÉLECTRIQUES 
ALTERNATIFS ET CONTINUS 
SUR LA LUMINESCENCE 
DE CERTAINS SULFURES DE ZINC EXCITÉS 
PAR UN RAYONNEMENT «. 


Par A. CocHE et R. HENGK, 


Laboratoire de Chimie Nucléaire, 
Faculté des Sciences de Strasbourg. 


Le renforcement par un champ électrique de la lumui- 
nescence de certains sulfures de zinc a été observé par 
plusieurs auteurs dans différentes conditions : avec des 
champs alternatifs, sur des cellules où le scintillateur, 
excité par un rayonnement X, est dispersé dans un 
diélectrique (Destriau [1]), avec des champs continus 
sur des cellules constituées par des films minces de 
scintillateur obtenus par vaporisation et excités par un 
rayonnement ultra-violet, (Cusano [2]) ; Mattler [3] a 
signalé un renforcement de la luminescence provoquée 
par irradiation, aux rayons &, d’un sulfure de zinc et 
de cadmium, activé au manganèse.et à l’argent. 

Nous avons examiné l'influence de champs élec- 
triques alternatifs et continus sur la luminescence de 
sulfures mixtes de zinc et de cadmium (1), activés au 
manganèse en proportion variable (A : 10-? g Mn/g 
SZnCd, B : 10-% g Mn/g SZnCd, C : 104 g Mn/g 
SZnCd) et excités par les particules &« de ?1Po. Les 
cellules sont constituées par une couche (épaisseur 
environ 100 u) de scintillateur dispersé dans un poids 
égal d’araldite, déposée sur une lame de verre conduc- 
teur. La deuxième électrode est formée d’une fine pelli- 
cule d’argent (obtenue par vaporisation), à travers 
laquelle se fait Pirradiation &. Dans toutes les expé- 
riences, la luminance est mesurée par l’intermédiaire 
du courant d’anode du photomultiplicateur 53 AVP 
sur la cathode duquel repose la cellule. 

L'application d’une même tension alternative, infé- 
rieure au seuil d’électroluminescence, à diverses cel- 
lules excitées par des rayons «, provoque une variation 
de l'intensité de la lumière émise, différente suivant les 
scintillateurs utilisés, comme le montre la figure 1. Sur 
cette dernière, on a porté, en fonction du temps #, 
pour trois cellules réalisées respectivement avec les 
sincillateurs À, B et C, le rapport de la luminance L,,y 
(sous l’action simultanée des rayons & et de la ten- 
sion V) à la luminance Z, (sous la seule action des 
rayons a). Le rapport (L,+v)/L, représente le coeffi- 
cient d'amplification © tant que la contribution de 
l’électroluminescence Ly est nulle. Dans le domaine de 
tensions supérieures au seuil d’électroluminescence (de 
l’ordre de 200 volts), © est défini par la relation : 

BE. (Lo+v F2 Ly) [La é 

Le coefficient d'amplification n’est supérieur à 1 que 
pour les scintillateurs A et B, le scintillateur C donnant 
un effet extincteur. 


(1) Fournis’par les Établissements Massiot-Fluor. 


La variation de p en fonction de la tension V est 
représentée sur la figure 2. Pour le scintillateur B, on a 
étudié séparément les effets produits dans les deux 
bandes d'émission de ce scintillateur, en interposant, 
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entre la cellule et la photocathode, des filtres appro- 
priés. On constate que le champ électrique provoque un 
effet renforçateur prononcé, augmentant constamment 
avec la tension, dans la bande d’émission de plus 
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grande longueur d’onde (bande jaune) et au contraire 
un effet extincteur dans la bande bleue. Gobrecht et 
Gumlich ont observé un comportement analogue dans 
le cas d’une excitation par des rayons X [4] et des 
rayons cathodiques [5]. Pour le scintillateur A, qui ne 
possède qu’une seule bande d’émission, le coefficient 
d'amplification p tend vers une valeur constante pour 
les tensions élevées. Ce résultat est à rapprocher de 
celui obtenu en irradiant le scintillateur par des 
rayons X [1]. 

En appliquant aux scintillateurs À, B et C des ten- 
sions continues, on constate également des modifi- 
cations durables de la luminescence excitée par un 
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rayonnement «. La figure 3 montre, par exemple, en 
fonction du temps {, la variation de l'intensité de la 
lumière émise, consécutive à l’application d’une ten- 
sion continue V à la cellule A. 

La diminution continuelle du renforcement observé 
provient de la formation progressive d’un champ de 
polarisation [6] réduisant la valeur du champ efficace. 
A la coupure de la tension, le champ de polarisation 
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formé subsiste et agit de la même manière qu’un 
champ extérieur, produisant une nouvelle stimulation 
de la luminescence. 

Les variations de l'intensité de la luminescence 
excitée par des rayons «, obtenues par application de 
tensions continues, sont qualitativement les mêmes 
que celles observées avec des tensions alternatives. 
Ainsi le scintillateur À donne toujours un effet de sti- 
mulation, tandis que le scintillateur C ne donne qu'un 
effet extincteur. 

Destriau [7] a signalé récemment un effet de mémoire 
lié à l'effet renforçateur produit par des champs alter- 
natifs, de la luminescence excitée par les rayons X. 
Nous avons trouvé un effet analogue en irradiant 
par des rayons « les scintillateurs À et B, soumis à des 
champs électriques tant alternatifs que continus, 
comme le fait ressortir la figure 4. Le scintillateur est 
sensibilisé par action simultanée d’un rayonnement 
et d’un champ électrique, pendant un intervalle de 
temps t, — t,. À l'instant {, on supprime l’excitation & 
et on coupe la tension. Si on maintient le scintillateur 
dans l’obscurité pendant un intervalle de temps #5 — t, 
(qui peut dépasser plusieurs heures) et si on l’excite 
ensuite à nouveau avec le même faisceau de rayons « 
qu’initialement, la luminance prend des valeurs supé- 
rieures pendant plusieurs dizaines de minutes à la 
valeur initiale Z,. L'effet est plus marqué avec une 
tension continue (courbe I) qu'avec une tension alter- 
native de valeur; efficace égale à celle de la tension 
continue (courbe IT). 

L'application d’une tension continue à un scintil- 
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lateur À ou B, en l’absence de tout rayonnement œ, 
provoque également une sensibilisation de celui-ci. 
La figure 5 montre les résultats obtenus avec le scintil- 
lateur B (dans la bande d’émission jaune). La tension 
continue appliquée au scintillateur pendant un inter- 
valle de temps {, — {, ne produit aucune émission de 
lumière par électroluminescence, Mais si l’on soumet 
le scintillateur à une excitation «, quelques instants 
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après la coupure de la tension, l’intensité de la lumière 
émise passe par un maximum, dont la valeur dépend de 
celle de la tension appliquée précédemment, de la durée 
de son application et de l’intervalle de temps entre la 
suppression de la tension et le début de l’irradiation «. 

Des résultats d'ensemble sur ces différents effets 
seront publiés ultérieurement. 


Lettre reçue le 15 juin 1959. 
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EFFETS MAGNÉTIQUES 
D'IMPURETÉS DE TRANSITION 
DANS L’ALUMINIUM 


Par M. Michael Alan TAyLor, 
Cavendish Laboratory, Cambridge, 


et MM. Jean-Paul BurGEr et Jules WUucHER, 
Laboratoire Pierre-Weiss, Institut de Physique, Strasbourg. 


Le manganèse, le chrome et le vanadium forment 
avec l’aluminium des solutions solides, à la tempé- 
rature de l’eutectique, jusqu’à de faibles concen- 


15 AX 
OS 0323% at.Cr 
OS ee 


TIR 
& 100 200 300 
105 O - OC,078 Ve at V 
-10 
O O192 % at.V 
EN AE 


trations de 0,67, 0,37 et 0,20 % atomique respecti- 
vement [1]. Les divers alliages ont été préparés au four 
à induction dans le vide à partir d’aluminium 99,99 
de manganèse 99,99, de chrome 99,90 et de vana- 
dium 99,50. Maintenus pendant 12 heures quelques 
degrés au-dessus du point de fusion, puis recuits 
pendant 12 heures à la température de l’eutectique, 
les divers alliages ont été trempés dans l’eau. 
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En vue des mesures magnétiques, effectuées à la 
balance de translation Foëx-Forrer, les culots 
d’alliages ont été tournés sous forme de cylindres de 


- dimensions et poids sensiblement égaux, ceci pour 


ue une précision relative aussi grande que pos- 
sible. 

Les résultats de l’étude thermomagnétique effec- 
tuée entre 90 0K et 290 ©K sont rassemblés dans le 
tableau suivant. 

Mesurés dans trois champs, les divers échantillons 
préalablement lavés dans NaOH et HCI, puis séchés 
dans l’éther, n’ont montré à température ambiante 
aucune dépendance de y avec le champ. Dans la 
figure 1 nous avons porté en fonction de la tempé- 
rature la valeur Ay = 7m — za (où xm est la suscep- 
tibilité moléculaire de l’alliage et y4 la susceptibilité 
atomique de l’aluminium pur) dans le cas du vana- 
dium et du chrome. L’addition de vanadium fait 
décroître le paramagnétisme de l’aluminium d’une 
quantité indépendante de la température entre 
100 ©K et 290 oK. Rapportée au vanadium seul, 
cette décroissance correspondrait à un y diama- 
gnétique de — 950. On pourrait attribuer cet effet 
à une valeur élevée du diamagnétisme de Landau. 
Cet effet est à rapprocher du fort diamagnétisme 
mesuré à température ambiante par Auer [2] dans le 
cas du Li (— 293) Cu (— 146) Ag (— 93) Mg (— 155) 
Zn (— 156) en solutions diluées dans l’aluminium. 

Vogt [3] a signalé que les alliages dilués AlMn ont 
une forte valeur de Ay, indépendante de la tempé- 
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rature. Dans la figure 2 nous avons porté les Ay 
déduits de nos mesures pour quatre concentrations 
de manganèse. En première approximation ces Ay 
sont constantes avec la température et leurs valeurs 
moyennes portées en fonction de la concentration de 
l'élément ajouté (fig. 3) sont en excellent accord avec 
les valeurs de Vogt. Rapportée au, Mn seul, cette 
augmentation correspondrait à un y paramagnétique 
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TABLEAU à | 

Aluminium TRUR 98 123 142 169 188 209 226 253 271 290 

XA 19,20 18,67 18,58 18,26 17,99 17,86 19209 17,20 16,92 L672 
AIl-Mn RO 104 120 141 180 197,5 218 240 250 271 290 
0,473. at. % Mn 4m 20,68 20,55 20,09 19,63 19,50 19,34 19247 19,01 18,79 18,55 
Al-Mn T'°oK 98 148 138 157 181 196 210 233 250 289 
0,362 at. % Mn Xm 23,89 23,61 28:51 23,45 23,07 22,76 DATE) 22,14 21,86 21,42 
Al-Mn T °K. 94 110 134 167 190 210 234 259 270 287 
0,542 at. % Mn XM 26,15 26,00 25792 25,08 DE CAL 25,02 24,75 24,04 24,43 24,32 
AIl-Mn T °0K 94 141 160 180 194 208 22) 241 263 290 
0,684 at. % Mn Xm 28,35 27,61 27P 27 27223 26,96 20,72 26,58 26,52 26,32 26,04 
Al-Cr TOC 93 411 127 150 169 185 220 290 289 
0/095ata0710r Xm 19,09 18,91 182 18,52 il a 18,16 17,76 17,47 17,06 
Al-Cr TOR 102 120 145 A7 195 218 242 268 285 
DATES ACT Xm 19,45 19,19 18,81 18,67 18,35 18,24 1814 17,68 17252 
Al-Cr TK 94 1416 146 170 198 292 256 289 
023 "at 1Cr M 20,11 19,87 19,60 19,23 19,04 18,70 18:52 18,10 
AI-V T'oK 95 125 170 193 208 224 248 269 290 
0,078 at .V Xm 18,62 18,18 17,84 17,47 RS dE 7 17,00 16,73 16,29 
AI-V TROIKS 94 109 130 162 202 233 250 288 
DA are GAM Xm 18,08 17,88 AE 129) 16,92 16,52 16,44 15,74 


de “+ 2 650. Contrairement au cas du vanadium 
les Ay ne sont pas constants à la précision des 
mesures (l'erreur sur Ay est estimée à + 0,1), mais 
présentent autour de 130 °K un maximum étalé. 

Remarquons que la forte valeur de Ay exclut la 
simple additivité, le ya du manganèse étant d’envi- 
ron + 540 dans cet intervalle de température, et que 
le maximum en y, 7 observé [4] dans Mn & pur au 
voisinage de 125 ©K serait dans ce cas noyé dans les 
erreurs de mesures. 

Les résultats concernant le chrome sont moins nets. 
L’addition du chrome a pour effet d'augmenter la 
susceptibilité de laluminium d’une grandeur plus 


mesuré 
calculé par 
additivité 


Hrcs: 


faible que dans le cas du manganèse, mais cependant 
supérieure à Ce qu’on obtiendrait par simple addi- 
tivité (fig. 3). Les Ay, pratiquement constantes entre 
100 ©K et 200 0K augmentent légèrement avec la 
température (fig. 1). Dans la figure 3 nous avons 
porté les Ay mesurées à 290 OK. 


Une étude théorique [5] du couplage magnétique 
entre impuretés à grande distance a permis de rendre 
compte notamment du comportement des alliages à 
base de métaux nobles (Cu, Ag, Au) avec les métaux 
de transition (Mn, Cr). Les résultats décrits ci-dessus 
parlent aussi en faveur de l’existence d'interaction 
antiferromagnétiques à grande distance des atomes Mn 
et Cr à faible concentration dans l'aluminium. 

Des mesures du coefficient de Hall compléteront 
cette étude thermomagnétique. Toutefois pour l’alliage 
à 0,36 % de Mn, le coefficient de Hall coïncide à la 
précision des mesures (1 %) avec celui de l’aluminium 
pur. 
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MESURE CALORIMÉTRIQUE 
DE LA PÉRIODE DE %Y. 


Par J. ROBERT, 
Laboratoire Curie. 


_Nous avons mesuré la période de %Y avec notre 
dispositif de microcalorimétrie adiabatique à compen- 
sation automatique précédemment décrit [1]. 

L’yttrium 90 décroît par émission 8 de 2,2 MeV. 
Sa période a été mesurée par de nombreux auteurs qui 
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ont utilisé des compteurs : les principales valeurs 
obtenues sont, par ordre chronologique, les suivantes : 


Stewart, Lawson et Cork (1937) T = 60,5 +2h [2] 


% 


Bothe (1946) T=61 +1h [3] 
Chetham-Strode et Kin- : 

derman (1954) T = 64,60 + 0,43 h [4] 
Volchok et Kulp (1955) T = 64,24 + 0,30 h [5] 
Schmeling (1955) T = 64,2 +0,2 h [6] 
Wright and al. (1957) T = 64,4 +0,4 h [7] 


Les quatre derniers auteurs ont mesuré pendant 8 
à 10 périodes un certain nombre d'échantillons de #Y 
obtenu par séparation à partir de %%Sr. 

Nous avons suivi pendant une semaine dans un 
calorimètre.en or pesant 66 g une source de ŸY obtenue 
par irradiation à la pile EL 3 sous forme d'oxyde Y,0. 
L'intensité de la source atteignait, au début des 
mesures, environ 150 millicuries ; l’élévation de la 
température horaire était de l’ordre de 200 milli- 
degrés. 

La source a été mesurée de nouveau, quinze jours 
après la fin de cette étude : aucune activité anormale 
n’a été décelée. 

L'analyse des points expérimentaux par la méthode 
des moindres carrés a conduit à la valeur de la période : 


T = 64,3 + 0,4 heures. 


Ce résultat est en accord avec les précédentes déter- 
minations bien qu’obtenu par une méthode entiè- 
rement différente. 
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ÉTUDE DE LA DIFFUSION RÉSONNANTE 
DES PHOTONS PAR LE MAGNÉSIUM 


Par M. LANGEvIN et A. BussiÈère DE NERCY, 


Laboratoire de Physique Nucléaire, 
Faculté des Sciences, Orsay. 


La diffusion résonnante des photons par le magné- 
sium a été étudiée avec un faisceau de photons de 
freinage produit par le Bétatron du Laboratoire de 
Synthèse Atomique du Centre National de la Recherche 
Scientifique. MERS 

Le dispositif expérimental est semblable à celui qui 
a été déjà employé pour les expériences de diffusion 
par le niveau de 15,1 MeV de 12C [1], [2]. Le faisceau 
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de photons est défini par un collimateur primaire en 
plomb et un collimateur secondaire en fer de 80 em 
de long. La cible est placée à 4 m de la sortie du colli- 
mateur ; une chambre d’ionisation à parois épaisses de 
graphite [3] permet de mesurer l’intensité du faisceau. 
Un mur de béton de 2 m sépare le bétatron de la salle 
d'expériences et le trajet des photons depuis le colli- 
mateur jusqu’à la cible se fait dans le vide ; on a placé 
un aimant permanent destiné à éliminer les électrons 
du faisceau et un piège en béton derrière la chambre 
d'ionisation. Les photons diffusés sont détectés par un 
cristal de Nal(Tl) de 12,5 em de diamètre et de 15 cm 
de longueur protégé par une épaisseur de plomb 
de 10 cm ; les impulsions du détecteur sont analysées 
par un sélecteur d'amplitude à 50 canaux. Pour réduire 
le mouvement propre du détecteur, le sélecteur 
n’admet les impulsions que pendant le temps de pas- 
sage du faisceau (20 us) ; le rayonnement détecté en 
l'absence de cible est ainsi rendu négligeable. [angle de 
diffusion est de 1359 pour diminuer l'effet des diffusions 
électroniques prépondérantes vers l’avant et un filtre 
de 10cm d’épaisseur de plexiglass placé devant le cristal 
permet d’arrêter les électrons et d’éviter l’empilement 
des impulsions correspondant à des photons de faible 
énergie. La diffusion non résonnante subsistante peut 
être éliminée par comparaison avec une cible de numéro 
atomique voisin. 

Dans ces conditions on a pu observer la diffusion 
résonnante par le niveau de 15,1 MeV de 1?2C déjà 
étudiée par ailleurs [1], [2]. La valeur de la section 
efficace intégrée obtenue, 2,28 + 0,45 MeV.mb est en 


_ bon accord avec les résultats des autres auteurs. 


N 


150 


100 


50 


Prctr 


Le même phénomène a été observé en étudiant la 
diffusion du faisceau par une cible de magnésium. Le 
spectre d’impulsions obtenu avec un faisceau d’énergie 
maximum 12,5 MeV et une cible de magnésium de 
5,72 gJem? est représenté (fig. 1). L'énergie du pic 
observé a pu être fixée à 10,3 Æ 0,5 MeV par compa- 
raison avec les photons de 15,1 MeV de T?C et ceux 
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de 4,43 MeV d’une source de Po-Be. L'intensité du 
phénomène a été étudiée en fonction de l’énergie 
maximum des photons. Au voisinage de l’extrémité du 
spectre, le nombre de photons d'énergie Æ varie en 


V/Em E 


———— 


1 L 
A1 


l'énergie maximum du spectre [4]. Le diagramme 
obtenu en portant le produit (/r.Em)?, 17 étant linten- 


RE: 


bonne approximation comme Em étant 


10,50+0,14 MeV 
Free? 


sité de la raie de résonance, en fonction de En (fig. 2), 
permet par extrapolation linéaire de fixer l’énergie du 
niveau excité à 10,50 + 0,14 MeV. 

Le calcul de la section efficace intégrée donne 


Î os dE — 1,93 + 0,4 MeV.mb, sans tenir compte de 


lautoabsorption des photons de 10,5 MeV dans la 
cible. Si on effectue les corrections nécessaires [2], en 
prenant pour section efficace d'absorption à la réso- 
nance 68 = 674? l',[l'= 66 barn, en supposant [, #7, 


on obtient 1: 65 dE = 2,24 + 0,44 MeV .mb. Il paraît 
raisonnable, par élimination, d’attribuer ce niveau 
6 dE .2 roi — 22 + 4eV 


comparée à la largeur théorique l,(H;,) — 0,021.Æ$ 
— 24,3 eV donne comme élément de matrice de la 
transition, (ONE — 0,89 en bon accord avec les 
autres transitions M, connues [5], [6]. Ce niveau est 
probablement produit par excitation dipolaire, ce qui 
conduit à lui attribuer le spin 1+. Un niveau à 
10,46 + 0,10 MeV a d’ailleurs été déjà observé dans 
l'étude de la réaction 23Na(d, n)2Mg [7]. Le seuil de la 
réaction AMg(Y, æ«)2Ne, 9,3 MeV, étant inférieur à 
l’énergie du niveau excité de 10,5 MeV, l’observation 
de cette résonance indique que la désexcitation de ce 
niveau par émission & doit être très peu probable, On 
peut donc lui attribuer le spin isobarique T = 1. Ce 


à 2Mg. La largeur l, — 


N°“ 


niveau du 2#Mg serait donc semblable au niveau 
de 45,1 MeV de 1?C déjà étudié par la même méthode. 
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CRYOMÉTRIE DE SELS ALCALINS 
A ANIONS OXYGÉNÉS 
DANS LE FLUORURE DE SODIUM 


Par M. Ali Seyyeot et M. G. PETIT, 
Laboratoire de Physique (Enseignement), Sorbonne. 


Nous avons montré dans une précédente note [1] 
que certains oxydes, particulièrement de métaux 
positifs, étaient solubles dans le fluorure de sodium. 
On pouvait penser que la solubilité des sels alcalins 
d’anions oxygénés serait plus prononcée. Cette étude 
ayant trait alors, à la stabilité des anions dans le 
solvant choisi. 

Nous ferons ici tout d’abord, deux remarques 
intéressantes : 


a) Le clivage possible d’un anion oxygéné n’est pas 
cryométriquement décelable, lorsque nous avons le 
schéma de dissociation : 

AMOA" = M Our 07 (1) 


et lorsque l’oxyde n’est pas coupé dans le solvant 


considéré, car deux particules initiales donnent deux 


particules tronçons, que l’anion soit clivé ou non. 
Ainsi, dans le cas des métaborates, aluminates, méta- 
phosphates et leurs analogues dans la classification de 
Mendelejeff, l’étude de la structure de l’anion à l’état 
dissous doit être précédée de l’étude préalable de 
l’oxyde correspondant dans le même solvant. Le cas 
de ces anions est en effet, cryométriquement remar- 
quable. Pour l’anion métaborique nous avons en effet 
a priori un type possible de coupure de l’anion 
(BO.)- suivant le schéma général (1) non cryomé- 
triquement décelable. 

Mais on trouve expérimentalement que : 
19 1(BO;)- donne une particule dans FNa :; 
2° 1 B,0; est clivé en deux particules dans FNa [7]. 

Ceci élimine donc la coupure de type (1) qui 
donnerait trois particules tronçons pour deux parti- 
cules initiales. 

Done la non coupure expérimentale révèle que la 
coupure possible de (BO,)- en (BO)}+ + O= ne peut 
se réaliser (éventuellement) qu’à des dilutions trop 
élevées pour pouvoir être suivies par des eryométries 
où le métaborate est le corps dissous. 

(Mais si loxyde B,0, avait donné une particule 


no 
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dans un certain solvant le problème eût été alors 


_ indéterminé.) 


Ainsi l'étude de la coupure d’un oxyde d’un élé- 
ment trivalent ou pentavalent est préliminaire à 
l'étude de la coupure de l’anion correspondant. 

Par contre les orthophosphates et homologues, les 
sulfates et homologues donnent des résultats Cr yo- 
métriques directement interprétables de ce point de 
vue, car n’obéissant pas au schéma n° 4. 


(0) 50 100 150 - 200 250 


F1G. 1. — v représente, à la limite de dilution infinie, le 
nombre de particules étrangères au solvant donné par le 
corps dissous. 


b) L’idéalité, traduite pratiquement par l’horizon- 
talité de la courbe cryométrique dans la représenta- 


tion de Raoult aux hautes dilutions (aux erreurs 


d'expérience près et au voisinage de la limite) pour 
l’anion orthophosphorique, implique que pour des 
anions tétraédriques (et «a fortiori pour les plus 
petits) l’entropie correspondante de forme et dimen- 
sion n'intervient pas et que seule donc la théorie 
d’Arrhenius rend compte de la « réalité » en milieu 
sel fondu électrolyte fort [6]. 

D'ailleurs dans ce cas et conformément aux résul- 
tats expérimentaux de Zarzycki [8] on peut avec 
Doucet [9] expliquer ce fait en se basant sur la 
structure très lacunaire de FNa, dont les lacunes 
logeraient les tétraèdres anioniques sans perturbation 
notable de la stéréostructure du solvant. Ceci est 
donc une nouvelle preuve théorique séduisante du 
bien fondé de l’hypothèse fondamentale que nous 
avançons plus haut. 

Ceci étant posé, nous avons obtenu les résultats 
suivants : 

De nombreux sels sont insolubles ou instables dans 


le fluorure de sodium, tels les aluminates insolubles, 
les arséniates, vanadates, stannates, tellurates, ins- 
tables. 


Les métaborates révèlent un non clivage de 
l’anion (BO;)-, tout au moins dans les limites de 
dilution accessibles à l'expérience, et compte tenu de 
la remarque générale précédente. 

Les tétraborates (B,O,)Na, donnent une coupure 
en quatre particules à la limite, donc vraisembla- 


blement 3(BO;)- + (BO)*. 
Nous avons étudié alors les phosphates (ortho, 
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pyro, méta) et les anions complexes dont les atomes 
centraux appartiennent à la 6€ colonne. 


Phosphates, — L’anion orthophosphorique (PO,)= 
n’est pas clivé à 995 oC. Nous trouvons donc un 
résultat identique à celui trouvé par Harrison à plus 
basse température [3]. L’idéalité de la solution a 
motivé notre seconde remarque faite plus haut. 

. Les pyrophosphates donnent deux particules à la 
limite : on peut donc penser à un clivage en (PO,)= 
et (PO). Le type de clivage possible de (PO;)- 
n'étant pas démontrable cryométriquement par suite 
de la première remarque générale et de la très grande 
difficulté d’étude d’une solution d’anhydride phos- 
phorique. 


Dérivés de la 6e colonne. — Les sulfates et chro- 
mates dans FNa ont été étudiés par Mlle Bour- 
lange [4] qui a trouvé une non coupure pour 
les (S0,)=, une légère coupure pour les chromates. 
Les molybdates et tungstates ont été étudiés par 
Harrison [3] et Mergault [5] dans CINa et dans la 
cryolithe. Comme ces auteurs nous avons trouvé dans 
le fluorure de sodium une légère coupure de l’anion 
complexe avec une chute lente vers une particule 
(voir fig. n° 1) et pour les tungstates une courbe 
pratiquement horizontale pour une particule tra- 
duisant une non coupure 


Silicates. — Le métasilicate de sodium donne une 
courbe légèrement inclinée et tendant vers une parti- 
cule à la limite. Il y a là une divergence avec les 
résultats obtenus par G. Petit [6] dans la cryolithe 
où cet auteur trouvait un clivage en silice et ion 
oxygène. Sans doute faut-il voir là une manifes- 
tation de l’affinité des particules aluminiques de la 
cryolithe pour l’ion oxygène. 

En résumé, le fluorure de sodium est un solvant 
relativement médiocre, beaucoup plus voisin à ce 
point de vue du chlorure de sodium que de la cryo- 
lithe malgré l'identité des températures de fusion des 
fluorures simples et complexes. 

Ceci révèle que les propriétés solvantes d’un sel, 
sont, toutes proportions gardées, plus fonction de sa 
structure que de sa température de fusion. 
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ÉTUDE COMPARATIVE 
DES SPECTRES D'ABSORPTION INFRAROUGE 
ET DE RÉSONANCE MAGNÉTIQUE NUCLÉAIRE 
(R. M. N.) 
DE SOLUTIONS D'EAU 
DANS DES CORPS ORGANIQUES 


Par G. MAVEL, 


Laboratoire de Spectroscopie hertzienne à la Sorbonne. 


Introduction. — Nous avons complété certaines de 
nos mesures et repris celles d’autres auteurs [1, 2], 
sur la R. M. N. du proton de l’eau dans des corps orga- 
niques ; nous nous proposons de comparer ici ces 
résultats à ceux de É. Greinacher, W. Lüttke, 
R. Mecke [3], relatifs à l’absorption infra-rouge. Les 
solvants étudiés sont : CCl,, CS,, benzène, nitro- 
méthane, acétonitrile, acétone, dioxane, éther, tétra- 
hydrofuranne (T. H. F.), pyridine, cités dans l’ordre 
de AV — Vyaz — Volution Croissant [3]. Nous avons 
ajouté à cette liste : chloroforme, morpholine, pipé- 
ridine, « et 8 picolines, corps pour lesquels nous avons 
également mesuré le déplacement infra-rouge. Nous 
avons trouvé dans un travail assez ancien de 
W. Gordy [4] des éléments de comparaison intéressants 
avec D,0. 

Un parallèle de ce genre a déjà été établi, par 
W. G. Schneider, H. J. Bernstein, J. A. Pople [5] sur 
divers composés hydrogénés, mais avec des données 
infrarouges d’une grande diversité d’origine. Nous 
pensons que l’homogénéité de nos données donnera 
plus d'unité à notre travail. Nous nous proposons de 
comparer Av — Vyaz — Vol. /nesuré sur les bandes 
fondamentales vom (moyenne de YS et Vr) et 
AS — Dyaz — do. (0.01. étant, lorsque le solvant n’agit 
pas par un groupement polaire contenant un H, la 
position limite du déplacement à pour une dilution 
infinie de l’eau). Comme précédemment [1], nous uti- 
lisons Ôgaz — 4,37 en tenant compte de la température 
à laquelle ont été faites les mesures, 22 °C 2 [5]. 


Partie expérimentale. —— Pour les corps se dissolvant 
dans l’eau en toutes proportions, nous avons établi les 
courbes Ô(x) ; nous donnons (fig. 1), celles qui corres- 
pondent à l’acétonitrile, Ia morpholine et la pipéridine 
(pour laqueile nous avons complété des mesures de M. 
et R. Freymann) ; on notera que pour ces deux derniers 
corps, qui contiennent des groupements NH, il y a 
triple échange [6] entre les protons de H,0, du com- 
plexe et du NH. Dans ce cas, la courbe aboutit à la 
position du NH pur, et non à Ô.a, que l’on obtient 
par un calcul [6]. 

Lorsque la solubilité de l’eau est faible, nous avons 
pris une solution saturée (CS,, benzène) ou à 0,5 % 
d’eau, en volume. 

Les mesures, en R. M. N., ont été faites avec un 
appareil Turb-Taüber (25 MHz; échantillon tournant ; 
bandes latérales). La figure 1 donne la valeur de à. 
pour les corps étudiés ; pour la morpholine et la pipé- 
ridine, le calcul fournit AS = 3,6 + 1,3 et 3,3 + 1,2. 
On remarquera la position aberrante du sulfure de car- 
bone ; on peut la rapprocher de celle du benzène, dont 
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le désaccord est moins net, et qui est cependant cohé- 
rente avec les mesures de Cohen et Reid [24]. 

: re 

Les mesures complémentaires, dans l’infrarouge, ont 
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été effectuées (*) avec un appareil Perkin-Elmer 21, à 
prisme de NaCI, identique à celui de Greinacher, 
Lüttke et Mecke [3]. 


Discussion. — La figure 2 donne, dans les axes Av 
en cm1 et AS en parties par million (p. p. m.), les. 
points expérimentaux : on constate que la relation est 
sensiblement linéaire entre les résultats des deux 
méthodes. On constate également que parmi les corps 
s’écartant de la droite, l’un des plus significatifs est le 
benzène ; ce désaccord est difficile à interpréter. On 
notera que l’infra-rouge ne donne aucune différence 
entre les trois corps contenant Of (dioxane, éther, 
T. H. F.), alors que la R. M. N. les sépare nettement ; 
il faut noter, à ce propos la faible dispersion du spec- 
trographe infra-rouge dans la région 3 000-5 000 em-1. 

Le tableau I fournit les valeurs de : AS en p. p. m. ; 
Avon (fréquence de valence), ASox (fréquence de défor- 
mation), pour H,0, en cm-1; Avon (valence) pour 
D,0 [4]; P. I, potentiel d’ionisation du solvant, 
donnant une évaluation de sa capacité de donneur [7] ; 
moment dipolaire de la molécule de solvant et moment 
de son orbitale associante («lone pair dipole ») [8]. 

Dans le cas du chloroforme comme donneur [9], la 
comparaison des résultats de R. M. N. avec le P. I. et 
avec le moment de l’orbitale associante avait été satis- 
faisante, dans l’ensemble. Nous voyons que, pour les 
corps susceptibles de s'associer à l’eau, en agissant 
comme accepteur, il ÿ a une corrélation qualitative, les 
picolines mises à part. Mais, vis-à-vis de la pipéridine, 
de la morpholine et de l’eau elle-même, Peau agit, au 
moins partiellement, comme accepteur et c’est l’affi- 
nité électronique du solvant qu’il faudrait plutôt consi- 
dérer, Dans les autres cas, le P. I. décroit lorsque AD 
croit. Le moment dipolaire de l’orbitale associante 
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: Vgaz - Vo ;: révèle des positions anormales, en particulier celles de 
‘ cm” Ê Picoline / l’acétone et surtout de l’acétonitrile, ce qui avait déjà 
| 300 CO & Picoline signalé en infra-rouge [3]. 
, Pyri dine Un dernier point intéressant à remarquer est Pécart 
Pipéridine / plus ou moins grand par rapport à la désassociation 
NP MU CAE — parfaite, réalisée pour l’eau à l’état vapeur [5]. Cet 
| / Eau liquide 
À / 
; / D_1 ! 
6 É / 
/ Dr pl 
/ 
; - 7 Morpholine ! 
PURES" ASE A —©------—- ! ba 
, / 0,5 Acétonitrile H,0 liquide 
74 es © Acétone (©) 
! 
k AREA Fthes: . Nitrométhane © 
© © OTHF O4 ! Fer 
/ 1 
/ | Éther 
Acétonitrile © Acétone [Chloroforme ® K. B.M. 
- f ee ] Morpholi 
x PS RNE 
| 100 “ rire Pipéridine 7 L. @ Picoline 
J Honzone ve © Witrométhane de © 7 0 CG) 
l © £ ’ ! CCL, TH © f Picoline 
@ HCL; LUE 
/CCL4 
/ 
i ( 1 ê 4 LE 
44 vapeur < Ô gaz = OÜcosy (p.p.m.) ‘ 
| rc 
TABLEAU I 
; A AVorx AVop AVox Par MomEenT pre. MOMENT DIr, 
L- : p. p. m. cm1 Cm Cie eV TOTAL ORB. ASSOC. 
| | (3) (4) (3) (7) (8) 
| H,0 vap. (0) (0) (0) (0) == = Æ 
DO CC 0,622a 45 es 0 11,5 0 je 
4 CS» æ— 0,3?) A9 = 0 11,1 0 D 
CHU 0,75?a DIE == — 11,4 1,02 _— 
: 0,29 70 — 0 9,2 0 — 
Benzène 0,122 
Nitrométh. 1,88 78 — — 28 — 3,25 — 
| Acétonitrile 1,22 115 138 ESA ae 3,9 3,74 
; 2,02 127 — — 36 O7) 2,88 3,03 
Acétone | 2 62b 
Dioxane 1298 153 —— — 43 9,9 0,4 2,78 
Éther A1 153 — — 43 9,5 1,22 2,78 
TH. F. 2,89 153 = en = 1,68 2,78 
Morpholine 3,6 205 * — — (=) = L 
Pipéridine 3,8 270 315 = (—) — — 
B picoline 4,30 305 — — 9,71 —— 3,90 
æ picoline 4,45 305%* 27) — 9,66 — 5 20 
H1,0 liquide 4,37 276 a eu (12,6) 8 
Pyridine 4,50 292 269 — 55 9,2 3.05 


* Nos mesures infra-rouge. 


- 
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écart est dû à l'effet décrit par Kirkwood [10], Bauer 
et Magat [11]; il est possible, comme dans linfra- 
rouge d'utiliser un diagramme AG fonction de 
(D —1)/(2D + 1), nous le donnons (fig. 3). On y 
constate que l’acétonitrile se conduit bien comme un 
solvant inerte, le nitrométhane étant plus actif. Les 
autres corps Ont un comportement normal, compte 
tenu des propriétés que nous leur connaissons. 


(*) Par les soins de Mme Dumain, du Laboratoire Central 
des Poudres, que nous tenons à remercier. 


Lettre reçue le 11 juillet 1959. 


Note ajoutée à la correction. — D’après Pople (J. Chem- 
Phys., 1958, 28, 601) le déplacement AS est proportionnel 
D —1 
2D Eu 
est proportionnel au champ lui-même. Le diagramme mo- 
difié en conséquence ne montre pas de différence essen- 
avec celui donné ici (fig. 3). 


au carré du champ créé par le solvant, alors que 
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VALEUR APPROCHÉE 
DE QUELQUES MOMENTS 
QUADRIPOLAIRES ÉLECTRIQUES 
DE NOYAUX DE MASSE A > 218 


par R. Foucxer 


Laboratoire de Physique Nucléaire 
Orsay (Seine et Oise). 


te 
2 existence de moments quadrupolaires électriques 
liés à une déformation permanente des noyaux de 
N et Z éloignés des nombres magiques, même si leur 
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spin est nul, est bien établie, mais leur relation avec 
les moments d'inertie correspondants reste un pro- 
blème des théories du noyau et leurs déterminations 
expérimentales dans la région des masses atomiques 
A > 210 par spectrométrie optique ou par excitation 
coulombienne sont peu nombreuses et peu précises 
[1] [2] [3]. Les atténuations observées dans les corré- 
lations angulaires &-y (noyaux pair-pair) ont été 
expliquées par le couplage de ces moments avec les 
gradientsde champ électriques cristallins [4] et per- 
mettent de prévoir la valeur relative des vies moyennes 
de quelques niveaux de rotation. Les vies les plus 
longues ont été mesurées directement par la méthode 
des coïncidences différées (221Ra, 226Ra) [5][6]. On en 
déduit les moments quadrupolaires à partir de la for- 
mule de Bohr et Mottelson [7], si on connait les coeffi- 
cients de conversion interne [8]. La vie moyenne du 
niveau 2+ de 241 MeV du ?* Rn, noyau situé dans la 
zone de transition vibration-rotation [9], a été mesurée 
[10] et le moment Q, qui en est déduit se place bien 
sur la courbe extrapolée Q, (A) [1] ; ceci amène à 
penser que la formule de Bohr et Mottelson est uti- 
lisable dans la région À > 218 pour prévoir une 
valeur approximative des Q, et vies moyennes + des 
niveaux à partir des corrélations angulaires &-y des 
220Rn 2*2Rn, 224Ra, 226Ra, 2#4U [11] avec comme seules 
hypothèses supplémentaires que le noyau de recul 
est arrêté lorsqu'il émet le rayonnement y (temps 
d'arrêt — 10-12 sec) et que les gradients de champ 
électriques q des chlorures et nitrates de ces éléments 
sont identiques. 

Dans ces conditions les fréquences d'interaction 
quadrupolaire sont toutes de la forme «, — & qQin 
et les vies moyennes de la forme +, — b,/Q8, ; nous 
obtenons donc les valeurs relatives des Q, et 7, à 
partir des produits ©, +, tirés des facteurs d’atténua- 
tion mesurés et si nous choisissons le 226Ra comme 
noyau de référence (ot & 1, tr — 9.10-19 sec) ces 
valeurs sont : 


N Gÿ, Do Tn(10— 10 sec) Qon barn 
220Rn 0,2 RUES ÉRUS 
222Rn 0,4 6,3 LR 
24Ra 1,3 12 6,5 
226Ra 1 9 7 5 
At 0, 5,6 8,6 


La similitude des Qm et 7, obtenus ainsi et par 
mesure directe ds vies moyennes pour les :2hRn 
**‘Ra, ?6Ra nous amène à estimer que Qm et +, 
(“2Rn et #4U) sont correctes. Il est certainement 
très approximatif de supposer constant le produit 
E,X Q6* A°/3, où E, est l’énergie du niveau 2 du noyau 
pair-pair de masse À; cependant nous l’admettrons 
pour a noyaux voisins du ?*%Rn dans la mesure où 
nous retrouvons ainsi une valeur voisine 
pour < Q, > (???Rn); alors < Q, (PP Rn) 2 ar 
et Q(Ra) — 6,6 barns. Enfin à partir de QT h) et 
@o(**Th) voisins de 7,5 et 8 barns et Q(AÀ > 238) 
inférieur à 14 barns, nous pouvons estimer la vie 


< 


Ts 1 
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moyenne des premiers niveaux correspondant. En 
résumé nous avons le tableau suivant : 


N E; t(10— 10 sec) Qo(barn) 
tn 325 12,2 RSR 
220Rn 241 4,5 OR 
22Rn 187 6,3 she 
222Ra 111 7,8 6,6 
224Ra 84 12 6,5 
CCR 68 9 GE 
226Th 72 6,8 7,5 
228Th 57 3,5 8 
232Th 50 8,5 
2341] LS 5,6 8,6 
2361] 45 10,3 

10-2258 > 40 > 1 < 1% 


Les valeurs de Q, et + indiquées ne sauraient être 
qu’approchées et leur détermination par la mesure 
directe des vies moyennes et l’étude des corrélations 
angulaires &-Y est en cours. Elles permettent néanmoins 
de fair: déjà plusieurs remarques sur les vies moyennes 
des premiers niveaux 2+ des noyaux pair-pair de 
masse À > 218 et sur leurs moments quadrupo- 
laires : 

1) Les vies moyennes expérimentales des premiers 
niveaux 2+ de cette région passent, lorsque À croit, 
par un maximum au ?**Ra puis restent supérieures 
à 10—19 sec ou de cet ordre. La valeur estimée pour 
T(?26Th) ne tient pas/compte du fait que dans ce cas 
E,(22Th) > E,(??SRa) alors que l'énergie E, de 
l’isobare de Z le plus élevé est normalement la plus 
faible. 

2) L’addition de deux neutrons fait varier le mo- 
ment quadrupolaire moyen des Radons de 1 barn 
environ ; celui des radium varie moins. Par contre les 
moments varient approximativement de 2,5 à 3 barns 
lorsqu'on ajoute 2 protons au *#%Rn et au ?#Rn (Dans 
la région des terres rares cette variation brusque de 
@ à lieu pour N = 88). 

Lettre reçue le 24 juin 1959. 
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POUVOIR SÉPARATEUR ET EFFICACITÉ 
D'UN SPECTROMÈTRE B A ORBITES SPIRALES 


Par Mme Adèle BuRDET, 
Institut de Physique Nucléaire, Université de Lyon. 


Le principe de ce spectromètre, imaginé par 
G. Miyamoto [1], est basé sur les propriétés des tra- 
jectoires des particules chargées dans un champ magné- 
tique de révolution. 


Caractéristiques générales d’un champ magnétique 
de révolution. — Nous représenterons la variation de 
l'intensité d’un champ magnétique de révolution à 
l’aide des coordonnées cylindriques r, ®, z ; où l’axe 
des z est l’axe de symétrie du champ, et le plan z — 0 
le plan équatorial de l’entrefer. Nous bornerons cette 
étude à ce plan. La variation de l’intensité en fonction 
de la distance à l’axe présente deux particularités essen- 
tielles : au voisinage de r — 0, la courbe A(r) est tan- 
gente à une droite horizontale, le champ étant continu 
en chaque point ; pour r très grand devant le rayon des 
pièces polaires l’électro-aimant se comporte comme un 
doublet magnétique, H(r) est en 7 fr. 

On en déduit la courbe rA(r) sur la figure 1. 


Par suite la courbe rA = | r.H(r).dr sera assi- 


milable à une parabole au voisinage de l’origine, puis le 
maximum de 74 indique un maximum de pente 
pour rA et, enfin, à l'infini rA se comporte 
comme — kfr, hyperbole équilatère asymptote à une 


re Courbes trHetrrA—= [ r.Har. 


droite horizontale. Il y a donc pour 7 A un maximum 
de concavité, ceci quel que soit la topographie de 
H entre 0 et l'infini (fig. 1). 


Équations du mouvement et forme des trajectoires. 
_— Üne particule chargée est soumise à la force de 
Lorentz F — e (w/\ H). Le potentiel vecteur À se 
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réduit à sa composante Ao d’après le choix des axes. 
L’équation des trajectoires décrites par ces particules 
est donnée par : 


r rA [pH 
PPI ES 1 ROIS PU ES dr 
o rVr?— (rA/pH)? 


ee 


pour une particule issue du centre [2]. La forme des 
trajectoires dépend des racines du dénominateur de la 
fonction à intégrer. S'il admet une racine double la 
trajectoire correspondante est asymptote à l'orbite 
stationnaire instable relative à l’énergie considérée. 
Par contre pour une énergie inférieure la courbe décrite 
sera bornée, alors que la particule s’éloignera à Pinfini 
si son énergie est plus grande. Il y a accumulation des 
particules sur l’orbite stationnaire (cercle de conver- 
gence). 

Ainsi le champ magnétique va séparer les particules 
suivant leur énergie. Il pourra donc permettre de cons- 
truire un spectromètre, et 1l suflira de placer un détec- 
teur sur le cercle de convergence pour recueillir les 
particules de 0/7 donné soit 6 — B. 


Pouvoir séparateur, — Représentons le détecteur par 
un segment de 1 em de longueur : une particule est 
détectée si sa trajectoire coupe ce segment. On peut 
voir sur la figure 2 que certaines des particules de po 


- 


D .. 


F1G. 2. — Définition de l’angle utile. 


différents de B atteignent le détecteur, il est donc 
nécessaire de faire une étude du pouvoir séparateur. 
Pour cela définissons l’angle utile du spectromêtre pour 
un bp donné : c’est l’angle Aw(0/1) à l’intérieur duquel 
doit se trouver la vitesse initiale de la particule consi- 
dérée pour qu’elle soit détectée (fig. 2). 

L'étude de la variation de Ao en fonction de 9H va 
nous permettre de déterminer le pouvoir de résolution, 
défini par: 0(0/7)1B, (d(oH) étant la largeur de la raie). 
Nous nous donnons la probabilité pour qu’une parti- 
cule d'énergie correspondant à LA soit émise : n(pH). 
Pour qu’elle soit détectée il faut que Ao(o/1) ne soit 
pas nul, le nombre de particules détectées sera done : 


N = f n(oHl).Ao(oH).d(pH). 


En supposant une raie gaussienne théorique à 
l’entrée du spectromètre, on peut en déduire le profil 
de la raie qu’il restitue. La forme de Aw(oH) montrée 


sur la figure 3 laisse prévoir une raie fortement dissy- 
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métrique présentant un pied élevé du côté des pH plus 
orand. que B (fig. 4). Il faut chercher une position du 
détecteur en le décalant vers l’extérieur du spectro- 
mètre ou au contraire vers le centre. 


js 


2H 


B20 006 oH 


Fic. 3. — Variation de A autour de B. 


100 


B-20 B ch 


Frc. 4 — Profil de raie. 


le-L ls 


F1G. 5. — Variation de + en fonction de r pour eH > B. 
D’après la variation de + en fonction de r, la difté- 


rence ®7 — @,1 est plus faible quand on se rapproche 
du centre (fig. 5), il faudra donc mettre le détecteur le 
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plus au centre possible, c’est-à-dire de façon que le 
bord externe coïncide avec le cerele de convergence. 

Avec une topographie de /{(r) donnant une orbite 
stationnaire de rayon 15 em pour une particule d’éner- 
gle correspondant à po — 2000 gauss.cm, on a 
obtenu un pouvoir séparateur de 0,6 °/,. 


Généralisation. — Nous allons étendre ces résultats 
à une source circulaire utilisée en réalité. On peut la 
considérer comme constituée par une infinité de sources 
élémentaires ponctuelles, situées à la distance 7, du 
centre. Pour de telles sources, pour une énergie donnée 
il y a plusieurs types de trajectoire suivant l’angle de 
tr « de la particule [3]. Pour & — À la trajectoire est 
la spirale, si œ« < À la trajectoire est à branche infinie, 
et six > À On a affaire à une trochoïde (fig. 6). 


Fic. 6. — Trajectoires pour « = X, « < À et &« > À. 


Soit u(x),x l'arc sur lequel doit se trouver la source 
élémentaire pour que la particule d'énergie déterminée 
atteigne le détecteur. Pour un pH donné on peut 
tracer la courbe u(x), la surface délimitée par cette 


courbe représente le nombre de particules recueillies ; 


pour la source ponctuelle ce nombre était représenté 
par Av, ainsi de u(x) dx généralise la notion d’angle 
utile. 

Notion d’efficacité. — Nous appellerons efficacité du 


spectromètre relative à la quantité de mouvement p/7 
le rapport : 


nombre de particules d'énergie pA recueillies 
-_ nombre de particules de cette énergie émises” 


En 


Ce rapport est aussi égal au rapport de la surface 
délimitée par u(x) à celle du rectangle 27 X tr, chaque 
source ponctuelle émettant dans x, et il ÿ a 27 sources 
ponctuelles. 

Pour une source circulaire de diamètre 5/10 mm 
nous avons trouvé une efficacité de 40 %, tandis que 
l’on peut atteindre 80 % pour une source de 2/10 mm 
de diamètre. Pour la première le pouvoir séparateur 
est devenu 0,65 %, soit légèrement plus faible que 
pour une source ponctuelle. 


Influence du centrage de la source. — On a également 
montré que si la source s'éloigne du centre Pefficacité 
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diminue très rapidement et varie suivant les positions 
relatives de la source et du détecteur. 

Le spectromètre & à orbites spirales peut atteindre 
une efficacité de 80 % et un pouvoir séparateur 
de 0,6 %. Il faut que la source soit parfaitement 
centrée, ce qui impose un champ rigoureusement de 
révolution et une détermination précise de la position 
de l’axe de symétrie. Cet appareil de précision, sera 
utile surtout pour l’étude de sources dont l’activité 
Aie est très faible (grand pouvoir de transmis- 
sion). 


Lettre reçue le 17 juillet 1959. 
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CALCUL DES VIBRATIONS FONDAMENTALES 
DES MOLÉCULES CO, ET CS, 
COMPTE TENU DE L'ANHARMONICITÉ 


Par Itaru GAmo, 


Laboratoire de Chimie, 
Université de Waseda, Tokio, Japon. 


Nous avons calculé les vibrations normales de CO, 
et de CS, à l’aide de la fonction potentielle la plus 
générale, en tenant compte de l’anharmonicité des 
vibrations. Les fréquences normales (fréquences 
d'ordre zéro) adoptées dans ce calcul sont celles esti- 
mées par Dennison [1] pour CO, et par Wentink [2] 
pour CS,. Le tableau suivant donne leurs nombres. 
d'onde en em-!, avec les fréquences fondamentales 
correspondantes observées. 


FRÉQUENCES FRÉQUENCES TYPE DE 
NORMALES FONDAMENTALES VIBRATION 
(w) (o) 
+ 
1 351,2 1 388,4 se 
CO; 672,2 667,3 TE, 
2 396,4 2 349,4 D 
674,4 658,0 El 
CS 398,6 396,8 IT: 
\ 1 551,9 1 532,5 D 


Les constantes de force ont été calculées à partir de 
la fonction potentielle suivante : 


2V = k,{ (AR)? + (AR)? } + k5(A0)? + Æ:2(AR3) (AR:) 
où AR,, AR, et A6 sont les variations des dis- 
tances C-O ou C-S et de l'angle dC-OPVoU 
S-C-S à partir de la position d'équilibre. 


On a admis que les angles ont pour valeur 1802 pour 
les deux molécules. 
On a les relations suivantes : 


ta ki + Kio 
QE 
DA 

ARS 

ONE 

1 2 
CE té VS 
E = F = ei 12) 


où À = 472 c2 w2/N ; cest la vitesse de la lumière, 
N le nombre d’Avogadro, m; la masse atomique du 
carbone, m, celle de l'oxygène ou du soufre. 

Pour ces valeurs, on a adopté respectivement 
12,003844,  16,000000 et  31,982274. 
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c — 2,997902 x 101 erm/sec,  N = 6,02544 X 10% 
R est la distance d'équilibre de la liaison C-O ou C-S. 
On a calculé les valeurs suivantes des constantes de 
force (en dynes par centimètre . 10°). 


k: ki2 ks [R? 
CO» 15,980 2,441 0,5807 
CS, 78279  0,6615  0,2352 


Lettre reçue le 16 juillet 1959. 
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REVUE DES LIVRES 


Pauzr (W.), Théorie de la relativité. (1 vol. relié 
45,5 X 25 cm, de 241 p., Pergamon Press, Londres et 
N. Y., 1958. Prix 35 shillings) (en anglais). 

En 1921, le jeune Pauli, âgé de 21 ans donnait un livre 
aujourd’hui classique : Relatiouäts theorie. Trente-cinq 
ans plus tard, paraît une traduction anglaise de cet ouvrage. 
Mais, à la fin, 26 pages de notes donnent un sommaire des 
progrès effectués. Signalons deux notes importantes sur la 
différentiation covariante dans un espace à connexion 
affine, et sur le problème cosmologique. Enfin dans la der- 
nière note (9 pages à elle seule), le vieux et illustre phy- 
sicien (il devait disparaître peu après) donne son opinion 
sur les tentatives de théorie unitaire. 11 examine la théorie 
d’Einstein et Kaufmann avec gix et lu dissymétriques,et la 
théorie de Kaluza dans l’espace à cinq dimensions. Il pré- 
sente de sévères objections et ne cache pas son scepticisme. 
Même l’unification de la gravitation et de l’électromagné- 
tisme est un échec, pour ne pas parler de l’unification de la 
Relativité générale et des quanta, irréalisable du point de 
vue Bohrien, qui est celui de M. Pauli. 

J. WINTER. 


KrisaNAN (R.S$.), Progrès en physique cristalline. (Progress 
in erystal physies.) 1 vol., 19 X 25 cartonné ; 198 pages. 
Central Art Press, « Aeton Lodge », Chetput, Madras 31. 
L’Académie Indienne des Sciences a été depuis plus de 

quinze ans à l’avant-garde des travaux de physique du 
solide et tout particulièrement de l’état cristallin. Les très 
nombreux travaux sur le spectre Raman, les propriétés 
élastiques, photoélastiques, thermiques, magnéo-optiques 
et diélectriques des cristaux, font que ces laboratoires 
occupent une position unique dans ces domaines de 
recherches. 

Le Docteur R. $S. Krishnan aidé des nombreux spécia- 
listes de l’Académie Indienne des Sciences (Vedam, Sun- 
dara Rao, Ramasesham, Narasunham, etc...) présente ici 
une monographie d’un intérêt considérable pour les cher- 
cheurs qui abordent les problèmes de physique cristalline. 

Les septs chapitres de cette monographie sont : 40 Dis- 
location thermique. 2° Conductibilité thermique. 3° Pro- 
priétés élastiques. 4° Propriétés photoélastiques. 5° Com- 
portement thermo-optique. 6° Effet Faraday. 7° Compor- 
tement diélectrique. 


Chaque chapitre contient une brève étude théorique, 
une description des méthodes expérimentales de mesures, 
les résultats obtenus et une discussion. \ 

Des tableaux très complets de constantes numériques | 
suivent chaque chapitre, ainsi qu’une liste très abondante 
de références bibliographiques, à jour jusqu’à décembre 
1957. 

Les cristaux métalliques et les matériaux polycristallins 
sont exclus de cette monographie. 

La présentation du livre est agréable et pratique ; 
l’ouvrage est bien imprimé et très aéré. 

S. PAUTHIER. 


Doucer (Y.), Techniques modernes et applications de la 
cryométrie. (1 vol. relié, 17 X 11 cm, 225 pages. Mono- 
graphies Dunod, Paris, 1959.) 

L’étude des lois de Raoult régissant l’abaissement de 
température de congélation des solutions organiques ou 
non fait partie depuis longtemps des programmes du bacca- 
lauréat. Malheureusement, les ouvrages élémentaires qui 
traitent de la cryométrie n’indiquent comme application. 
que la détermination « approchée » des masses imolaires, 
ce qui laisse penser qu’il s’agit là d’une méthode imprécise 
aux applications limitées. 

Or, il n’en est rien. Depuis Raoult (c’est-à-dire 1880-1900) 
les techniques se sont notablement améliorées. En parti- 
culier, la méthode dite d'équilibre dans laquelle on cherche 
la concentration d’une solution, dont le point de fusion- 
congélation est stable et bien déterminé, donne des résul- 
tats d’une grande précision. 

Aussi les applications actuelles sont-elles nombreuses et 
variées. L’auteur en décrit une quinzaine après avoir soi- 
gneusement passé en revue tous les perfectionnements qui 
ont été apportés récemment aux appareils utilisés, aux 
méthodes physiques de mesure de la concentration, à l’équi- 
libre adiabatique et à la technique de la prise de tempé- 
rature. L 
. Ouvrage qui sera extrêmement utile à tous ceux : Phy- 
siciens, Chimistes, Biologistes, Médecins, Ingénieurs agro- 
nomes, qui utilisent la cryométrie comme technique de 
recherches. 


P. Rouanrp. 


